Skip to main content

Nationwide evaluation of malaria infections, morbidity, mortality, and coverage of malaria control interventions in Madagascar

Abstract

Background

In the last decade, an important scale-up was observed in malaria control interventions. Madagascar entered the process for pre-elimination in 2007. Policy making needs operational indicators, but also indicators about effectiveness and impact of malaria control interventions (MCI). This study is aimed at providing data about malaria infection, morbidity, and mortality, and MCI in Madagascar.

Methods

Two nationwide surveys were simultaneously conducted in 2012–2013 in Madagascar: a study about non-complicated clinical malaria cases in 31 sentinel health facilities, and a cross-sectional survey (CSS) in 62 sites. The CSS encompassed interviews, collection of biological samples and verbal autopsies (VA). Data from CSS were weighted for age, sex, malaria transmission pattern, and population density. VA data were processed with InterVA-4 software.

Results

CSS included 15,746 individuals of all ages. Parasite rate (PR) as measured by rapid diagnostic tests was 3.1%, and was significantly higher in five to 19 year olds, in males, poorer socio-economic status (SES) quintiles and rural areas. Long-lasting insecticidal nets (LLIN) use was 41.7% and was significantly lower in five to 19 year olds, males and wealthier SES quintiles. Proportion of persons covered by indoor residual spraying (IRS) was 66.8% in targeted zones. Proportion of persons using other insecticides than IRS was 22.8%. Coverage of intermittent preventive treatment during pregnancy was 21.5%. Exposure to information, education and communication messages about malaria was significantly higher in wealthier SES for all media but information meetings. The proportion of fever case managements considered as appropriate with regard to malaria was 15.8%. Malaria was attributed as the cause of death in 14.0% of 86 VA, and 50% of these deaths involved persons above the age of five years. The clinical case study included 818 cases of which people above the age of five accounted for 79.7%. In targeted zones, coverage of LLIN and IRS were lower in clinical cases than in general population.

Conclusions

This study provides valuable data for the evaluation of effectiveness and factors affecting MCI. MCI and evaluation surveys should consider the whole population and not only focus on under-fives and pregnant women in pre-elimination or elimination strategies.

Background

Malaria control was marked in the last decade by a steep and almost constant increase in international funding. Domestic funding for malaria control increased as well in the same period, although less sharply. The successes of control programmes resulting from availability of resources have led some countries towards pre-elimination, including Madagascar in 2007 [1]. Unfortunately resources available for malaria control are reaching a plateau, especially in low-income countries that depend heavily on international funding [2]. Those countries also bear the highest burden of malaria [3]. In order for malaria control programmes to reduce the impact of this lethal disease, two major options exist: i) a new control tool that performs better at lower cost, or ii) ensuring that interventions deployed perform accurately. To evaluate such performance, individual effectiveness, coverage and community effectiveness need to be monitored [4]. In 2012–2013, a project aimed at measuring the coverage and effectiveness of malaria control interventions in Madagascar in terms of infection, morbidity and mortality, and the factors associated with their effectiveness as measured by both quantitative and qualitative approaches was set up. This project was named MEDALI (an acronym for Mission d’Etude des Déterminants de l’Accès aux Méthodes de Lutte antipaludique et de leur Impact), and took place in Madagascar. Herein are presented i) the design and methodology of the quantitative surveys of the MEDALI project; ii) a nationwide assessment of malaria infections, morbidity and mortality; and, iii) an evaluation of the coverage of malaria control interventions in all population groups.

Methods

A nationwide survey in 2012–2013 in 31 zones representing all transmission patterns of Madagascar was conducted. This survey included a cross-sectional study and a concurrent study about clinical cases. Clinical data, sociodemographic data and exposure to malaria control interventions have been collected.

Interventions studied

In Madagascar, malaria control relies on the following interventions: i) long-lasting insecticidal nets (LLIN) distribution; ii) indoor residual spraying (IRS) campaigns; iii) intermittent preventive treatment of pregnant women (IPTp); iv) testing of fever cases by rapid diagnostic tests (RDT) or microscopy; v) treatment of non-complicated malaria cases by artemisinin-based combination therapy (ACT); and, vi) information, education and communication (IEC) campaigns.

IPTp was introduced in 2004. In 2006, the national policy switched from chloroquine to ACT as treatment for uncomplicated malaria, together with the use of RDT. Artesunate/amodiaquine combination therapy is the first-line treatment, while artemether/lumefantrine is the second-line treatment with oral quinine as the alternative. In 2008, universal coverage with LLINs was scaled up at national level. While there was a half-century long history of IRS in Madagascar using dichlorodiphenyltrichloroethane (DDT), IRS with pyrethroids was deployed in 2008 in the Central Highlands and Fringe transmission patterns [5]. In 2009–2010, IRS was extended to some South and Western districts through the Global Fund to Fight AIDS, Tuberculosis and Malaria and the US President’s Malaria Initiative; while it was down-scaled to focalized IRS in the Central Highlands. Similarly to other interventions, IEC campaigns were reinforced in the 2000s [6].

Selection of study zones

Districts of Madagascar are divided in five main operational zones that correspond to the transmission patterns of Madagascar (Figure 1) [7]. The two coastal regions exhibit highest endemic patterns with a transmission lasting all year long in the East and more than six months per year in the West. In the central highlands, the transmission is unstable, and episodic or epidemic. In the Fringe areas, i.e., at intermediate altitudes, the transmission pattern is seasonal, lasting from November to May (rainy season). In the South, the period of transmission is short and episodic. Fringe, Central Highlands and the South are prone to outbreaks.

Figure 1
figure 1

Malaria transmission patterns in the districts of Madagascar, study sites and their population density.

The selection of study zones was based on a network of sentinel health centres (SHC) that were established in order to cover all the ecosystems of Madagascar for surveillance of fever-associated diseases [8]. Each location, where at least one SHC existed in 2012, was included in the study design defined 31 study zones. All malaria transmission patterns were represented: 13 study zones were located in the Western transmission pattern, seven in the East, five in the Fringe, four in the Central Highlands and two in the South (Figure 1). These patterns encompass, respectively, 21.0, 27.5, 13.7, 31.9, and 5.9% of Malagasy population.

Cross-sectional survey

The cross-sectional survey (CSS) aimed at studying the prevalence rates of malaria infection and exposure to factors that could be related to malaria, including malaria control interventions. The sample is made of clusters, stratified by transmission patterns.

Selection of study sites

For each study zone, two study sites, were defined to make a total of 62 study sites. To achieve this, two villages or fokontany, the smallest administrative subdivision of Madagascar, were randomly selected in a list of the fokontany surrounding the SHC: one close (i.e., centroid within a 2-km radius from the SHC) and the other further (i.e., centroid within a 10-15-km radius from the SHC), so as to include rural areas in the sample. In a given fokontany, investigators recorded the approximate number of households in each settlement and they included a number of households in each settlement proportionate to the total number of households in the whole fokontany. Households within a given settlement were selected according to a random pathway: a random point on a map of the settlement indicated the first household, and investigators then visited the second household on the right when stepping out from the previous household. Study teams were told to include at least 225 people of all age categories belonging to a minimum of 50 households in each of the 62 study sites, to make a total sample size of 13,950 individuals. If the study team could not include the expected number of individuals or households, or if they faced security and/or accessibility troubles, they investigated the closest fokontany in the list of randomized fokontany. Four teams of six surveyors and one supervisor investigated approximately eight study zones each, from September 2012 to January 2013. Coastal areas were investigated first, at the end of the dry season, in order to ensure accessibility to the study zones and to lessen the effect of seasonal transmission in the other transmission patterns. The latter were investigated at the beginning of the rainy season, i.e., just before or at the very beginning of the malaria transmission season.

Sample size calculation

The primary endpoint for sample size calculation was the proportion of the population infected by any Plasmodium species as diagnosed by RDT, and the proportion exposed to malaria control interventions. In a total sample size of 13,950, a two-sided 95% confidence interval for an expected proportion of population exposed to interventions of 50% and for an expected proportion of population infected by a Plasmodium of 5% will extend 0.8% and 0.4% from the observed proportion, respectively. In a sample size of at least 225 by site, a two-sided 95% confidence interval for an expected proportion of population exposed to interventions of 50% and for an expected proportion of population infected by a Plasmodium of 5% will extend 6.5% and 2.3% from the observed proportion in the sites, respectively.

Inclusion and data collection

Criteria for inclusion of individuals in the CSS were: i) member of or visitor in a household selected at random for inclusion in the CSS and whose head of household or his representative signed the household informed consent; ii) age ≥ six months; iii) signed individual informed consent including agreement for blood sampling; and, iv) ability to take a per os treatment in case of positive RDT. Geographic coordinates of each household were obtained using global positioning system (GPS).

The head of household or his representative answered a questionnaire about age and sex of all members of households and people having slept in the household the previous night (referred as ‘visitors’), housing features, living conditions, asset ownership, access to i) drugs retailer, ii) community health worker (CHW), iii) trained health staff, details of last IRS, bed net ownership, characteristics of the bed nets and use previous night, and verbal autopsy (VA) of household members deceased the previous year. Some characteristics of the bed nets (brand, positioning above sleeping places, and physical integrity) were checked visually by investigators. In one out of four households, the interviewee was asked to answer an additional questionnaire about social data: religion, ethnicity, social network, socio-economic self-perception, income, expenses, level of isolation, and migration. All included individuals were asked to answer a questionnaire about sociodemographic features, any fever episode during previous three months and case management in case of fever, bed net use over previous three months, anti-malarial drug intake over last year, mobility during previous three months, exposure to IEC messages about malaria (media, delay since last event), and, for females ≥15 years old, reproductive history, outcome (life/death) of all children, management of current or last pregnancy (antenatal care, IPTp).

In one out of four households, all individuals were asked about their professional activity and marital status, and an individual ≥15 years old was randomly selected to answer a knowledge, attitude and practice (KAP) questionnaire about health, pregnancy, fever, malaria, and malaria control interventions.

The VA questionnaire was designed from the World Health Organization’s international standard VA questionnaires [9], with added questions so as to collect all InterVA-3 indicators [10]. Criteria for answering the VA questionnaire were: specific informed consent, and delay since reported death between two and 12 months in order to respect a minimal mourning period while avoiding memory biases [11].

All answers were collected on notebooks with the help of a data entry form build on Microsoft Access® (Microsoft Office Professional Plus 2010, Redmond, WA, USA). Investigators were allowed to use paper questionnaires in case of software or hardware failure.

Collection of biological samples

All included individuals were punctured on finger or heel in order to perform a RDT (CareStart® Malaria, Access Bio Inc, Monmouth Junction, NJ, USA), to perform thick and thin blood smears, and to collect 100–500 μL capillary blood in an appropriate tube (MiniCollect® Z Serum Sep, Greiner Bio-One, Kremsmünster, Austria). Results of RDT were disclosed to individuals or their guardian after the interview, and people with a positive RDT result were proposed to take an anti-malarial treatment according to national guidelines. Thin smears were fixed with absolute methanol and serum tubes were centrifuged 10 min at 2,800 g on a daily basis. Serums tubes were kept at < 22°C until their transfer to the Institut Pasteur de Madagascar, where they were stored at -20°C.

Definition of variables

The individual coverage of LLINs was estimated by two means: i) whether the individual was listed among household members having slept under a LLIN the previous night according to the head of household or his representative, and, ii) according to the answer to the question “How many nights did you spend under a mosquito net during the last 3 months?” Answers to the latter question were categorized as following: “every night” versus any other answer (“1-6 time(s)/week”, “1-4 time(s)/month”, “less than once a month or never”). LLIN coverage was calculated for all individuals and not among individuals in household which owned a net only. In this manuscript the expression “coverage” refers to “use coverage” and not to “ownership coverage”.

A household was considered as covered by IRS if the head of household or his representative answered positively to the question “Did someone come in your household in the last 12 months to spray walls in order to fight against malaria?” and reported a delay since last IRS campaign ≤20 months so as to include only spraying during last campaign. For visitors, the IRS status at home was considered. Since the IRS campaign was late in the 2012–2013 rainy season, no household had been sprayed in the study sites at time of survey and no subset of households recently sprayed could be defined. In order to catch properly households covered by IRS, local designations (e.g., DDT) were used instead of, or in complement with, the term “spraying walls in order to fight against malaria”.

A woman was considered as covered by IPTp if she mentioned having taken at least two doses of sulphadoxine-pyrimethamine (SP) in the course of antenatal care during a ≥ four-month ongoing pregnancy or during her last pregnancy if she delivered within the previous 12 months. Males, females not of childbearing age (15–50 years old), nulliparous women, women pregnant < four months, and women who delivered more than 12 months previously were excluded from the calculation of IPTp coverage.

The following media were considered for the evaluation of exposure to IEC messages about malaria: radio (≤one month, >one to ≤ four months, or > four months since last message heard about malaria), poster (≤one month, >one to ≤ four months, or > four months since last poster seen), mobile video unit (MVU) (≤one year or > one year since last seen), television (≤four months or > four months since last message seen), leaflet or written press article (≤one year or > one year since last message read), and other media/exhibitions: whatever the delay spent since last event (e.g., hiragasy (traditional malagasy theatre), puppets, theatre, etc.).

A fever episode was defined as the self-report of any clinical sign of fever (e.g., hot body feeling, shivering) within the previous three months, whether or not accompanied by other clinical signs. A fever case management was defined as ‘appropriate’ as regards malaria if: i) either the patient underwent malaria diagnosis (RDT or microscopy, whether the result was negative or positive) by a trained health staff or a CHW (i.e., he/she was supposed to be treated by ACT in case of diagnosis of Plasmodium infection), or ii) he/she took ACT (whether prescribed or self-administered) as presumptive treatment when he/she could not have access to malaria diagnosis. A recent uptake of anti-malarial drugs was defined as the treatment of last fever episode with anti-malarial drugs or self-report of uptake of anti-malarial drugs during three months before the interview.

Household socio-economic status (SES) was calculated using the first principal component score based on variables related to house size, living conditions (toilet, waste, water, lighting, cooking fuel), asset ownership, and household utilities. The scores were converted to wealth quintiles for analysis (Additional file 1). Housing construction materials (walls, roof, floor) and presence of holes in housing structure (i.e., lack of mosquito-proof housing conditions) were kept apart from the principal component analysis. The head of household or his representative was also interviewed about domestic use of other insecticides than IRS inside the household, i.e. aerosol sprays, coils, or agricultural insecticides.

Data processing and analysis

GIS data processing

A raster image of the density of population in Madagascar was obtained from WorldPop/AfriPop database [12]. WGS 84 image was converted to UTM 38S so as to reconstruct a 100*100 m grid by resampling closest pixel from WGS 84. Then, for each household and the centroid of each study site of the CSS, the population value was obtained. Cut-offs to classify the population density of a given location as high, intermediate or low were determined by comparison of population densities and landscape characterized by surveyors when conducting CSS (‘urban’, ‘town or suburban’ or ‘rural’). All pixel values corresponding to the territory of Madagascar were counted in order to determine the proportion of the population in the country living in each of these three categories of population density.

Weights and sample design

In order to estimate standardized proportions representative of the situation in the whole country, counts and proportions and their confidence intervals were weighted to account for the complex sample design and for the adjustment on i) gender and age categories, according to data collected during the previous Demographic and Health Survey in 2008–2009 [13]; ii) malaria transmission pattern according to population data provided by the UN Office for the Coordination of Humanitarian Affairs [14]; and, iii) population density categories as described above, using the survey package in R version 3.0.2 [15]. Significance of the differences between proportions was calculated using χ2 test as it applies to both crude and weighted proportions, but using Rao and Scott corrections in the latter situation [15], and considering a 5% significance threshold.

VA data processing

Although the VA questionnaire was based on InterVA-3, VA data were processed with the InterVA-4 software, given the improvements made since the previous version [16]. This software addresses up to three most probable causes – or an indeterminate diagnostic – to each death in the database and the likelihood associated with each cause. Deaths having multiple causes were attributed as fractional deaths for particular causes, proportional to cause-associated likelihoods, as suggested [16].

Clinical malaria cases study

This clinical case study (CCS) aimed at evaluating the population affected by non-complicated malaria episodes. The design consisted of recruiting non-complicated malaria cases in health facilities belonging to the SHC network. As in the CSS, the sample of the CCS is made of clusters, stratified by transmission patterns.

Sample size calculation

The primary endpoint for sample size calculation was the exposure to malaria control interventions among patients presenting with clinical malaria due to any Plasmodium species. In a total sample size of 800, a two-sided 95% confidence interval for the proportion of population exposed to interventions will extend 3% from the observed proportion for an expected proportion of 50%.

Inclusion and data collection

Each of the 31 study zones was centred on one SHC. All these 31 SHCs were proposed to participate in the CCS. In the participating SHCs, all patients presenting with clinical malaria cases or their tutors submitted a short one-page questionnaire about socio-demographic data and exposure to malaria control interventions: LLIN, IRS and IPTp. Inclusion criteria were i) fever, i.e., axillary temperature ≥37.5°C [17] or self-reported symptoms of fever; ii) RDT or microscopy positive for any Plasmodium species; and, iii) informed consent of the patient or his/her tutor. Cases were retained in the database if they came from the same commune as individuals included in the CSS. Data were collected from September 2012 to August 2013.

Definition of variables

Exposure to malaria control interventions was defined the same way as for coverage in the CSS. The only exception was that a woman was exposed to IPTp if she mentioned having taken at least two doses of SP in the course of antenatal care during a ≥ four-month ongoing pregnancy or during last pregnancy if she delivered within the previous month (instead of “last 12 months”). For the comparison between cases and controls, the same definition was used in the CSS sample.

Ethical considerations

Both CSS and CCS followed ethical principles according to the Helsinki Declaration. Informed consent was obtained from the individuals or the parents/tutors of the children before inclusion. The CSS and CCS were approved by the National Ethic Committee of the Ministry of Public Health of Madagascar (approval #CNE 57/MSANP/CE of 24 July, 2012).

Results

Cross-sectional survey

Population sampled

Among the 62 fokontany selected at random, 49 could be directly investigated; among the 13 other fokontany, three were rejected for security reasons, five could not be accessed by survey teams, e.g., if a river was flooding, in two cases the authorities refused to cooperate, e.g., because of their beliefs towards blood sampling, in other two cases a misunderstanding caused the team to investigate the wrong fokontany, and the last one was actually located in another district than the one being investigated. All 13 fokontany were replaced by other randomly selected fokontany. In the 62 study sites, the population of three hamlets refused to be investigated: one because of an ongoing event (funeral), another because of local taboos (fady) and the last one refused collection of blood without further explanation. In these cases, only the rest of the fokontany was investigated.

Overall, 4,356 households were selected at random for interview. In 325 (7.5%), the household remained vacant during the survey. In 455 of the 4,031 remaining (11.3%), the head of household or his representative refused to participate (Figure 2). In 388 (85.3%) of the refusals, a reason was reported; most cited reasons for refusal were: not having time (35.3%), blood sampling refusal without further explanations (24.5%), blood sampling refusal for religious or taboos (fady) reasons (16.8%), and feeling healthy (i.e., they “don’t need medical investigation” 16.5%). All households that declined to participate were replaced with other households randomly selected from the same hamlet.

Figure 2
figure 2

Flow diagram of cross-sectional survey.

A total population of 19,921 individuals were members or visitors of the 3,572 households kept for analysis. Of the 16,234 eligible individuals, 561 (3.5%) either refused to be interviewed or sampled; in both cases, only sociodemographic data were collected. The final sample was 15,746 individuals.

Crude male/female ratio was 0.76; weighted male/female ratio was 0.98. Males were mainly under-represented in the crude sample in 20–39 year olds (Additional file 2). Crude proportion of visitors was 2.6% and weighted proportion was 2.5%. Regarding transmission patterns, the west was the most oversampled since 43.0% of the sample was located in this pattern while weighted proportion was 21.0%; Central Highlands and the East were the most subsampled since crude proportion were 12.8% and 22.3%, respectively, while weighted proportions were 32.4 and 28.1%. In the Fringe and the Southern patterns, crude proportions were 15.7 and 6.2% respectively, and weighted proportions were 13.1 and 5.4%.

Cut-offs of population density were set at two and 140 habitants/10,000 sq m. Considering this, 62.9% of the sample was estimated to live in rural areas (≤two habitants/10,000 sq m), 29.6% in small towns and suburban areas (>two and ≤140 habitants/10,000 sq m), and 7.6% in urban areas (>140 habitants/10,000 sq m). Weighted figures for these proportions were 59.1, 26.6 and 14.3%, respectively.

Low SES households were over-represented in the sample since the weighted proportions are inferior to the crude proportions in the two first quintiles, and the opposite is observed in the two last quintiles (Additional file 3). In the Central Highlands, high SES quintiles are prominent while the opposite is observed in the West, where more than half of the population belongs to the first (lowest) SES quintile. Wealthiest SES quintile is more frequently observed in settings where population density is intermediate or high (weighted proportions: 32.6 and 61.2%, respectively) than in low population density areas (weighted proportion: 9.6%,Additional file 4). Conversely the poorest quintile is almost absent from high-density areas (weighted proportion: 0.7%) while it represents 6.5 and 25.0% in intermediate-density and low-density areas, respectively (Additional file 4).

The largest group in terms of level of education among people ≥15 years old had accomplished primary school only (weighted proportion: 40.3%) and this situation was found in all transmission patterns (range of weighted proportions: 37.8-49.3%) except in the West (weighted proportion: 23.1%, Additional file 5). In the Western transmission pattern, the proportion of individuals ≥15 years old without any formal education was prominent (weighted proportion: 58.1%) while in other patterns this category was minor (range of weighted proportions: 4.9-18.1%, Additional file 5). The proportion of people ≥15 years old that had either no education or primary grade only was more important in areas of low population density (weighted proportion: 63.3%) than in other areas (weighted proportions: 41.3 and 32.9% in intermediate and high population density respectively, Additional file 4).

Among the 15,476 persons tested, 552 had a positive RDT result. Among these 552 samples, six slides were damaged and unreadable. Among the 546 remaining, all were examined by microscopy: 228 were found positive for Plasmodium falciparum, one positive for Plasmodium vivax and 317 negative. A sample of 459 slides from people with a negative RDT result was examined by microscopy: two were positive for P. falciparum. Overall, crude proportion of positive RDTs – or parasite rate (PR) – was 3.6% (95% CI [3.3-3.9]) and weighted PR was 3.1% (95% CI [2.1-4.5]). Weighted PR2–10 or PR between two and ten years old, an indicator for the intensity of transmission [18], was 4.2% (95% CI [2.7-6.6]). Another key indicator [19] is the proportion of children aged six to 59 months with malaria infection: the weighted value of this proportion was 2.8% (95% CI [1.8-4.4]) in the present study. The PR increased with age up to the age of nine then decreased progressively (p < 0.0001, Table 1). Infections in age groups between five and 19 years accounted for 57.8% of all Plasmodium infections. The weighted PR was higher in males (3.6%, 95% CI [2.3-4.9]) than in females (2.6%, 95% CI [1.5-3.7], p = 0.0004, Table 1). This difference in gender was observed in all age categories, except in under two year olds (Additional file 6).

Table 1 Weighted proportions of individuals with a positive rapid diagnostic test by sociodemographic characteristics

Weighted PR was highest in the Eastern (4.9%) and Western (4.6%) transmission patterns. In the Southern transmission pattern, PR was 3.1%; in the Central Highlands and Fringe transmission patterns, it was 1.5 and 0.8%, respectively. These differences in PR between transmission patterns were significant (p = 0.0119). PR were higher in poorer SES quintiles (p = 0.0001) and in lower education levels although not significantly in the latter case (p > 0.1, Table 1). In rural areas (i.e., low population density), PR were higher than in the other two categories but this difference was not significant (p = 0.10, Table 1).In the Western and Eastern transmission patterns, the PR was the highest among five to nine year olds (8.0 and 8.4%, respectively, Figure 3). In the South, PR was the highest among ten to 14 year olds (7.5%). In the Central Highlands and Fringe, PR was high among under one year olds (2.4 and 2.6%, respectively) but also higher among 15–19 year olds (2.3 and 0.9%, respectively) than among other age groups (Figure 3). There was a shift in higher PR toward the oldest in the lowest malaria transmission conditions.

Figure 3
figure 3

Weighted parasite rates according to age and transmission patterns. Weighted proportions of individuals with a positive RDT result by age categories centered on median age for the category, and by transmission pattern.

Coverage of malaria control interventions

Bed nets

The vast majority of bed nets were LLINs (5,461/6,267, 87.1%). Only nine bed nets were described as having been treated with insecticide at home. Most bed nets were used by one to three persons the previous night (65.8%) and a quarter (25.5%) remained unused. An important proportion of bed nets had been received from mass distribution (65.2% of all bed nets, 74.7% of LLIN). Figures of coverage according to transmission pattern or population density were similar whether measured as having slept under LLIN or non-impregnated bed net (NIBN) the previous night, or every night during the previous three months, although the proportion of every-night LLIN users was consistently lower than previous-night LLIN users (Additional file 7). In the Central Highlands, where no mass distribution of LLIN occurs, LLIN use was lower (9.4%, 95% CI [3.4-15.3]) than in other transmission patterns (57.2%, 95% CI [52.0-62.4], p < 0.0001, Table 2 and Figure 4).

Table 2 Weighted proportions of long-lasting insecticidal net users on the previous night by sociodemographic characteristics
Figure 4
figure 4

Proportions of long-lasting insecticide net users, non-impregnated bed net users and no bed net users on the previous night. Weighted proportions of individuals having used a LLIN, a NIBN or no bed net on the previous night, by transmission pattern.

Crude proportion of LLIN users the previous night was 51.3% (95% CI [50.5-52.1]) and weighted proportion was 41.7% (95% CI [33.2-50.7]). LLIN use was highest among under-five age groups, and among adults over 19 years old (p < 0.0001, Table 2). Males used LLIN less frequently (40.4%) than females (42.9%, p = 0.007, Table 2). Proportions of LLIN users were higher in low SES categories and in low education levels, but these differences were significant only in the latter situation (p > 0.1 and p = 0.04, respectively, Table 2). Differences of coverage between categories of population density were not significant (p > 0.1, Table 2).

Key indicators of LLIN coverage [19] included: i) LLIN ownership, i.e. the proportion of households with at least one LLIN: 64% (weighted, 95% CI [52.9-73.8]); ii) the proportion of households with at least one LLIN for every two people: 22.7% (weighted, 95% CI [17.7-28.8]); iii) the proportion of population with access to a LLIN within their household, i.e., the proportion of individuals who could sleep under an LLIN if each LLIN in the household is used by two people [19]: 42.7% (weighted); iv) the proportion of children under five years old who slept under a LLIN the previous night: 49.1% (weighted, 95% CI [39.9-58.4 ]); v) the proportion of pregnant women who slept under an insecticide-treated net the previous night: 38.4% (weighted, 95% CI [29.0-48.8]); and, vi) the proportion of existing LLIN used the previous night: 64.9% (crude, 95% CI [63.7-66.1]).

The proportion of NIBN users was largely inferior to the proportion of LLIN users: crude proportion was 6.9% (95% CI [6.6-7.4]) and weighted proportion was 6.7% (95% CI [4.9-9.0]). The differences in age groups and sex were the same as with LLIN users, but proportions of NIBN users increase with SES quintiles, level of education and population density, contrary to what was observed with LLIN users (Additional file 8). Proportions of NIBN users differed also between transmission patterns, but the figure is different from proportions of LLIN users (Additional file 8 and Figure 4).

IRS

The crude proportion of persons covered by IRS the previous 12 months was 27.6% (95% CI [26.9-28.3]) and weighted proportion was 30.2% (95% CI [21.4-40.8]). When considering only the zones targeted by IRS campaigns, i.e., Fringe, Southern and parts of Highlands and Western transmission patterns, the weighted proportion is 66.8% (95% CI [61.9-71.7]). Since all people living in a household having received IRS were considered as covered, there was no gradient between age groups. Differences between infants and older children or adults were significant (p = 0.01). There was no significant difference between sex (p > 0.1, Table 3). Urban areas, wealthiest SES quintiles and high levels of education had lower coverage proportions (p = 0.04, p = 0.003, and p = 0.02, respectively, Table 3).

Table 3 Weighted proportions of sample covered by indoor residual spraying in the previous 12 months by sociodemographic characteristics

Key indicators of vector control coverage [19] also included: i) the proportion of households with at least one LLIN and/or sprayed by IRS in the previous 12 months: 78% (weighted, 95% CI [69.7-84.5]); and, ii) the proportion of households with at least one LLIN for every two people and/or sprayed by IRS within the previous 12 months: 47.5% (weighted, 95% CI [40.0-55.1]).

Domestic use of insecticides

The crude proportion of persons living in households where they used other insecticides than IRS was 17.7% (95% CI [17.1-18.3]) and weighted proportion was 22.8% (95% CI [17.9-28.5]). Upper SES quintiles, high levels of education and central transmission patterns (Highlands, Fringe) displayed much higher proportions of domestic use of insecticides (p = 0.0003, p < 0.0001 and p < 0.0001, respectively, Table 4).

Table 4 Weighted proportions of sample in households reporting domestic use of insecticides by sociodemographic characteristics

IPTp

The crude proportion of women having received ≥ two doses of IPTp among pregnant women over four months of gestational age or women of childbearing age (i.e., 15–49 years old) having delivered in previous 12 months was 23.0% (95% CI [20.3-26.0]), and weighted proportion was 21.5% (95% CI [17.2-26.5]). When considering only the zones targeted, i.e., all transmission patterns but Central Highlands, the weighted proportion was 25.3% (95% CI [20.0-30.5]); in the South, coverage was also lower than in other transmission patterns (p = 0.0005, Table 5). Age, SES, education, and population density were not significantly associated with IPTp coverage (p > 0.1, p > 0.1, p > 0.1, and p = 0.07, respectively, Table 5). The proportion of women having received ≥ three doses of IPTp among women of childbearing age having delivered in previous 24 months was 8.4% (weighted, 95% CI [6.7-10.6]).

Table 5 Weighted proportions of women receiving two doses of intermittent preventative treatment in pregnancy by sociodemographic characteristics

IEC campaigns

Most cited medium for IEC messages about malaria was the radio (weighted proportion, 37.1%, 95% CI [32.5-41.8]), and second most cited was posters (19.6%, 95% CI [14.8-25.5]). Other media include television (13.4%, 95% CI [10.0-17.9]), information meetings (6.1%, 95% CI [4.8-7.6]), reading documents and/or articles (3.5%, 95% CI [2.6-4.7]), MVU (2.4%, 95% CI [1.6-3.5]), and other exhibitions such as hiragasy, i.e., traditional theatre (1.4%, 95% CI [0.8-2.2]). Exposure to radio IEC messages about malaria is significantly lower in younger ages, lower education levels and urban areas (p < 0.0001, p < 0.0001 and p = 0.004, respectively, Table 6); differences between sex and transmission patterns are not significant (p = 0.06 and p > 0.1, respectively, Table 6). Exposure to IEC messages was significantly higher in upper SES for most media (p = 0.01, p < 0.0001, p = 0.001, p = 0.007 for radio, TV, documents, and MVU respectively), non-significantly higher in upper SES for posters and ‘others’ (p = 0.1 and p = 0.05, respectively), and non-significantly higher in lower SES for information meetings (p = 0.07, Table 6). Exposure to information meetings was significantly higher in areas of low population density (p = 0.04, Table 6).

Table 6 Weighted proportions of individuals having received information, education and communication messages about malaria by sociodemographic characteristics

Fever case management

The weighted proportion of people who declared having presented a fever in previous three months was 12.1% (95% CI [10.5-14.0]). Data about 1,380 fever case managements was collected. Among these, 37.6% sought advice or treatment from a professional health worker or a CHW, 13.8% had a finger or heel stick for malaria diagnosis, and 17.6% received a treatment against malaria. Among people having consulted at a health centre or hospital, 71.5% resorted to the public sector. Among persons having received an anti-malarial, 31.4% received ACT. Overall, whether febrile or not, 2.7% of the individuals interviewed declared intake of a treatment against malaria in the previous three months. The weighted proportion of fever case managements considered as appropriate as regards malaria was 15.8% (95% CI [12.7-19.5]). The crude proportion was 17.5% (95% CI [15.6-19.7]). Proportions of appropriate case managements (ACM) were significantly higher in younger ages and in transmission patterns with higher PRs (p = 0.008 and p = 0.046, respectively, Table 7), but differences between sex, SES quintiles, level of education, and population density were not significant (all p > 0.1, Table 7).

Table 7 Weighted proportions of appropriate case managements by sociodemographic characteristics

Data about case management in 258 children under five years old who had fever in the previous two weeks were collected. Key indicators of case management [19] also included: i) the proportion of children under five years old who had fever in the previous two weeks and who had a finger or heel stick: 14.2% (weighted, 95% CI [8.3-23.1]); ii) the proportion of children under five years old who had fever in the previous two weeks for whom advice or treatment was sought from a professional health worker or a CHW: 38.8% (weighted, 95% CI [32.1-45.9]); and, iii) the proportion receiving ACT, among children under five years old who had fever in the previous two weeks who received any anti-malarial drugs: 37.2% (weighted, 95% CI [16.3-64.4]).

Mortality attributable to malaria

Among the 3,516 households interviewed for the occurrence of death in the household during the previous 12 months, 114 deaths were noted. Since 18,597 individuals were included in these 3,516 household, crude mortality rate is evaluated to be 6.1% (95% CI [5.1-7.3]). No explicit refusal to conduct VA was noted, but in eight cases, data were missing. In 20 cases, the death occurred outside the two to 12 months’ delay before the date of interview and no VA was conducted. The VA data from 86 deaths was collected, most of them occurring in elders or in under-fives (Figure 5).

Figure 5
figure 5

Age and sex profile of 86 deaths investigated by verbal autopsy.

After analysis by InterVA-4 software, the cause-specific mortality fraction (CSMF) of malaria proportional to cause-associated likelihood was 14.0% (95% CI [7.7-23.5]). Malaria ranked second CSMF after cardiovascular diseases and was attributed as most probable cause of death in 12 cases and as second most probable cause of death in two cases (Table 8). Apart from infants and neonates, there were deaths attributable to malaria in all age groups, but CSMF were higher in the one to four and five to 14 years age groups (p < 0.0001, Table 8). CSMF was 43.6% (95% CI [22.3-67.4]) before the age of five, and 8.7% (95% CI [3.5-18.8]) after that age. Only half of the malaria-attributed deaths occurred before the age of five (7/14).

Table 8 Age distribution of deaths and deaths attributable to malaria

Clinical case study

From September 2012 to August 2013, 6,771 malaria cases were observed in all 34 SHCs (Figure 6). Among those cases, 4,223 individual sheets (62.4%) were filled. Twenty-eight SHCs participated in the CCS and sent 1,585 filled questionnaires, among which 923 (58.2%) were individuals coming from same communes as individuals included in the CSS. Among those, 818 questionnaires (88.6%) were correctly filled, from 25 SHCs.

Figure 6
figure 6

Monthly number of malaria cases detected in the fever sentinel network. Monthly number of malaria cases detected in SHCs during the period of data collection, per transmission pattern.

In transmission patterns with high incidence, the age of clinical cases was younger than in transmission patterns with lower incidence (p < 0.0001, Table 9). People above the age of five accounted for 79.7% of all clinical cases. Males were more affected than females (p = 0.01), and there were no significant differences between transmission patterns (p > 0.1, Table 9). The proportion of LLIN users on the previous night was 52.1% (95% CI [48.6-55.5]), and the difference between transmission patterns was close to significance threshold (p = 0.0503, Table 9). Excluding Central Highlands, where no LLIN distributions occurred, proportion of LLIN users among clinical cases was 52.4% (95% CI [48.9-55.9]). The overall proportion of clinical cases living in households having received IRS the previous year was 11.6% (95% CI [9.4-13.9]), and this proportion in IRS-targeted zones was 48.4% (95% CI [41.0-55.8]). There were 11 women pregnant (second or third trimester) or having delivered the previous month; three of whom had received IPTp (Table 9).

Table 9 Distribution of sentinel health centres and clinical malaria cases per transmission patterns

Discussion

These results mainly confirm figures found in Malaria Indicator Survey (MIS) in Madagascar or other surveys, e.g., PRs are higher in Eastern and Western transmission patterns, PR increases with age up to the age of five, and coverage of malaria control interventions remain suboptimal, especially IPTp [7, 20]. These results also confirm that SES, education and environmental factors, such as population density, are all related to malaria infection and are closely intertwined [21]. Despite better LLIN coverage of individuals belonging to lower SES or educational level, social inequalities regarding malaria infection are striking. These inequalities can be explained by better access to prophylaxis (e.g., purchase and use of NIBN, as results show) and treatment, better housing quality and better nutritional status [21], and deserve specific investigation. It cannot be excluded that the slightly higher rates of appropriate fever case management in poorer SES quintiles are due to indication bias since higher ACM rates are observed in high-transmission areas, where higher proportions of low SES quintiles were observed. The lower IRS coverage in high-density areas is related to the absence of IRS campaigns in the centres of main cities in the Central Highlands.

PR were inferior to figures observed during MIS in 2011 and 2013 (Table 10). Two main factors could explain this difference: i) in the present study, the sampling took place between September and January, i.e., at the end of the dry season and beginning of rainy season, while MIS took place in November-December [7] or February-March [20], i.e., during the rainy season when incidence of malaria is higher; ii) the present study covers the whole country, including high altitude areas and big city centres that are supposed to be malaria free and not investigated in MIS [7, 20]. Seasons and inclusion of all areas also explains the lower LLIN coverage that was observed (Table 10), since LLIN use may depend on density of mosquitoes according to qualitative results [Mattern C, Raboanary E, Kesteman T, Rogier C, Piola P, Pourette D: Malaria is not a Problem. Consequences on Fever Case Management and Use of Bed Nets in Madagascar. Manuscript submitted] and drops in the Central Highlands where there is no distribution of LLIN. This result also mitigates the optimism arising from MIS results as regards LLIN coverage since a seasonal decrease in LLIN use can contribute to Plasmodium transmission during the dry season in areas where Anopheles breeding sites could persist, e.g., in rice fields and permanent water collections. These results regarding IRS, IPTp and fever case management stand between or are close to MIS results (Table 10).

Table 10 Main indicators found from the Malaria Indicator Survey and the present study

At first glance, this study may look like a MIS. Nevertheless, this study presents noteworthy differences. First, the population as a whole was investigated, not excluding geographical regions as it is done in MIS in Madagascar [7, 20], and above all not focusing on under-fives and women. Therefore it appeared that individuals between five and 19 years old presented higher PR, and represent probably a major reservoir and should be the main target of malaria control interventions in the context of pre-elimination in Madagascar. Concurrently, this age group presented lower coverage of LLIN, indicating a way to reduce this reservoir if coverage with LLIN was uniform between age groups. This result confirms that control programmes and evaluation surveys such as MIS should drop target groups (e.g. under-fives and pregnant women) and encompass the whole population when it comes to reducing transmission. While this policy is stated in international guidelines [22], national or local policies may stick to target groups [1]. In areas other than Madagascar, characterized by higher transmission, the reservoir may involve younger age groups [23]. It was evidenced that males of nearly all age groups were more frequently infected by Plasmodium as measured by RDT. This phenomenon could be explained by a biological susceptibility to infection, e.g., since lower serum iron concentration in females might protect them from Plasmodium infection [24], but also by differences in behaviour between gender that could deserve special attention as a target for IEC messages. Big city centres in Central Highlands and areas located at high altitude were included for two main reasons: i) SHCs in these areas detect few, but some, malaria cases that are not imported [personal communication, Randrianasolo L, Piola P, and Boyer S, Institut Pasteur de Madagascar]; and, ii) severe outbreaks occurred in these areas in the late 1980s and may occur again [25]. Recent detection of locally acquired clinical cases of malaria in the surroundings of the capital Antananarivo during the dry season of 2014 [personal communication, Randrianasolo L, Piola P, and Boyer S, Institut Pasteur de Madagascar] confirms that a priori any region of the country is permanently malaria-free.

Another major difference from MIS is the collection of data about malaria morbidity and mortality. As far as we know, this study is the first to encompass investigation of Plasmodium infection, control interventions, clinical cases, and malaria-attributable mortality. Results show that, similar to age profiles of PR, clinical cases occur in under-fives in 20.3% of cases, and that half the casualties could involve people over the age of five. Data from the CCS showed that in target zones coverage of LLIN and IRS were lower in clinical cases (52.4 and 48.4%, respectively) than in the overall population (57.2 and 66.8%, respectively). Adjusted ORs and protective effectiveness of malaria control interventions will be provided in a separate publication.

Results from a concurrent qualitative survey in four of the 31 settings investigated in the present survey suggest that the effect of IEC campaigns might be mitigated by predominance of informal channels in what people retain about malaria [Mattern C, Raboanary E, Kesteman T, Rogier C, Piola P, Pourette D: Malaria is not a Problem. Consequences on Fever Case Management and Use of Bed Nets in Madagascar. Manuscript submitted], which is consistent with the apparent low coverage of conventional media that was observed in quantitative results, and with findings in other countries [26]. The same qualitative survey also suggest that mosquitoes density affect LLIN use [Mattern C, Raboanary E, Kesteman T, Rogier C, Piola P, Pourette D: Malaria is not a Problem. Consequences on Fever Case Management and Use of Bed Nets in Madagascar. Manuscript submitted], which can explain that LLIN use in MIS is higher than in the present study (Table 10), since the presence of mosquitoes varies with seasons.

In Madagascar, hospital mortality attributable to malaria ranked fourth in 2011 [27]; in the present study it ranked second in overall mortality. In 2011, hospital mortality attributable to malaria was 8% among under-fives and 2% above this age; in the present study, overall mortality attributable to malaria in these age groups were 43.6 and 8.7%, respectively. These differences are, at least partially, explained by the fact that most malaria-attributable deaths occur at home and not in hospital [28], and because causes of death from VA are obtained differently from inpatient death certificates [11]. Nevertheless, the sample for measuring malaria-attributable mortality is limited since malaria was attributed as most probable cause by InterVA-4 software in only 12 of the 86 deaths analyzed, and as second most probable cause in two additional deaths. The performances of the previous version of this software has been suboptimal, especially regarding malaria [2931]. The current version of InterVA seems to perform better, similar to the assignment of causes of death by physicians [16, 3234], but all tools for interpretation of VA results are poorly accurate in the identification of death attributable to malaria [34]. In the absence of training dataset for the assignment of causes of death, InterVA is the only automated option or analysing VA data [16]. Because of the small size of the sample, and because vital registration system of Madagascar is not comprehensive, results were not weighted. As far as the authors know, no study about causes of death has ever been conducted across Madagascar. This study, with the limitations cited above, thus provides baseline results for CSMF in the general population of Madagascar.

Finally, this study provides information about factors affecting Plasmodium infection that were not collected in MIS. For example, domestic use of insecticides is widespread in Madagascar since 22.8% of the population live in households where insecticides other than IRS are used. This may be one of the factors explaining socio-economic differences in terms of Plasmodium infection, together with already documented factors such as access to prophylaxis and treatment, housing quality and nutritional status [21].

Numerous KAP surveys regarding malaria have been conducted, including MIS, but the exposure to IEC messages and analysis of media through which those messages are disseminated have less been investigated [27, 35, 36]. Results show that people belonging to upper SES and the more educated people are more exposed to, or better retain IEC messages about malaria from all media, except information meetings. Remarkably, most media broadcast better in areas of low population density, i.e., rural areas where malaria transmission is the highest, than in high-density areas, i.e., urban areas. As expected, television was an exception but also, more surprisingly, MVU that is designed to reach remote settings [37]. This result suggests that the design of IEC messages should pay more attention to reduce socio-economic inequities. The association between exposure to IEC messages and better coverage of malaria control interventions remain to be investigated, and will be presented later in a separate publication.

This study also presents various shortcomings. First, only RDT results and no thick smear results are available. Since 58.1% (317/546) of slides from individuals with a positive RDT results were found to be negative by microscopy, RDT result was selected for the assessment of Plasmodium infections. This high proportion of negative slides can be explained either by a lack of sensitivity of microscopy in the context of low Plasmodium transmission and low parasite density in this season, or by a problem in the preparation of slides. This second explanation is more plausible because the agreement between the RDTs used in the present study and microscopy has been demonstrated earlier in Madagascar [38], and because the preparation of the slides were performed in the households of participants, thus often on uneven surfaces. Separate results for non- P. falciparum infection were not presented because 99.6% of slides positive by microscopy exhibited P. falciparum, and only one slide was found positive for P. vivax. It was thus estimated that a positive RDT result was a reasonable proxy for P. falciparum infection.

Since clusters were chosen in the vicinity of health facilities belonging to fever sentinel network, the sample cannot be considered as strictly representative of the country. Additionally, numerous absences, refusals and/or impossibilities to investigate places and households (21.0% of fokontany and 17.9% of households) were noted, which could introduce a bias in the sampling. These pitfalls were overcome by applying weighting by age, sex, transmission pattern, and population density to the data, as described in Methods. As a consequence, the weighted figures observed were in rather good accordance with previous and subsequent studies (IRS, IPTp, and case management, Table 10). This establishes the proof-of-concept of the methodology used in the present study for the evaluation of nationwide indicators.

Another shortcoming is the collection of data about clinical cases in SHCs that belong to the public sector. The population of clinical cases is thus blinded for consultants in private sector, which may introduce a bias although this proportion was less than 30% in the CSS sample. Additionally, no excessive proportions of under-fives or pregnant women was observed in clinical cases, although they benefit from free care in the public sector. It was thus assumed that, if there is any bias, its importance may be limited.

Conclusions

The original design of this study provided concurrent results about Plasmodium infection, malaria morbidity and mortality attributable to malaria at a nationwide level. Results largely confirm facts described elsewhere, but some deserve to be underlined. First, the variation in LLIN coverage may be largely overestimated if only one season is considered while collecting indicators. Second, when malaria control interventions not only try to reduce burden but also aim to reduce Plasmodium transmission, other subgroups than under-fives and pregnant women, or even the entire population, should be considered as targets for control interventions and evaluation surveys. In addition, there was concern that people above the age of five in Madagascar may bear an important burden of malaria in terms of infection, morbidity and mortality. Third, important socio-economic inequalities have been found, in the prevalence of infection as described previously, but also in exposure to IEC messages and in domestic use of insecticides. The authors plead for more attention to these issues while public health policy makers design and implement malaria control interventions or conduct evaluation surveys.

Abbreviations

ACM:

Appropriate case management

ACT:

Artemisinin-based combination therapy

CCS:

Clinical case study

CHW:

Community Health Worker

CSMF:

Cause-specific mortality fraction

CSS:

Cross-sectional survey

DDT:

Dichlorodiphenyltrichloroethane

GPS:

Global positioning system

IEC:

Information, Education and Communication

IPTp:

Intermittent preventive treatment of pregnant women

IRS:

Indoor residual spraying

KAP:

Knowledge, attitude and practice

LLIN:

Long-lasting insecticidal nets

MEDALI:

Mission d’Etude des Déterminants de l’Accès aux Méthodes de Lutte antipaludique et de leur Impact

MIS:

Malaria Indicator Survey

MVU:

Mobile video units

NIBN:

Non-impregnated bed net

OR:

Odds ratio

PR:

Parasite rate

RDT:

Rapid diagnostic tests

SES:

Socio-economic status

SHC:

Sentinel health centres

SP:

Sulphadoxine-pyrimethamine

VA:

Verbal autopsy.

References

  1. Ministère de la Santé Publique: Plan Stratégique de Lutte Contre Le Paludisme 2007–2012 - Du Contrôle Vers L’élimination Du Paludisme À Madagascar. 2007, Antananarivo, Madagascar: Ministère de la Santé Publique, 1-54.

    Google Scholar 

  2. Global Malaria Programme WHO: World Malaria Report 2013. 2013, Geneva: WHO

    Google Scholar 

  3. Murray CJL, Rosenfeld LC, Lim SS, Andrews KG, Foreman KJ, Haring D, Fullman N, Naghavi M, Lozano R, Lopez AD: Global malaria mortality between 1980 and 2010: a systematic analysis. Lancet. 2010, 379: 413-431.

    Article  Google Scholar 

  4. Lengeler C, Snow R: From efficacy to effectiveness: insecticide-treated bednets in Africa. Bull World Health Organ. 1996, 74: 325-332.

    PubMed Central  CAS  PubMed  Google Scholar 

  5. Ratovonjato J, Randrianarivelojosia M, Rakotondrainibe ME, Raharimanga V, Andrianaivolambo L, Le Goff G, Rogier C, Ariey F, Boyer S, Robert V: Entomological and parasitological impacts of indoor residual spraying with DDT, alphacypermethrin and deltamethrin in the western foothill area of Madagascar. Malar J. 2014, 13: 21-10.1186/1475-2875-13-21.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Ministère de la Santé Publique, Programme National de Lutte contre le Paludisme: Revue Du Programme Paludisme À Madagascar. 2011, Antananarivo, Madagascar: Ministère de la Santé Publique

    Google Scholar 

  7. Institut National de la Statistique, Programme National de Lutte contre le Paludisme, ICF International: Enquête Sur Les Indicateurs Du Paludisme À Madagascar 2011. 2012, Antananarivo, Madagascar: Institut National de la Statistique

    Google Scholar 

  8. Randrianasolo L, Raoelina Y, Ratsitorahina M, Ravolomanana L, Andriamandimby S, Heraud J-M, Rakotomanana F, Ramanjato R, Randrianarivo-Solofoniaina AE, Richard V: Sentinel surveillance system for early outbreak detection in Madagascar. BMC Public Health. 2010, 10: 31-10.1186/1471-2458-10-31.

    Article  PubMed Central  PubMed  Google Scholar 

  9. World Health Organization: Verbal Autopsy Standards - Ascertaining and Attributing Cause of Death. 2007, Geneva: WHO Press

    Google Scholar 

  10. Fottrell E, Kahn K, Ng N, Sartorius B, Huong DL, Van Minh H, Fantahun M, Byass P: Mortality measurement in transition: proof of principle for standardised multi-country comparisons. Trop Med Int Health. 2010, 15: 1256-1265. 10.1111/j.1365-3156.2010.02601.x.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Fottrell E, Byass P: Verbal autopsy: methods in transition. Epidemiol Rev. 2010, 32: 38-55. 10.1093/epirev/mxq003.

    Article  PubMed  Google Scholar 

  12. Worldpop. [http://www.worldpop.org.uk/. Accessed: 2014-06-13. (Archived by WebCite® at http://www.webcitation.org/6QIczygsd)]

  13. Institut National de la Statistique, ICF Macro: Madagascar - Enquête Démographique et de Santé 2008–2009. 2010, Antananarivo, Madagascar: Institut National de la Statistique, 1-279.

    Google Scholar 

  14. UN OCHA ROSA: Madagascar Population per Commune (Landscan 2007). 2012, Johannesburg, South Africa: UN OCHA ROSA

    Google Scholar 

  15. Lumley T: Analysis of complex survey samples. J Stat Softw. 2004, 9: 1-19.

    Article  Google Scholar 

  16. Byass P, Chandramohan D, Clark SJ, Ambruoso LD, Fottrell E, Graham WJ, Herbst AJ, Hodgson A, Hounton S, Kahn K, Krishnan A, Leitao J, Odhiambo F, Sankoh OA, Tollman SM: Strengthening standardised interpretation of verbal autopsy data: the new InterVA-4 tool. Glob Health Action. 2012, 5: 1-8.

    PubMed  Google Scholar 

  17. Oyakhirome S, Profanter K, Kremsner PG: Assessment of fever in African children: implication for malaria trials. Am J Trop Med Hyg. 2010, 82: 215-218. 10.4269/ajtmh.2010.09-0419.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Smith DL, Guerra C, Snow RW, Hay SI: Standardizing estimates of the Plasmodium falciparum parasite rate. Malar J. 2007, 6: 131-10.1186/1475-2875-6-131.

    Article  PubMed Central  PubMed  Google Scholar 

  19. MEASURE Evaluation, MEASURE DHS, President’s Malaria Initiative, Roll Back Malaria, UNICEF, World Health Organization: Household Survey Indicators for Malaria Control. 2013, Geneva: WHO, 1-54.

    Google Scholar 

  20. Institut National de la Statistique, Programme National de Lutte contre le Paludisme, ICF International: Enquête Sur Les Indicateurs Du Paludisme À Madagascar 2013 - Rapport Préliminaire. 2013, Antananarivo, Madagascar: Institut National de la Statistique

    Google Scholar 

  21. Tusting LS, Willey B, Lucas H, Thompson J, Kafy HT, Smith R, Lindsay SW: Socioeconomic development as an intervention against malaria: a systematic review and meta-analysis. Lancet. 2013, 382: 963-972. 10.1016/S0140-6736(13)60851-X.

    Article  PubMed  Google Scholar 

  22. WHO Global Malaria Programme: Malaria Elimination: A Field Manual for Low and Moderate Endemic Countries. 2007, Geneva: WHO, 1-85.

    Google Scholar 

  23. Guillebaud J, Mahamadou A, Zamanka H, Katzelma M, Arzika I, Ibrahim ML, Eltahir EA, Labbo R, Druilhe P, Duchemin J-B, Fandeur T: Epidemiology of malaria in an area of seasonal transmission in Niger and implications for the design of a seasonal malaria chemoprevention strategy. Malar J. 2013, 12: 379-10.1186/1475-2875-12-379.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Oppenheimer SJ: Iron and its relation to immunity and infectious disease. J Nutr. 2001, 131: 616-635.

    Google Scholar 

  25. Jambou R, Tombo M-L, Raharimalala L, Rakotonjanabelo A, Rabe T, Laventure S, Boisier P: Le paludisme à Antanarivo: évaluation d’une situation post-épidémique. Cah Santé. 1998, 8: 257-264.

    CAS  Google Scholar 

  26. Ahmed SM, Hossain MS, Kabir M: Conventional or interpersonal communication: which works best in disseminating malaria information in an endemic rural Bangladeshi community?. PLoS One. 2014, 9: e90711-10.1371/journal.pone.0090711.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Roll Back Malaria: Rapports Pays N°7 - Focus Sur Madagascar. 2013, Geneva: Collection, 1-80.

    Google Scholar 

  28. Breman JG: The ears of the hippopotamus: manifestations, determinants, and estimates of the malaria burden. Am J Trop Med Hyg. 2001, 64 (1–2 Suppl): 1-11.

    CAS  PubMed  Google Scholar 

  29. Lozano R, Freeman MK, James SL, Campbell B, Lopez AD, Flaxman AD, Murray CJ: Performance of InterVA for assigning causes of death to verbal autopsies: multisite validation study using clinical diagnostic gold standards. Popul Health Metr. 2011, 9: 50-10.1186/1478-7954-9-50.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Rankin JC, Lorenz E, Neuhann F, Yé M, Sié A, Becher H, Ramroth H: Exploring the role narrative free-text plays in discrepancies between physician coding and the InterVA regarding determination of malaria as cause of death, in a malaria holo-endemic region. Malar J. 2012, 11: 51-10.1186/1475-2875-11-51.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Leitao J, Desai N, Aleksandrowicz L, Byass P, Miasnikof P, Tollman S, Alam D, Lu Y, Rathi SK, Singh A, Suraweera W, Ram F, Jha P: Comparison of physician-certified verbal autopsy with computer-coded verbal autopsy for cause of death assignment in hospitalized patients in low- and middle-income countries: systematic review. BMC Med. 2014, 12: 22-10.1186/1741-7015-12-22.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Ndila C, Bauni E, Nyirongo V, Mochamah G, Makazi A, Kosgei P, Nyutu G, Macharia A, Kapesa S, Byass P, Williams TN: Verbal autopsy as a tool for identifying children dying of sickle cell disease: a validation study conducted in Kilifi district. Kenya. BMC Med. 2014, 12: 65-10.1186/1741-7015-12-65.

    Article  PubMed  Google Scholar 

  33. Desai N, Aleksandrowicz L, Miasnikof P, Lu Y, Leitao J, Byass P, Tollman S, Mee P, Alam D, Rathi SK, Singh A, Kumar R, Ram F, Jha P: Performance of four computer-coded verbal autopsy methods for cause of death assignment compared with physician coding on 24,000 deaths in low- and middle-income countries. BMC Med. 2014, 12: 20-10.1186/1741-7015-12-20.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Murray CJL, Lozano R, Flaxman AD, Serina P, Phillips D, Stewart A, James SL, Vahdatpour A, Atkinson C, Freeman MK, Ohno SL, Black R, Ali SM, Baqui AH, Dandona L, Dantzer E, Darmstadt GL, Das V, Dhingra U, Dutta A, Fawzi W, Gómez S, Hernández B, Joshi R, Kalter HD, Kumar A, Kumar V, Lucero M, Mehta S, Neal B: Using verbal autopsy to measure causes of death: the comparative performance of existing methods. BMC Med. 2014, 12: 5-10.1186/1741-7015-12-5.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Owusu Adjah ES, Panayiotou AG: Impact of malaria related messages on insecticide-treated net (ITN) use for malaria prevention in Ghana. Malar J. 2014, 13: 123-10.1186/1475-2875-13-123.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Bowen HL: Impact of a mass media campaign on bed net use in Cameroon. Malar J. 2013, 12: 36-10.1186/1475-2875-12-36.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Koenker H, Boulay M, Ciloglu A: Rural communication activities increase net use in Tanzania. Commun Impact. 2010, 25: 1-2.

    Google Scholar 

  38. Ratsimbasoa A, Randriamanantena A, Raherinjafy R, Rasoarilalao N, Ménard D: Which malaria rapid test for Madagascar? Field and laboratory evaluation of three tests and expert microscopy of samples from suspected malaria patients in Madagascar. Am J Trop Med Hyg. 2007, 76: 481-485.

    PubMed  Google Scholar 

Download references

Acknowledgements

This research was supported by the Global Fund to Fight AIDS, Tuberculosis and Malaria via the Initiative 5% programme initiated by the French Ministry of Foreign Affairs and France Expertise Internationale (grant #12INI109). We especially thank the population of Madagascar who participated the study. We thank those who facilitated the survey, i.e., the fokonolona, heads of fokontany, local administration authorities and health authorities from Ministry of Health and National Malaria Control Programme. We thank the survey teams and members of the Monitoring Committee: P Tafangy, B Ramarosandratana, P Rajaonarison, A Finlay, L Rabarijaona, S Zigirumugabe, A Ratsimbasoa, R Andrianasolo, A Lamina, J D Rahariniaina, S A Rafalimanantsoa, E E Rambolamanana, M C Rakotoarivony, H Razafindranaivo, M Rabary, S Goedefroit, T Andrianatoandro, N Elissa, J Rakotonjato, J A Ravelonarivo, F Rakotomanana, P L Lesage, F Sandron, M Gaimard, C Bellas-Cabane, and B Meslet. Last but not least, A Boy deserves special thanks for her thorough English proofreading.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Thomas Kesteman.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CR, TK, CM, ER and DP designed the CSS questionnaire, TK, MR and VR conducted the CSS, and TK analysed the data from the CSS. CR, TK, and LR designed the CCS questionnaire, TK, LR, PP, and MR conducted the CCS, and TK and VR analysed the data from the CCS. FG designed and applied process of GIS data. CR designed the MEDALI project. TK coordinated the MEDALI project overall and drafted the manuscript. All authors conceived the study, contributed to the interpretation of data, and read and approved the final manuscript.

Electronic supplementary material

12936_2014_3655_MOESM1_ESM.xlsx

Additional file 1: Wealth component scores and measures across quintiles. Table displaying wealth component scores and measures across quintiles in the CSS. (XLSX 14 KB)

12936_2014_3655_MOESM2_ESM.xlsx

Additional file 2: Age distribution of males and females in the cross-sectional survey. Table displaying crude and weighted numbers and proportions of male, female and total according to age categories. (XLSX 11 KB)

12936_2014_3655_MOESM3_ESM.xlsx

Additional file 3: Distribution of socio-economic status quintiles in transmission patterns. Table displaying crude and weighted numbers and proportions of SES quintiles according to transmission patterns. (XLSX 12 KB)

12936_2014_3655_MOESM4_ESM.xlsx

Additional file 4: Distribution of socio-economic status quintiles and levels of education according to population density. Table displaying crude and weighted numbers and proportions of SES quintiles and levels of education according to population density. (XLSX 12 KB)

12936_2014_3655_MOESM5_ESM.xlsx

Additional file 5: Distribution of levels of education according to transmission patterns. Table displaying crude and weighted numbers and proportions of levels of education according to transmission patterns, among people ≥15 years old. (XLSX 10 KB)

12936_2014_3655_MOESM6_ESM.xlsx

Additional file 6: Age distribution of parasite rates. Table displaying weighted parasite rates (PR) in males, females and total according to age categories. (XLSX 10 KB)

12936_2014_3655_MOESM7_ESM.xlsx

Additional file 7: Distribution of long-lasting insecticidal nets, non-impregnated bed nets and no bed net users, according to definition, transmission pattern and population density. Table displaying weighted proportions of LLIN users, NIBN users and people not using bed nets, whether users are defined according to use the previous night, or every night during previous three months, according to transmission patterns and population density. (XLSX 10 KB)

12936_2014_3655_MOESM8_ESM.xlsx

Additional file 8: Sociodemographic characteristics of non-impregnated bed nets users the previous night. Table displaying weighted numbers and proportions of NIBN users the previous night, according to sociodemographic characteristics. (XLSX 13 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kesteman, T., Randrianarivelojosia, M., Mattern, C. et al. Nationwide evaluation of malaria infections, morbidity, mortality, and coverage of malaria control interventions in Madagascar. Malar J 13, 465 (2014). https://doi.org/10.1186/1475-2875-13-465

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-13-465

Keywords