In July 2020, 2283 households were randomly selected to form the base sample for the enrolment of children for the post-intervention cohort follow-up. There were 3129 eligible children, of whom 1829 were randomly selected and 1806 gave their consent, with 52.2% children under 5 years of age and a sex ratio of 1.07. In the study area, 36.7% of the surface was fine medium soil (Table 1).
Table 1: Baseline characteristics of the children enrolled in the study at the recruitment
Characteristics
|
All
|
Py-PPF-LLIN
|
Py-CFP-LLIN
|
Standard Py-LLIN
|
Number of clusters
|
60
|
20
|
20
|
20
|
Number of children
|
1806
|
604
|
601
|
601
|
Proportion of female children
|
868/1806(48.1%)
|
285/604 (47.2%)
|
291/601 (48.4%)
|
293/601(48.8%)
|
Proportion of children under 5 years
|
942/1806 (52.2%)
|
316/604 (52.3%)
|
304/601 (50.6%)
|
322/601(53.6%)
|
Proportion of <50% household with high SL)
|
32/60 (53.3%)
|
12/20 (60.0%)
|
10/20 (50.0%)
|
10/20 (50.0%)
|
Net usage the night before
|
1787/1806 (98.9%)
|
598/604 (99.0%)
|
593/601 (98.7%)
|
596/601 (99.2%)
|
Abbreviations: Py-PPF: pyrethroid-pyriproxyfen; Py-CFP: pyrethroid-chlorfenapyr; SL: standard of living
|
Between August 2020 and March 2022, there were 2112 malaria cases in cohort children in the CoZo health zone (491, 732, and 889 cases in the Py-CFP LLIN arm, Py-PPF LLIN arm, and standard LLIN arm, respectively). The highest number of malaria cases (228) occurred in August 2021 and the lowest (23) in January 2022. The monthly average malaria incidence rate was 6.65 (95% CI: 06.37-06.94) malaria cases per 100 child-month at risk for the period. The lowest malaria incidence rate of 01.41 (95% CI: 00.96-02.08) cases per 100 child-month was reached in April 2021. This may be due to the deworming treatment administered to children in March 2021. In the CFP-LLIN arm, we observed a relatively low malaria incidence rate compared with the other arms, with an average of 04.61 (95% CI: 04.22-05.03) cases per 100 child-month (Table S2-appendix).
Hotspot analysis
Malaria incidence-related maps, classified by cluster and month, were used to identify the hotspots of households most affected by malaria in the study region. We observed 13 significant hotspots (P-value < 0.05) out of a total of 113 malaria-susceptible hotspots over the whole study period (Figure 1).
Malaria incidence in the significant hotspots varied between 0.86 and 11.37 cases per 100 child-months, with a relative risk (RR) of 6.05 (P-value<0.001) and 7.78 (P-value=0.05) by the estimated risk outside the cluster, respectively (Table S1-appendix).
Stability analysis
Heterogeneity by season and by trial arm
Stability analysis of the significant clusters showed that the northern clusters were the most malaria-prone, with a malaria stability score of at least 20%, compared with some of the southern clusters, which maintained a stability score of less than 20% throughout the study period. We also observed variability in levels of stability across seasons. There was a decrease in the number of clusters with a stability score of at least 20% between the first and second rainy and dry seasons. It should also be noted that there was a change in the areas showing stability in malaria incidence between the two dry seasons (Figure 2). Similarly, northern clusters showed high stability (>50%) in malaria incidence during the two rainy seasons and the first dry season. In addition, the degree of stability in the incidence of malaria was distributed across the study arms.
Heterogeneity by season and soil type
Evaluation of the stability score by soil type shows a variation in the score depending on the soil type in the cluster. It was observed that clusters with the same soil type have almost the same characteristics in terms of malaria incidence stability. Clusters with the same geographical characteristics (soil type=medium fine) tend to have a high incidence area, regardless of the season. The opposite effect is observed in clusters geographically differentiated by fine soil type (Figure 3).
Temporal and risk factors analysis
Regarding the type of LLIN used, clusters in the Py-CFP LLIN arm had the lowest incidence rates compared to the other arms, regardless of the month. However, we observed a similar trend in the incidence rate in the other study LLIN groups (Py-PPF and standard LLIN groups) (Figure 4).
Figure 5 presents the effect of meteorological variables on the average malaria incidence. We observed a moderate positive correlation between rainfall and malaria incidence (correlation coefficient =0.58) and a moderate negative correlation with temperature (correlation coefficient =-0.56). Significant shifts related to malaria incidence were observed with temperature (p < 0.001) and rainfall (p < 0.001). We observed significant lags of up to two months for rainfall and one month for temperature.
We observed a significant association between average house altitude, soil type, temperature, and type of bed net used and malaria incidence (Table 2).
Table 2. Factors associated with malaria incidence in cohort children living in the Cove Zagnanado and Ouinhi (CoZO) health area, Benin, 2020-2022
Variables
|
Univariate analysis
|
Multivariate analysis
|
|
IRR (IC: 95%)
|
P-value
|
IRR (IC: 95%)
|
P-value
|
|
Soil type
|
Medium Fine
|
Reference
|
|
Reference
|
|
|
Fine
|
0.78 (0.61 - 1)
|
0.048
|
0.54 (0.39 - 0.75)
|
p<0.001
|
|
Standard of living
|
<50% HH have a high SL
|
Reference
|
|
Reference
|
|
|
> 50% HH have a high SL
|
0.80 (0.62 - 1.04)
|
0.094
|
0.91 (0.70 -1.20)
|
0.516
|
|
Study LLIN
|
Standard-LLIN
|
Reference
|
|
Reference
|
|
|
Py-PPF LLIN
|
0.81 (0.62 - 1.07)
|
0.15
|
0.77 (0.60 - 1)
|
0.048
|
|
Py-CFP LLIN
|
0.54 (0.42 - 0.70)
|
p<0.001
|
0.56 (0.45 - 0.71)
|
p<0.001
|
|
Temperature (°)
|
0.69 (0.66-0.73)
|
p<0.001
|
0.91 (0.84 - 0.98)
|
0.013
|
|
Rainfall (mm)
|
1.14 (1.12 - 1.15)
|
p<0.001
|
1.09 (1.05 - 1.12)
|
P<0.001
|
|
Altitude of Household
|
1 (1 -1)
|
0.6
|
1 (0.99 - 1)
|
P<0.001
|
|
Abbreviations: HH: Household; SL: Standard of living, Py-PPF: pyrethroid-pyriproxyfen; Py-CFP: pyrethroid-chlorfenapyr
|
|