Visceral leishmaniasis cyclical trends in Bihar, India - implications for the elimination programme

Modelling spatiotemporal patterns of visceral leishmaniasis incidence in two endemic states in India using environment, bioclimatic and demographic data, 2013-2022

BACKGROUND

As of 2021, the National Kala-azar Elimination Programme (NKAEP) in India has achieved visceral leishmaniasis (VL) elimination (<1 case / 10,000 population/year per block) in 625 of the 633 endemic blocks (subdistricts) in four states. The programme needs to sustain this achievement and target interventions in the remaining blocks to achieve the WHO 2030 target of VL elimination as a public health problem. An effective tool to analyse programme data and predict/ forecast the spatial and temporal trends of VL incidence, elimination threshold, and risk of resurgence will be of use to the programme management at this juncture.

METHODOLOGY/PRINCIPAL FINDINGS

We employed spatiotemporal models incorporating environment, climatic and demographic factors as covariates to describe monthly VL cases for 8-years (2013-2020) in 491 and 27 endemic and non-endemic blocks of Bihar and Jharkhand states. We fitted 37 models of spatial, temporal, and spatiotemporal interaction random effects with covariates to monthly VL cases for 6-years (2013-2018, training data) using Bayesian inference via Integrated Nested Laplace Approximation (INLA) approach. The best-fitting model was selected based on deviance information criterion (DIC) and Watanabe-Akaike Information Criterion (WAIC) and was validated with monthly cases for 2019-2020 (test data). The model could describe observed spatial and temporal patterns of VL incidence in the two states having widely differing incidence trajectories, with >93% and 99% coverage probability (proportion of observations falling inside 95% Bayesian credible interval for the predicted number of VL cases per month) during the training and testing periods. PIT (probability integral transform) histograms confirmed consistency between prediction and observation for the test period. Forecasting for 2021-2023 showed that the annual VL incidence is likely to exceed elimination threshold in 16-18 blocks in 4 districts of Jharkhand and 33-38 blocks in 10 districts of Bihar. The risk of VL in non-endemic neighbouring blocks of both Bihar and Jharkhand are less than 0.5 during the training and test periods, and for 2021-2023, the probability that the risk greater than 1 is negligible (P<0.1). Fitted model showed that VL occurrence was positively associated with mean temperature, minimum temperature, enhanced vegetation index, precipitation, and isothermality, and negatively with maximum temperature, land surface temperature, soil moisture and population density.

CONCLUSIONS/SIGNIFICANCE

The spatiotemporal model incorporating environmental, bioclimatic, and demographic factors demonstrated that the KAMIS database of the national programmme can be used for block level predictions of long-term spatial and temporal trends in VL incidence and risk of outbreak / resurgence in endemic and non-endemic settings. The database integrated with the modelling framework and a dashboard facility can facilitate such analysis and predictions. This could aid the programme to monitor progress of VL elimination at least one-year ahead, assess risk of resurgence or outbreak in post-elimination settings, and implement timely and targeted interventions or preventive measures so that the NKAEP meet the target of achieving elimination by 2030.

Molecular surveillance of insecticide resistance in Phlebotomus argentipes targeted by indoor residual spraying for visceral leishmaniasis elimination in India

Molecular surveillance of resistance is an increasingly important part of vector borne disease control programmes that utilise insecticides. The visceral leishmaniasis (VL) elimination programme in India uses indoor residual spraying (IRS) with the pyrethroid, alpha-cypermethrin to control Phlebotomus argentipes the vector of Leishmania donovani, the causative agent of VL. Prior long-term use of DDT may have selected for knockdown resistance (kdr) mutants (1014F and S) at the shared DDT and pyrethroid target site, which are common in India and can also cause pyrethroid cross-resistance. We monitored the frequency of these marker mutations over five years from 2017-2021 in sentinel sites in eight districts of north-eastern India covered by IRS. Frequencies varied markedly among the districts, though finer scale variation, among villages within districts, was limited. A pronounced and highly significant increase in resistance-associated genotypes occurred between 2017 and 2018, but with relative stability thereafter, and some reversion toward more susceptible genotypes in 2021. Analyses linked IRS with mutant frequencies suggesting an advantage to more resistant genotypes, especially when pyrethroid was under-sprayed in IRS. However, this advantage did not translate into sustained allele frequency changes over the study period, potentially because of a relatively greater net advantage under field conditions for a wild-type/mutant genotype than projected from laboratory studies and/or high costs of the most resistant genotype. Further work is required to improve calibration of each 1014 genotype with resistance, preferably using operationally relevant measures. The lack of change in resistance mechanism over the span of the study period, coupled with available bioassay data suggesting susceptibility, suggests that resistance has yet to emerge despite intensive IRS. Nevertheless, the advantage of resistance-associated genotypes with IRS and under spraying, suggest that measures to continue monitoring and improvement of spray quality are vital, and consideration of future alternatives to pyrethroids for IRS would be advisable.

Social determinants of visceral leishmaniasis elimination in Eastern Africa

Visceral leishmaniasis is a vector-borne, protozoan disease with severe public health implications. Following the successful implementation of an elimination programme in South Asia, there is now a concerted endeavour to replicate these efforts in Eastern Africa based on the five essential elimination pillars of case management, integrated vector management, effective surveillance, social mobilisation and operational research. This article highlights how key social determinants (SD) of health (poverty, sociocultural factors and gender, housing and clustering, migration and the healthcare system) operate at five different levels (socioeconomic context and position, differential exposure, differential vulnerability, differential outcomes and differential consequences). These SD should be considered within the context of increasing the success of the five-pillar elimination programme and reducing inequity in health.

Human leishmaniasis vaccines: Use cases, target population and potential global demand

The development of vaccines against one or all forms of human leishmaniasis remains hampered by a paucity of investment, at least in part resulting from the lack of well-evidenced and agreed estimates of vaccine demand. Starting from the definition of 4 main use cases (prevention of visceral leishmaniasis, prevention of cutaneous leishmaniasis, prevention of post-kala-azar dermal leishmaniasis and treatment of post-kala-azar dermal leishmaniasis), we have estimated the size of each target population, focusing on those endemic countries where incidence levels are sufficiently high to justify decisions to adopt a vaccine. We assumed a dual vaccine delivery strategy, including a wide age-range catch-up campaign before the start of routine immunisation. Vaccine characteristics and delivery parameters reflective of a target product profile and the likely duration of the clinical development effort were considered in forecasting the demand for each of the four indications. Over a period of 10 years, this demand is forecasted to range from 300-830 million doses for a vaccine preventing visceral leishmaniasis and 557-1400 million doses for a vaccine preventing cutaneous leishmaniasis under the different scenarios we simulated. In a scenario with an effective prophylactic visceral leishmaniasis vaccine, demand for use to prevent or treat post-kala-azar dermal leishmaniasis would be more limited (over the 10 years ~160,000 doses for prevention and ~7,000 doses for treatment). Demand would rise to exceed 330,000 doses, however, in the absence of an effective vaccine for visceral leishmaniasis. Because of the sizeable demand and potential for public health impact, a single-indication prophylactic vaccine for visceral or cutaneous leishmaniasis, and even more so a cross-protective prophylactic vaccine could attract the interest of commercial developers. Continuous refinement of these first-of-their kind estimates and confirmation of country willingness and ability to pay will be paramount to inform the decisions of policy makers and developers in relation to a leishmaniasis vaccine. Positive decisions can provide a much-needed contribution towards the achievement of global leishmaniasis control.

Impact of IRS: Four-years of entomological surveillance of the Indian Visceral Leishmaniases elimination programme

Background

In 2005, Bangladesh, India and Nepal agreed to eliminate visceral leishmaniasis (VL) as a public health problem. The approach to this was through improved case detection and treatment, and controlling transmission by the sand fly vector Phlebotomus argentipes, with indoor residual spraying (IRS) of insecticide. Initially, India applied DDT with stirrup pumps for IRS, however, this did not reduce transmission. After 2015 onwards, the pyrethroid alpha-cypermethrin was applied with compression pumps, and entomological surveillance was initiated in 2016.

Methods

Eight sentinel sites were established in the Indian states of Bihar, Jharkhand and West Bengal. IRS coverage was monitored by household survey, quality of insecticide application was measured by HPLC, presence and abundance of the VL vector was monitored by CDC light traps, insecticide resistance was measured with WHO diagnostic assays and case incidence was determined from the VL case register KAMIS.

Results

Complete treatment of houses with IRS increased across all sites from 57% in 2016 to 70% of houses in 2019, rising to >80% if partial house IRS coverage is included (except West Bengal). The quality of insecticide application has improved compared to previous studies, average doses of insecticide on filters papers ranged from 1.52 times the target dose of 25mg/m² alpha-cypermethrin in 2019 to 1.67 times in 2018. Resistance to DDT has continued to increase, but the vector was not resistant to carbamates, organophosphates or pyrethroids. The annual and seasonal abundance of P. argentipes declined between 2016 to 2019 with an overall infection rate of 0.03%. This was associated with a decline in VL incidence for the blocks represented by the sentinel sites from 1.16 per 10,000 population in 2016 to 0.51 per 10,000 in 2019.

Conclusion

Through effective case detection and management reducing the infection reservoirs for P. argentipes in the human population combined with IRS keeping P. argentipes abundance and infectivity low has reduced VL transmission. This combination of effective case management and vector control has now brought India within reach of the VL elimination targets.

Visceral Leishmaniasis (VL), also known as kala-azar, is a major parasitic disease in South Asia (Indian subcontinent), with 85% of the disease incidence in India. Historically VL had been controlled and almost eliminated with Indoor Residual Spraying (IRS) using dichlorodiphenyltrichloroethane (DDT). However, reinitiating this approach in 2015 failed due to high insecticide resistance in the sand fly vector and poor IRS quality, meaning that VL elimination targets were not met. To improve this the National Vector Borne Disease Control Programme changed to an effective insecticide, alpha-cypermethrin and altered the mode of application to compression pumps. Sentinel sites were established to monitor the entomological indicators, these showed the positive impact of these changes from 2016 to 2019. During this period the overall incidence of disease has decreased, and India is now on track to reach it’s target incidence for VL of less than 1/1000 people at the sub-district (block) level.

Costs and outcomes of active and passive case detection for visceral leishmaniasis (Kala-Azar) to inform elimination strategies in Bihar, India

BACKGROUND

Effective case identification strategies are fundamental to capturing the remaining visceral leishmaniasis (VL) cases in India. To inform government strategies to reach and sustain elimination benchmarks, this study presents costs of active- and passive- case detection (ACD and PCD) strategies used in India's most VL-endemic state, Bihar, with a focus on programme outcomes stratified by district-level incidence.

METHODS

Expenditure analysis was complemented by onsite micro-costing to compare the cost of PCD in hospitals alongside index case-based ACD and a combination of blanket (house-to-house) and camp ACD from January to December 2018. From the provider's perspective, a cost analysis evaluated the overall programme cost of each activity, the cost per case detected, and the cost of scaling up ACD.

RESULTS

During 2018, index case-based ACD, blanket and camp ACD, and PCD reported 1,497, 131, and 1,983 VL-positive cases at a unit cost of $522.81, $4,186.81, and $246.79, respectively. In high endemic districts, more VL cases were identified through PCD while in meso- and low-endemic districts more cases were identified through ACD. The cost of scaling up ACD to identify 3,000 additional cases ranged from $1.6-4 million, depending on the extent to which blanket and camp ACD was relied upon.

CONCLUSION

Cost per VL test conducted (rather than VL-positive case identified) may be a better metric estimating unit costs to scale up ACD in Bihar. As more VL cases were identified in meso-and low-endemic districts through ACD than PCD, health authorities in India should consider bolstering ACD in these areas. Blanket and camp ACD identified fewer cases at a higher unit cost than index case-based ACD. However, the value of detecting additional VL cases early outweighs long-term costs for reaching and sustaining VL elimination benchmarks in India.

Visceral leishmaniasis: Spatiotemporal heterogeneity and drivers underlying the hotspots in Muzaffarpur, Bihar, India

BACKGROUND

Despite the overall decrease in visceral leishmaniasis (VL) incidence on the Indian subcontinent, there remain spatiotemporal clusters or 'hotspots' of new cases. The characteristics of these hotspots, underlying transmission dynamics, and their importance for shaping control strategies are not yet fully understood and are investigated in this study for a VL endemic area of ~100,000 inhabitants in Bihar, India between 2007-2015.

METHODOLOGY/PRINCIPAL FINDINGS

VL incidence (cases/10,000/year) dropped from 12.3 in 2007 to 0.9 in 2015, which is just below the World Health Organizations' threshold for elimination as a public health problem. Clustering of VL was assessed between subvillages (hamlets), using multiple geospatial and (spatio)temporal autocorrelation and hotspot analyses. One to three hotspots were identified each year, often persisting for 1-5 successive years with a modal radius of ~500m. The relative risk of having VL was 5-86 times higher for inhabitants of hotspots, compared to those living outside hotspots. Hotspots harbour significantly more households from the two lowest asset quintiles (as proxy for socio-economic status). Overall, children and young adelescents (5-14 years) have the highest risk for VL, but within hotspots and at the start of outbreaks, older age groups (35+ years) show a comparable high risk.

CONCLUSIONS/SIGNIFICANCE

This study demonstrates significant spatiotemporal heterogeneity in VL incidence at subdistrict level. The association between poverty and hotspots confirms that VL is a disease of 'the poorest of the poor' and age patterns suggest a potential role of waning immunity as underlying driver of hotspots. The recommended insecticide spraying radius of 500m around detected VL cases corresponds to the modal hotspot radius found in this study. Additional data on immunity and asymptomatic infection, and the development of spatiotemporally explicit transmission models that simulate hotspot dynamics and predict the impact of interventions at the smaller geographical scale will be crucial tools in sustaining elimination.

Visceral leishmaniasis cyclical trends in Bihar, India - implications for the elimination programme

Background: Visceral leishmaniasis (VL) is a vector-borne disease of public health importance in India, with the highest burden of disease in the states of Bihar, Jharkhand, West Bengal and Uttar Pradesh. The disease is currently targeted for elimination (annual incidence to less than one per 10,000 population) using indoor residual spraying, active case detection and treatment. Historically the disease trend in India has been regarded as cyclical with case resurgence characteristically occurring every 15 years.  Understanding this pattern is essential if the VL elimination gains are to be sustained. To better understand the cyclical trends, annual climatic indicators including rainfall, temperature and humidity over time were compared with annual VL case incidence data. 

Methods: Annual climate data (rainfall, average and maximum temperature and specific humidity) from 1956-2004 were used to identify potential factors influencing VL incidence.  Months relevant to the VL life-cycle were identified and defined (Monsoon, Sand-fly Peak, Pre-Sand-fly Peak and Annual) for analysis. The Kruskall-Wallis test was used to determine significant difference between categorical rainfall and VL incidence, whilst univariate negative binomial regression models were used to determine predictors of disease incidence.

Results: The negative binomial regression model showed statistically significant associations (p <0.05) for VL incidence and maximum temperature, and average temperature, when considering annual and pre-sand fly peak time periods. No other associations between humidity, rainfall or temperature and VL incidence were detected (all values p >0.05). 

Conclusion: The VL programme in Bihar has made significant progress in adopting best practices for improved treatment and vector control, with the aim to achieve VL elimination.  However, open access granular programme data for indoor residual spray activities and case detection is required to fully understand the role of climate in disease transmission and potential resurgence.