A prospective cohort study to assess the micro-epidemiology of Plasmodium falciparum clinical malaria in Ilha Josina Machel (Manhiça, Mozambique)

Background After the decrease in clinical malaria incidence observed in Mozambique until 2009, a steady resurgence of cases per year has been reported nationally, reaching alarming levels in 2014. However, little is known about the clinical profile of the cases presented, or the possible epidemiological factors contributing to the resurgence of cases. Methods An analysis of surveillance data collected between July 2003 and June 2013 in the high malaria-transmission area of Ilha Josina Machel (Southern Mozambique) through a paediatric outpatient morbidity surveillance system was conducted to calculate hospital-based clinical malaria rates, slide-positivity rates, and minimum community-based incidence rates (MCBIRs) and incidence rate ratios per malaria season in children younger than 15 years of age. Clinical malaria was defined as a fever ≥37.5 °C or a reported fever in the previous 24 h with a positive blood smear. Yearly mean age, geometric mean parasitaemia (GMP) and mean packed cell volume (PCV) were also described for all clinical malaria cases and compared between seasons using DID analysis or ANOVA tests. Results During the study period, the percentage of outpatient visits presenting with confirmed clinical malaria decreased from 51 % in the 2003–2004 season to 23 % in 2008–2009, followed by an increase back to 51 % in 2012–2013. The yearly mean age of cases significantly increased from 2.9 (95 % CI 2.8–3.0) in 2003–2004 to 5.7 (95 % CI 5.6–5.7) in 2012–2013, compared to non-malaria cases. An increase in mean PCV levels was also observed (p < 0.001), as well as in GMPs: from 5778 parasites/µL in 2002–2003 to 17,316 parasites/µL in 2012–2013 (p < 0.001) mainly driven by an increase in GMP in children older than 1 year of age. MCBIRs in infants decreased by 70 % (RR = 0.3, p < 0.001) between 2003–2004 and 2012–2013. Incidence diminished by a third among children 1- to 4-years between 2003 and 2007, although such drop was unsustained as observed in 2012–2013 (RR = 1.0, 95 % CI 0.9–1.0). Finally, the incidence among children 5–14 years was 3.8 (95 % CI 3.4–4.3) times higher in 2012–2013 compared to 2003–2004. Conclusion Since 2003, Ilha Josina Machel observed a significant reduction of clinical malaria cases which was followed by an upsurge, following the national trend. A shift in the age distribution towards older children was observed, indicating that the changes in the transmission intensity patterns resulted in a slower acquisition of the naturally acquired immunity to malaria in children.


Background
As part of the Millennium Development Goals (MDGs) established by the United Nations in the year 2000, all malaria-endemic countries concentrated their efforts on reducing malaria incidence globally by 2015. In the last decade, 55 of the 106 countries with ongoing malaria transmission reached the target set in Objective 6 to reduce their malaria burden by 75 %. However, these countries contributed only to 6 % of the global estimated cases, suggesting that countries with the most affected populations, most of them from sub-Saharan Africa, have yet to reach this target [1].
Mozambique is part of the ten countries with the highest malaria endemicity in the world. Transmission intensity varies widely within the country, with high transmission intensity occurring in the north and a decreasing pattern observed as provinces are closer to the south of the country [2][3][4]. However, high transmission hotspots can still be found within wider low transmission areas. Despite the increasing funding available for malaria control, Mozambique has experienced an increase in malaria incidence since 2007, according to national malaria surveillance data [1]. Although it is not clear what the main factors are driving these increasing malaria incidence trends, coverage of malaria control interventions is heterogeneous, and interventions are not delivered on a regular basis [2,5]. Indoor residual spraying (IRS) was first introduced in urban areas in 1994, whereby only ten districts in the country were selected. Insecticide-treated nets (ITN) deployment only started in 2004, primarily through donor funding. The south of Mozambique also took part in the Lubombo Spatial Development Initiative (LSDI), a tri-lateral initiative between the governments of Swaziland, South Africa and Mozambique that took place in the early 2000s, and officially finished its activities in 2011 [6,7].
Monitoring malaria trends in the different areas of the country and characterizing changes in clinical presentation and age distribution of malaria cases during shifts in transmission is crucial to inform policy. High-quality, health facility-based surveillance data, despite the potential selection bias due to differences in health-seeking behaviour and accessibility, allow finely characterizing of clinical malaria epidemiology of an area over time, offering accurate and detailed information that can complement information on trends observed through the national health management information system. Data from the Manhiça Health Research Centre (CISM) outpatient paediatric morbidity surveillance system, operating in the southernmost province of Mozambique, identified heterogeneous epidemiological patterns within its study area. This study aims to analyse the epidemiological trends of the outpatient clinical malaria cases from 2003 until 2013 in an area of high endemicity within the district if Manhiça, that experienced a drop in transmission after IRS, ITNs, and artemisinin combination therapy (ACT), followed by a later resurgence of clinical cases. Evidence from this analysis may offer the National Malaria Control Programme (NMCP) additional insight on the clinical malaria profile that could be observed in other areas of the country that have experienced similar changes in malaria transmission.

Study design
This is an analysis of surveillance data collected through the Ilha Josina Machel health post (IJMHP) outpatient paediatric morbidity surveillance system. All visits in children younger than 15 years old from 1 July 2003 to 30 June 2013 were included.

Study area and population
Ilha Josina Machel is a 188 km 2 river island in the confluences of the Incomati River, located about 50 km northeast of Manhiça village (Manhiça District, Maputo Province, Southern Mozambique) (Fig. 1). Approximately 15,000 people live currently in the area, 40 % of the population is under the age of 15 years, and 55 % are female. A third of the population is illiterate and the main occupation consists of subsistence agriculture. Houses are simple, with walls typically made of cane with thatched or corrugated roofs. Basic and preventive care is provided at the IJMHP, including antenatal and delivery services, but no admission facilities exist for sick patients, which are referred to Manhiça District Hospital in case of need. Additionally, as part of CISM's study area, Ilha Josina Machel is covered by a demographic surveillance system (DSS), explained in detailed elsewhere [8], that identifies all residents within the catchment area and provides accurate yearly-updated estimates of community denominators.
The region has two distinct seasons: a warm and rainy season from November to April, and a cooler and drier season the rest of the year. Rainfall trends varied throughout the study period, with wetter rainy seasons between 2002 and 2006, and slightly drier seasons between 2007 and 2012 (Fig. 2). Malaria transmission is perennial with seasonality, mostly due to Plasmodium falciparum, with a few sporadic cases of Plasmodium malariae, and with Anopheles funestus as the main vector. Malaria transmission was high in Ilha Josina Machel at the beginning of the study period. A cohort study to assess the P. falciparum re-infection rate was done in 2002, prior to the integration of Ilha Josina Machel in CISM's health demographic surveillance system (HDSS). The study included 288 children aged 6-18 months, who were followed up by active detection of P. falciparum infection during 6 weeks after curative therapy with sulfadoxine-pyrimethamine (SP). Eighty-seven percent of participants had been reinfected after 6 weeks, with an incidence rate of malaria infection of 4.35 episodes per child year at risk (95 % CI 3.34-5.69) (unpublished data).
IRS started in 2004 and was implemented once every year, except for 2011. IRS was conducted using carbamates (lambda-cyhalothrin) at first, and dichloro-diphe During this study period (July 2003-June 2013), firstline treatment for uncomplicated malaria was amodiaquine plus SP until September 2006, followed by ACT (initially artesunate plus SP and from 2009 onwards arthemether-lumefantrine, Coartem ® ).

Data collection
Since 2002, a passive case detection system was established at the IJMHP to cover all outpatient visits of children younger than 15 years, as part of the morbidity surveillance system run by CISM in all other health facilities in the DSS area and described in detail elsewhere [9]. Briefly, a standardized questionnaire gathering personal and demographic data (including the DSS unique permanent identification number, linking each patient to their geo-positioned household) as well as clinical information regarding the current illness is completed for each child attending the outpatient clinic. A qualified clinician records the physical signs and symptoms and their duration as referred by the child and/or accompanying guardian. A malaria blood smear (thick and thin films on the same slide) is performed from a finger prick for all children with measured fever (axillary temperature ≥37.5 °C), a history of fever in the preceding 24 h, or paleness, and blood is also collected into a microcapillary for packed cell volume (PCV) determinations. Final diagnosis/es using ICD-10 coding and treatment are recorded on the questionnaire by the clinician upon discharge or referral. Children fulfilling severity criteria are referred to the Manhiça District Hospital, where they are managed accordingly. Malaria cases are treated according to national guidelines.
Thin and thick blood films were read to quantify parasitaemia in the CISM laboratory according to standard procedures. Blood films were air-dried, Giemsa-stained, and examined to quantify parasitaemia using a light microscope fitted with a 100× oil immersion lens and a 10× eyepiece. Slides were declared negative only after counting 2000 leukocytes. Parasite numbers were converted to a count/µL by assuming a standard leukocyte count of 8000/µL for blood slides collected until the end of June 2009 and using a fixed volume method [10,11] from then onwards. Slides were first read at the health post using the 'plus system' [12] to guide case management. If the parasitaemia result was missing, results from the 'plus'-reading were used instead for the malaria case definition. PCV was measured using a micro-haematocrit centrifuge and a Hawksley reader (Hawksley & Sons Ltd, Lancing, UK) to assess levels of anaemia among cases.

Definitions
Clinical malaria was defined as the presence of fever (axillary temperature ≥37.5 °C) or history of fever in the preceding 24 h plus any asexual P. falciparum parasitaemia confirmed through blood smear microscopy. No age cut-off for the level of parasitaemia was used to increase the specificity of the definition, as this is the definition used at the outpatient clinic for all ages to diagnose and treat malaria.
The geometric mean parasite density (GMP) was calculated as the mean of the log-transformed parasite densities of clinical malaria cases per season. Mild anaemia was defined as a PCV between 25 and <33 %, moderate anaemia as a PCV of 15 to <25 % and severe anaemia as a PCV <15 %. Time was divided into 'seasonal' years (1 July-30 June of the following year) in order to capture the full transmission season and be able to account for any fluctuations of rain patterns or other environmental factors between years.

Data management and analysis
Outpatient questionnaires were double-entered into a database using a program written in Fox Pro (Microsoft Monthly rainfall data for every year in the study were accessed from the Climate Hazards Group InfraRed Precipitation with Station dataset [33]. Data were withdrawn for the coordinates of the district of Manhiça to obtain robust rainfall estimates. Mean rainfall estimates were calculated for every malaria season and compared through the study period Corp, Seattle, WA, USA) at CISM. Statistical analyses were performed using Stata 13.1 (Stata Corp, College Station, TX, USA).
The health post-based positive rate of clinical malaria (number of clinical malaria cases divided by total outpatient visits to IJMHP) and slide positivity rates (number of positive blood smears divided by total number of blood smears collected) were calculated for every seasonal year from 2003-2004 to 2012-2013.
Yearly mean age differences were calculated for clinical malaria cases and non-cases. A difference in difference (DID) analysis was conducted to compare the shifts in ages in time between malaria cases and non-cases. Mean age differences between years for malaria cases and noncases as well as the difference between those means were calculated using linear regressions with interaction terms between case category and season.
Yearly means of PCV and GMP by age group and their corresponding 95 % CI were also calculated to further explore the clinical and epidemiological trends of the disease through time. Analysis of variance (ANOVA) tests were conducted to compare means of normally distributed variables across seasonal years. PCV and GMPs were compared between age groups across the study period using multiple linear regressions and t tests. An increasing amount of outpatient visits with no PCV information was observed between 2007 and 2013 among clinical malaria cases, reaching 86 % of missing information for PCV among the malaria cases observed between 2012 and 2013. This was accounted for in the interpretation of results.
The DSS unique permanent identification number, which is recorded on the morbidity questionnaires when children attend the health post, was used to classify children as residents/non-residents in the study area. Minimum community-based incidence rates (MCBIRs) were calculated as the age-specific yearly number of malaria cases in children resident in Ilha Josina Machel who visited the outpatient clinic divided by the total child years at risk (CYAR) for that age group and year. CYARs were estimated from CISM's DSS databases and children did not contribute to the numerator or denominator for an arbitrary period of 28 days after each episode of malaria. They were also excluded when they were outside the study area or after death. Clinical malaria incidence rate ratios (IRR) were calculated through Poisson regression to further explore the relative shift in clinical malaria incidence over the years within each age group. IRRs were also calculated for every age group separately, which permitted better characterization of the age-specific trends of the disease through time.  (Fig. 3).

Results
The mean age of clinical malaria cases presenting to the IJMHP steadily increased for every seasonal year, shifting from 2.9 (SD 2.  (Fig. 4a). The DID analysis to compare age shifting patterns between malaria cases and non-cases showed a significantly larger increase in the mean age of patients with malaria than those with other pathologies. In general, the increase in malaria case ages was significantly higher for every year compared to the previous one except for 2007-2008 compared to 2006-2007 and in between 2009 and 2012, where the increase in age was similar in both groups (Fig. 4b).
Out of all clinical malaria cases attending the IJMHP between 2003 and 2013, 46.8 % were mildly anaemic, 13.4 % moderately anaemic and 0.7 % suffered from severe anaemia. A significant increase in the mean PCV among clinical malaria cases was observed, shifting from  Table 2 presents the clinical malaria MCBIRs and IRRs for every year of the study overall and by age group. Overall, children aged 1-4 years experienced the highest number of incident clinical cases per 1000 CYAR at any point in time except for the first year, and children aged 5-14 years were the age group with the lowest incidence during the majority of the study period.

Table 1 Summary characteristics of clinical malaria cases attending the IJMHP between 2003 and 2013
SPR slide positivity rate, SD standard deviation, PCV packed cell volume, GMP geometric mean parasite density a 86 % clinical malaria cases without PCV information   Presents the DID analysis that shows that the age difference from 1 year to the next between children with malaria (blue) was significantly larger than that for non-malaria cases (red). The green line represents the difference in mean age difference from 1 year to the next  [1]. Nevertheless, the national health management information system does not yet offer complete and accurate data that allows further characterization of the profile of malaria cases observed nationally, and data from Malaria Indicator Surveys are outdated (last performed in 2007).

Discussion
In Ilha Josina Machel, a sharper increase of clinical cases was observed particularly among children on their first to fifth birthday between 2007 and 2013, as well as those between the ages of 5 and 15 years. Overall, a shift in the age distribution of the clinical malaria cases towards older children was observed particularly among malaria cases, indicative of an alteration in the transmission intensity patterns in the area [13] and a potential consequent impact on the naturally acquired immunity (NAI) of the population.
The current analysis does not allow for any assessment of causality to evaluate which factors contributed to the decline and subsequent rise in transmission in the area. Rainfall patterns during the study period show a decrease in precipitation levels around the same years as malaria cases started to drop. This suggests that environmental factors may have had a considerable effect on the epidemiology of malaria in the past decades. Further spatiotemporal studies in the area shall be conducted to more accurately measure the impact of climatic factors on malaria through time. However, it is plausible that the introduction of new malaria control measures also played an important role in the initial decline in transmission. IRS was first deployed in Ilha Josina Machel in 2004-2005, according to district records, which correlates with the observed decrease in incidence between 2006 and 2009. However, an increase in the number of clinical malaria cases was later observed, which could be associated with a decrease in the IRS effectiveness due to low coverage (as reported between 2010 and 2013), or lack of continuity of IRS activities (such as in 2011). In fact, malaria resurgences have been observed in a number of different settings across the world since the beginning of the twentieth century, most of which are attributed to weakening of malaria control programmes due to resource constraints, and a subsequent inability to ensure continuity of IRS in places where it had already been delivered [14].
The the area at that time, considering that no changes were made in the reporting system and that malaria-testing rates remained constant throughout this period. The significant shift in the age distribution of clinical malaria cases towards older children observed in IJMHP supports this hypothesis, as in lower transmission settings young children will be exposed to less infectious bites and will take longer to build up a naturally acquired immunity, with older children having a relatively high incidence of clinical malaria [15,16].  [17][18][19][20][21]. The low-transmission period between 2007 and 2009 was followed by an immediate increase in the incidence of clinical malaria cases, particularly among 1-4 years olds and later also in the 5-14 years olds. During the lower malaria transmission period there was a decrease in the exposure to parasites in children (evident in the incidence rates among infants since 2007), delaying the acquisition of NAI in younger children and diminishing the premunition in older children that had already developed some NAI during the previous years [22]. Several studies, conducted to evaluate this question through actively interrupting exposure to the parasite using chemoprophylaxis at a very young age, have reported a similar rebound effect in the incidence of clinical malaria cases among children in the treatment arm, up to 1 years after treatment [15,[23][24][25]. The biological mechanisms driving the acquisition of NAI are still not well understood; however, evidence from mathematical models [26,27] and immunological studies [28,29] clearly identify age, exposure to the parasite and transmission intensity as important factors involved in this process.
The lower increase in the incidence rate of clinical malaria cases among infants may be partly explained by a drop in transmission; however other interventions such as the introduction of vector control activities targeting pregnant women and infants could have also contributed to an increased protection against parasite exposure among this age group.
The age-specific GMP trends observed in this study, with significant increases in parasitaemia in children older than 1 year after 2007 also show a slower acquisition of NAI with age and a lower 'anti-parasite' immunity [22]. In addition, the increase of mean PCV with age despite a  Fig. 6 Yearly incidence rates (in 1000 per CYAR) by age group, observed at the IJMHP. Yearly MCBIRs in 1000 per CYAR were calculated to assess clinical malaria incidence observed in infants (blue), 1-4 years olds (red), and 5-14 years olds (green), at the IJMHP between 2003 and 2013. CYARs were obtained from CISM's DSS, which is updated on a yearly basis raise in parasitaemia could possibly be a reflection of lower chronic infections or lower number of episodes per year, which cause less anaemia and does not allow the anti-parasite immunity to develop fast, resulting in the delay in the rate of NAI and the observed high parasitaemic infections.
Insecticide resistance may also be a contributing factor to the observed resurgence in the number of clinical malaria cases in the area. Mosquito resistance to pyrethroids was first reported in 2011 in a nearby sugar cane plantation of Maragra in the district of Manhiça [30]. Another study in a Northern province also reported wide resistance to pyrethroids and carbamates in 2007 [31]. Entomological data from CISM further supports these findings revealing less than 9.6 % sensitivity to pyrethroids among the A. funestus population of the district of Manhiça [32]. Insecticide resistance to pyrethroids has serious implications regarding the likely impact of LLINs towards reducing malaria transmission in the area, which since 2011 may have contributed to the increasing malaria incidence trends in the area.
The estimates obtained from this analysis are subject to a series of limitations. First, access to health care may have increased through the years. This could have caused an underestimation of the malaria burden in the area on the first years of the study, and a consequent overestimation of IRRs when comparing rates across years. Second, IJMHP data are also susceptible to the Hawthorne effect [16], as the presence of CISM, its impact on case management as well as research studies conducted, may affect the overall epidemiology of malaria in the study area of Manhiça. Particularly, Ilha Josina Machel was included as one of the sites where the RTS,S vaccine trials took place between 2002 and 2008. This analysis was also affected by missing data of PCV particularly in the last years of the study, which may have affected the estimates of severe and mild anaemia cases during these years.
Additionally, this study aims to provide a temporal description of the clinical malaria epidemiology in a health facility. Therefore, it did not conduct any sociodemographic or environmental risk factor analysis, which would enrich the interpretations withdrawn from the results obtained from this particular assessment. Separate studies should be conducted to respond to the hypothesis raised in this study, identify the main drivers of malaria transmission in the area and determine the social factors that may put a resident in the area at a higher risk of malaria than others.

Conclusion
This study shows a reduction in the clinical malaria incidence in a small area of initially high transmission in Southern Mozambique, and a subsequent resurgence of cases. Evidence from this study indicates that efforts to control malaria and reduce transmission need to be sustained over time in order to guarantee gains in the long term. Malaria control programmes with limited resources, such as the one in Mozambique, should consider malaria control interventions that ensure sustainability and target the most vulnerable populations, particularly infants and children under 5, at all times. Considering the costs and operational difficulty of this task, NMCPs should consider investing in innovative strategies to drastically interrupt transmission and eliminate the parasite reservoir from communities in a short period of time, and subsequently concentrate efforts to manage and eliminate remaining foci of transmission and imported infections.