Skip to content
ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Article

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

[version 1; peer review: 1 approved, 2 approved with reservations]
PUBLISHED 21 Feb 2018
Author details Author details

Abstract

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.

Keywords

Visceral leishmaniasis, kala-azar, climate, India, Bihar, elimination

Background

Visceral leishmaniasis (VL) is a tropical disease of public health importance in India. Caused by Leishmania donovani parasites and transmitted by the sand fly Phlebotomus argentipes, the disease is anthroponotic in India and endemic in four States: Bihar, Jharkhand, West Bengal and Uttar Pradesh (Addy & Nandy, 1992). It is estimated that 90% of VL cases in India originate from Bihar (Singh et al., 2010). In 2005, India with Bangladesh and Nepal, set a target to eliminate VL with incidence at less than one case per 10,000 population by 2015, which was subsequently readjusted to 2017. This is in line with the London declaration on Neglected Tropical Diseases (NTDs), which set the aim to eliminate VL as a public health problem by 2020.

The seasonal and interannual scale climate trends seen in vector-transmitted diseases have been widely reported (Connor et al., 1999). In India, malaria prediction was first conducted using climatic and socioeconomic data in the 1910s by Captain S. R. Christophers from the British Army (Gill, 1938). This system was implemented until the 1940s, when malaria was not seen as disease of public health concern within the Indian subcontinent (Gill, 1938). Despite the burden of VL disease in India, little work has been done to associate climatic indicators with case incidence. Historical trends have previously shown that there is a resurgence of VL every 15 years post control: whilst widely accepted, there is limited understanding about the cause of this pattern (Malaviya et al., 2011; Muniaraj, 2014).

Napier first proposed that VL in India was cyclical in nature in 1946, before the suspected epidemics of 1977 and 1991. Visceral leishmaniasis epidemics in Assam (1875 and 1950) however were thought to be caused by the influenza pandemic of 1918–1919, malaria, famine and earthquakes (McCombie Young, 1924; Napier, 1943; Rogers, 1908). Findings from McCombie Young (1924) in Assam, supported this and after the introduction of antimonial drugs, a dramatic decrease in case mortality from 90% to 10% was observed (Dye, 1992; McCombie Young, 1924; Muniaraj, 2014). More recently, Muniaraj attempted to elucidate the cause behind the 10–15 year cycle, attributing the VL case trend to vector control practices and the need for appropriate therapeutic strategies to reduce VL mortality (Muniaraj, 2014). However, Dye has referred to this cyclical trend of L. donovani VL as “part of the folk wisdom of tropical medicine”, suggesting that post-kala-azar dermal leishmaniasis (PKDL) patients serve as a parasite reservoir: chronically infective and available to sand flies (Dye, 1992; Rogers, 1908). Studies by Addy and Nandy have suggested that disease transmission during the suspected epidemic period of 1977 in Bihar, was most likely passed on person to person from active VL cases (Addy & Nandy, 1992).

Associations between climatic indicators, vector density and VL incidence in India have previously been detected when analysing district level data (Bern et al., 2010; Dhimal et al., 2015; Malaviya et al., 2011; Tiwary et al., 2013). Generalised associations for Bihar as a State, currently do not exist. Humidity, temperature, rainfall, soil temperature and moisture are all widely accepted factors that influence the development of Ph. argentipes (Bhunia et al., 2010). There are two seasonal peaks of Ph. argentipes in Bihar, one between March to June and another in October to November (Dinesh et al., 2001; Picado et al., 2010). Associations have also been made between VL incidence and air temperature, relative humidity and annual rainfall in the Gangetic plain, which includes regions within Bihar (Bhunia et al., 2010). Malaviya et al. noted seasonality of VL case trends, with peaks in March to April and a minor secondary peak in July (Malaviya et al., 2011). In addition, the onset of disease symptoms were recorded, in descending order, from April to June, June to September and October to December (Perry et al., 2013). This paper investigates the impact of annual climatic and disease trends on the cyclic nature of VL burden in Bihar.

Methods

Study area

The State of Bihar is located in the eastern part of India, extending over 94,163 km2 and averages 52.73m above sea level. It is entirely land-locked, bounded by Nepal in the north and the State of Jharkhand in the south. To the east lies the humid State of West Bengal and to the west, the sub-humid State of Uttar Pradesh. There are 38 districts in Bihar and the total population, according to the last census in 2011 census, was 104,099,452.

Data sources

For full details on data sources accessed for this research study, see supplementary Table 1.

Epidemiological data

Yearly case numbers from Bihar were compiled using open access data available from the National Vector Borne Disease Control Programme (NVBDCP), Government of India (NVBDCP n.d; NVBDCP n.d.b; Kala-azar or Visceral Leishmaniasis) and published literature quoting historical case data from Government of Bihar (GOB) and Government of India records. Case numbers obtained from the NVBDCP originate from primary health care facilities and hospital records within Bihar, where diagnostic testing, followed national guidelines (NVBDCP n.d.). Annual incidence was calculated using population totals obtained from Bihar population data, obtained from decennial census surveys conducted by Ministry of Home Affairs, Government of India (Directorate of Economics & Statistics, 2011). Annual population growth statistics from the Directorate of Economics and Statistics, Patna, Bihar (1956–1959) and the World Bank Databank (1961–2014) was used to produce yearly population figures between census surveys ( by multiplication of annual growth rate and population for previous year). Census numbers for Bihar were taken directly where available. Visceral leishmaniasis incidence was calculated per 10,000 for all years where case and total population data were available.

Indoor residual spraying (IRS) data

Peer-reviewed literature searches were used to determine years when IRS was performed in Bihar between 1931–2014 (Date of last search: 31 November 2015). The search term strings included “visceral leishmaniasis”, “kala-azar”, “indoor residual spraying”, “IRS”, “India” and “Bihar” were used on PubMed and Google search engines to establish historical spray activities. Bibliographies from relevant peer-reviewed journals were used to identify additional information and data.

Rainfall data

Bihar State historical monthly rainfall data (1871–2013) was obtained at the time of analysis from the Indian Institute of Tropical Meteorology. Adapted Earth System Science Organisation (ESSO)-India Meteorological Department (MoES) long period average (LPA) rain fall categories were used to classify rainfall data for the selected time period into standardised groups: Deficient = rainfall below 80% LPA, Below Normal = rainfall 80–90% LPA, Normal = rainfall 90–110% LPA, Above Normal = rainfall 110–120% LPA and Excess = rainfall above 120%.

Temperature and humidity

The average and maximum temperature and specific humidity data corresponding to the coordinates of Bihar were extracted from satellite sources and digital databases through the International Research Institute for Climate and Society website. Data for Bihar, accessed through the website, originated from the National Centers for Environmental Prediction (NCEP) and the National Centre for Atmospheric Research (NCAR), as part of the NCEP/NCAR Climate Data Assimilation System (CDAS) Reanalysis Project. Daily data provided at a 1 degree spatial resolution were aggregated to provide a monthly output for the entire State. Temperature data was converted from Kelvin to degrees Celsius.

Climate periodicity

In order to aggregate monthly climate data into an appropriate annual indicator, literature searches were conducted and two clear time periods were identified as the most relevant for VL transmission: Monsoon Season (June-September), and Sand Fly Peak abundance (spanning months March - June and October - November) (Dinesh et al., 2001).

As the sand fly life cycle is estimated to take approximately one month, a one month window prior to sand fly abundance peaks (spanning months February-May and September-October), categorised as Pre-Sand Fly Peak, was also considered relevant in determining risk factors associated with VL transmission.

Finally, the Annual category for all variables was used to identify any trends between annual VL incidences and calendar year annual climatic variables.

Analysis

Monthly data from the specific climatic periods (Annual, Monsoon Season, Sand Fly Peak and Pre-Sand Fly Peak) were averaged to produce a single annual figure per time period. A non-parametric Kruskal-Wallis Test was performed to determine whether VL incidence varied between the different annual rain categories (deficient, below normal, normal, above normal, excess). Univariate negative binomial regression models were fitted to the VL incidence data for each of the climatic indicators (specific humidity, average and maximum temperature) using data from the periods 1956–1960 and 1977–2013. To account for multiple testing, the resulting p-values were adjusted using the false discovery rate (FDR) correction (Benjamini & Hochberg, 1995). Climate variables with a p-value less than 5% were significant, plus the Akaike Information Criterion (AIC) values were compared to determine which of the models was the best fit to the data. To account for long-term changes in climate and VL incidence, the analysis was repeated using temporally continuous data only (1977–2013) using climate anomalies (observed values – period average) and a linear temporal trend as predictors in the negative binomial regression models. Analysis was conducted in R 3.2.5, SPSS (Version 22) and Graphpad Prism (Version 6.07).

Results

Epidemiological data

The earliest record of case data available for Bihar was in 1934, however review of the historical literature also showed a disparity in the total number of cases reported between 1934 and 2014. A total of 1,190,166 cases were recorded when using GOB only records(NVBDCP n.d; NVBDCP n.d.b; Sanyal et al., 1979). Conversely, 1,461,963 cases were recorded when including GOB records and independent surveys (GOB&I) data conducted through National Institute of Communicable Diseases, Delhi (1977), and Government of India under the Malaria Department (1991).

Records of VL cases were unavailable between 1938–1955 and 1961–1976, most likely coinciding with public health priorities moving away from VL in Bihar. However, case records were available consistently from 1977 to present (Figure 1). As noted by Bora, “due to limited surveillance strategies” information on outbreaks before 1977–78 was patchy (Bora, 1999). Cases recorded were largely confirmed by tissue specimen microscopy, however since 2005 the NVBDCP guidelines advise use of the rk39 rapid diagnostic test. Other methods such as direct agglutination test, enzyme-linked immunosorbent assay, parasite culture and PCR may have been used at a smaller scale to diagnose cases presented in health facilities associated with research institutes.

0ae3d92f-5ddc-4bb3-8628-f840a9bdc524_figure1.gif

Figure 1. Historical Overview: Annual VL Incidence in Bihar and IRS activities.

The dark blue bars represent the annual VL incidence per 10,000, as reported by the VL Programme (GOB only). In the absence of historical population at risk data, the total population in Bihar was used to calculate incidence. Independent surveys conducted in 1977 by National Institute of Communicable Diseases, Delhi and GOB and 1991 by Government of India under the Malaria Department, reported VL incidence per 10,000 is shown in orange for those two years (GOB&I). Light blue shading represents years where IRS was ongoing, whilst white represents years where IRS was stopped. Years marked with an asterisk were used in analysis of climatic, case and IRS indicators.

GOB VL incidence data, collected solely through Primary Health Centres, showed the highest VL incidence was in 1992 (11.480 per 10,000 (75,523)) (Figure 1 & Table 3). The second highest VL incidence was recorded in 1991 (59614 cases (9.238 per 10,000). Bora reported that a further peak in VL cases occurred in 1974, however, as the programme assumed elimination had been achieved and stopped VL case data collection, no information was available to calculate incidence (Bora, 1999). The lowest incidences recorded were seen in 2014 (7615 cases - 0.705 per 10,000) and 2013 (1.005 per 10,000 (10730 cases) (Figure 1 & Table 3).

According to GOB&I data, the highest annual case reports were recorded in 1977 (100,000 cases – incidence of 20.720 per 10,000) and 1991 (250,000 cases – incidence of 38.741 per 10,000) (Figure 1). GOB data only was used in the analysis of climatic indicators.

IRS data

IRS was first adopted for VL elimination in Bihar through the National Malaria Eradication Programme (NMEP) in the 1940s (Thakur, 1984a). Activities were stopped 1962–76 in response to a rapid fall in VL case numbers, however resumed for two years again in 1977, with support from Government of India, to combat the outbreak declared that year (Muniaraj, 2014; Thakur, 1984a). Following a survey in 1991, an outbreak of VL was declared and IRS was conducted 1992–1995, after which IRS was stopped (Thakur, 2007). In 2005, with the launch of the “Kala-azar Elimination Programme”, IRS was restarted to reach the elimination target (less than one in 10,000 within a primary healthcare centre/block) by a recently revised deadline of 2020 (Mondal et al., 2009; WHO, 2016). No records of IRS activities were available during years 1934–1936, when there was high case burden, and for the VL peak of 1974 (Bora, 1999).

Rainfall

Sixty-two years of rainfall data, 1951–2013, were available to determine trends and association between rain time categories and VL incidence, however due to gaps in VL incidence data (1951–1955 and 1961–1976), a total of 42 matched years (1956–1960 and 1977–2013) were used to complete the analysis (marked with an * in Figure 1). The average (Annual) rainfall for the total 42 matched years was 1017.66mm. The average rainfall during Monsoon Season periods for the same matched years was 2517.92mm. When considering the sand fly-related time categories for these years, the average rainfall during Sand Fly Peak months was 586.46mm and 659.01mm during the Pre-Sand Fly Peak months.

After categorising the data into the five standardised groups, the rainfall during Annual and Monsoon Season time periods were found to be mostly Normal to the LPA (18/42 (42.86%) and 14/42 (33.33%) respectively). During the Sand Fly Peak period, rainfall was classified most frequently as Deficient to the LPA (13/42 (30.91%)). Conversely, when considering the Pre-Sand Fly Peak period, rainfall was found to be Excess to the LPA most frequently (12/42 (28.57%)).

Temperature and humidity

A total of 66 years of temperature and humidity data was available from data sources, however, given the limited availability of VL incidence data, a total of 43 years of temperature and specific humidity data (1956–1960 and 1977–2014) were used. As shown in Table 2, the mean average temperature over the different time groupings ranged between 21.74 and 25.36°C whilst the mean maximum temperature was 28.36–31.14°C. Furthermore, the mean specific humidity (x1000) ranged between 12.000 – 20.000kg/kg.

Climatic variables and VL incidence

Rainfall during 1992, the highest VL incidence according to GOB data, was classified as Deficient during all time periods (Table 3). During the year of the second highest VL incidence, rainfall was classified as Below Normal to the LPA when considering the Annual and Monsoon Season time periods (Table 3). During Sand Fly Peak and Pre-Sand Fly Peak time periods, rainfall was Deficient to the LPA.

The lowest VL incidence was recorded in 2014, however no rainfall data was available for this time frame. VL incidence was recorded as low in 2013 (1.005 per 10,000) during which, rainfall was classified as Deficient during the Monsoon Season, Excess during the Pre-Sand Fly Peak and Sand Fly Peak, and Normal for the Annual (Table 4).

When considering VL incidence and rainfall together, no statistically significant relationship was detected for any of the time periods (Annual (p=0.265), Monsoon Season (p=0.281), Sand Fly Peaks (p=0.602) and Pre-Sand Fly Peak (p=0.416)) (Table 4). Overall, average VL case incidence, irrespective of time period, was lower (2.172–3.133) during years of higher rainfall (Above Normal or Excess to the LPA).

The average VL incidence when considering all 22 IRS years, between 1956 and 2014 (Table 1), was 3.590. For the 21 years when IRS stopped, the average VL incidence was 3.000. The negative binomial regression model fitting identified statistically significant negative associations between both Average Temperature and Maximum Temperature and VL incidence (GOB data) during the Annual (RR= 0.682, p=0.0186 and RR=0.667, p=0.003 respectively) and Pre-Sand fly Peak (RR=0.796, p=0.0494 and RR= 0.814, p=0.0501 (close to statistically significant) respectively) time periods only (Table 2). However, after correcting for multiple testing using the FDR correction, only Annual Maximum temperature was significant (p=0.0360). Of these regression models, Annual Maximum temperature was the best fitting model as assessed using AIC. No other significant associations were detected for temperature and other time periods, or humidity and VL incidence (Table 2). Summaries of the models fitted to climate anomalies for the period 1977–2013 can be found in Supplementary Table 2. After adjusting for a decreasing temporal trend in incidence, only annual specific humidity anomalies were associated with incidence (p=0.0398) however after adjusting for multiple testing this association became non-significant (p=0.4776).

Table 1. Overview of IRS activities for VL in Bihar, India.

IRSYearsProgramme NameProgramme led by
11953–1962National Malaria
Eradication
Programme
National-run
Programme
21977–1979Kala-azar Control
Programme
National Institute
of Communicable
Diseases with UNDP
assistance
31992–1995Kala-azar Control
Programme
42005–PresentKala-azar Elimination
Programme
State-run Programme

Table 2. Temperature and Humidity descriptive for all time periods.

Negative Binomial model with VL incidence using data from GoB only.

AnnualMonsoon
Season
Sand Fly PeakPre-Sandfly
Peak
Average TemperatureMean (°C )21.7425.3623.7322.99
Range (°C )20.88-23.2524.05-26.7722.85-25.7821.95-24.60
s.d.0.5510.4460.6840.782
Coefficient0.68180.87080.83830.7959
95% CI0.5040 - 0.93550.5127 - 1.49190.6541 - 1.09140.6454 - 0.9890
p-value0.0186*0.51420.17880.0494*
Adjusted p-value0.11160.68560.35760.1503
AIC927932931929
Maximum TemperatureMean (°C )28.3628.7531.1430.18
Range (°C )27.30-29.7527.21-30.8129.89-32.8428.78-31.97
s.d.0.6450.6450.7960.866
Coefficient0.6670.94660.8650.8137
95% CI0.5055 - 0.88630.6656 - 1.35740.6793 - 1.10860.6660 - 0.9977
p-value0.003*0.7090.21740.0501¤
Adjusted p-value0.0360*0.85080.37270.1503
AIC926933931929
Specific Humidity (x1000)Mean (kg/kg )13131313
Range (kg/kg)12.000 - 13.00018.800 - 20.00011.000 - 13.00010.000 - 13.000
s.d.0.000330.000310.000560.00053
Coefficient0.66640.96640.98450.8296
95% CI0.4121 - 1.03130.6600 - 1.40130.7431 - 1.30110.5784 - 1.1716
p-value0.08130.85860.91650.2811
Adjusted p-value0.19510.91650.19650.4217
AIC930933933932

* P-value is significant at the 0.05 level. ¤P-value is nearly significant. (Akaike Information Criterion (AIC), Confidence interval (CI), standard deviation (s.d.))

Table 3. Rainfall, VL case incidence per 10,000 (GOB data) and Kruskal-Wallis analysis.

(Average temperature (Ave. Temp), Maximum temperature (Max. Temp), Specific Humidity (Specific Hum.))

YearVL
incidence
AnnualMonsoon SeasonSand Fly PeaksPre-Sand fly Peaks
Ave.
Temp
(°C)
Max.
Temp
(°C)
Specific
Hum.
(kg/kg)
RainfallAve.
Temp
(°C)
Max.
Temp
(°C)
Specific
Hum.
(kg/kg)
RainfallAve.
Temp
(°C)
Max.
Temp
(°C)
Specific
Hum.
(kg/kg)
RainfallAve.
Temp
(°C)
Max.
Temp
(°C)
Specific
Hum.
(kg/kg)
Rainfall
19919.23821.6428.340.013Below
Normal
25.6429.090.020Below
Normal
23.6231.470.012Deficient23.0530.590.012Deficient
199211.48021.1227.970.012Deficient25.3929.230.019Deficient23.4431.290.011Deficient22.0929.750.011Deficient
20131.00522.3829.380.013Normal22.3828.930.020Deficient24.6332.310.012Excess24.4231.920.012Excess
20140.70522.4629.750.013N/A22.4630.810.020N/A24.5632.830.011N/A23.5531.170.012N/A

Table 4. Annual and Maximum Temperature, Humidity and Rainfall during high and low VL incidence.

VL incidence based on data from GOB data. (Minimum (Min.), Maximum (Max.), Average (Ave.)

AnnualMonsoonSand fly PeaksPre Sand fly Peaks
Rainfall
Category
DescriptionRainfall
Range/
mm
#
Years
Incidence per 10,000P
value
Rainfall
Range/
mm
#
Years
Incidence per 10,000P
value
Rainfall
Range/
mm
#
Years
Incidence per 10,000P
value
Rainfall
Range/
mm
#
Years
Incidence per 10,000P
value
Min.Max.Ave.Min.Max.Ave.Min.Max.Ave.Min.Max.Ave.
Deficient<80% LPA670.250-
792.250
71.14811.4803.3050.2651663.750-
2010.250
91.00511.4803.2910.281255.333-
454.167
131.14811.4803.7040.602265.000-
524.166
91.70211.4804.4620.416
Below
Normal
80–90%
LPA
825.333-
887.500
51.7029.2384.3512038.750-
2237.250
71.2339.2383.492480.833-
524.167
61.6988.5233.489532.166-
574.000
91.1488.5233.193
Normal90–110%
LPA
926.833-
1089.916
181.0058.5233.7002340.500-
2764.500
141.1818.5234.034543.500-
627.833
72.0566.5873.672598.833-
720.000
72.1404.7982.961
Above
Normal
110-120%
LPA
1124.333-
1198.585
61.6305.9672.8892780.750-
2848.750
51.6235.9673.033660.833-
697.333
61.9914.7982.970737.000-
781.833
51.5606.5872.923
Excess>120% LPA1253.750-
1441.083
61.5603.9822.1723085.000-
3960.750
71.5603.1822.304723.000-
960.667
101.0058.4852.924792.833-
1029.333
111.0058.4853.133

Discussion

Associations between climatic indicators and vector-borne diseases have been established within many diseases, including VL (Chen et al., 2012; Dhimal et al., 2015; Malaviya et al., 2011; Perry et al., 2013; Thomson et al., 1999; Thomson et al., 2005; Tiwary et al., 2013). Predictors of Leishmaniasis spp. transmission cycles and sensitivity to meteorological and climatic variables is known to vary spatially; dependent on a range of factors including species composition, host competence, contact rates, vector competence, sensitivity to weather and other environmental stressors (Chaves et al., 2008; Lewnard et al., 2014). These factors however also remain essential for developing early warning systems to prevent epidemics. Due to the lack of spatial and temporal open access granular case data, there are a limited number of studies investigating the interaction of climatic variables and VL in Bihar. Previous studies have predominantly focused on district level climatic trends identifying various monthly or annual temperature, relative humidity and rainfall measurements as the key variables affecting VL transmission (Malaviya et al., 2011; Napier, 1926; Sivaramakrishnaiah & Ramanatham, 1967; Tiwary et al., 2013). This is the first study to look at historical data spanning 43 years and include climatic variables during specific month periods estimated to impact transmission, namely seasonal sand fly abundance and monsoon period (Picado et al., 2010).

In this study, four statistically significant negative associations were detected between Average Temperature and VL incidence (RR= 0.682, p=0.0186 (Annual) and RR=0.796, p=0.0494 (Pre-Sand fly Peak)) and Maximum Temperature and VL incidence (RR=0.667, p=0.003 (Annual) and RR= 0.814, p=0.0501 (near significant - Pre-Sand fly Peak), using GOB data (Table 2). However after considering the FD correction for multiple testing, only Maximum Temperature and VL incidence retained a significant association (0.0360). Overall this shows that risk is reduced with each unit increase of temperature. This could be due to human behavioural patterns, such as more people sleeping outside in warmer temperatures, unfavourable temperatures for optimal sand fly emergence or other trends associated with the vector and transmission dynamics. This is the first time that an association with temperature has been suggested to VL incidence for the Pre-Sand fly Peak time period. Annual humidity and VL incidence when fitted to climate anomalies (1977–2013) and adjusted for a decreasing temporal trend in incidence, also showed some association (p=0.0398) however only when corrections for multiple testing were not considered. This suggests that Annual Humidity could also play a similar key role in VL transmission to Annual Maximum Temperature: influencing human sleeping behaviour.

Historical research indicated that the factors favouring the development of an epidemic would include: altitude below 609.60m above sea level, greater than 1270mm annual rainfall, greater than 70% mean humidity, presence of alluvial soil, maximum temperature below 37.78°C, with diurnal variation less than -6.666°C, abundant vegetation and subsoil water, and a rural environment (Napier, 1926; Sivaramakrishnaiah & Ramanatham, 1967). Bihar is estimated to be 52.73m above sea level, and maximum temperature was below 37.78°C across all time periods (Table 2), however rainfall during years of high VL incidence was considerably lower than previously suggested (Table 4). Other criteria such as “abundant vegetation and subsoil water” cannot be quantified to form a comparable variable within this study. Given the changes in global climate, these historical guidelines on optimal conditions for a VL epidemic may no longer be suitable.

More recent research conducted by Malaviya et al. (Malaviya et al., 2011) showed associations between monthly VL incidence, rainfall and specific humidity, but contrary to this study, found no clear association between incidence and average or maximum temperature. A study conducted in the Gangetic plain, which includes VL endemic West Bengal and Bangladesh, suggested optimal meteorological factors for VL incidence included air temperature of 25.0–27.5°C, relative humidity of 66–75% and annual rainfall between 1000 and 1600 mm (Bhunia et al., 2010). Due to a lack of available open access data, only annual rainfall data can be compared to findings by Bhunia et al. (Bhunia et al., 2010): rainfall during the periods of high VL incidence fell below the suggested optimal range (1992: 834.667mm, 1991: 670.250mm).

The term epidemic is generally referred to as a sudden increase in the number of cases, beyond what is normally expected within the population of a specific area (CDC, 2016). The underlying assumption when identifying an epidemic is that the data collection methods adopted are comparable in accuracy and precision. When including GOB&I data, the highest VL incidences were in 1977 and 1991, and historically they have been referred to as “epidemics” (Bora, 1999; Thakur, 1984b; Thakur, 2007; Thakur et al., 2008). Data forming the basis for these suspected epidemics were collected through active case detection surveys. Such surveys typically detect higher levels of incidence than those recorded through normal passive reporting channels and highlight the need for better case reporting mechanisms. Comparisons between GOB and GOB&I data suggest that case numbers recorded through standard reporting channels in 1977 and 1991 were under-reported (81.41% and 76.15% respectively). Under-reporting of VL in Bihar has been documented since 1977, with most recent reports up to 2006 (Bora et al., 1994; Singh et al., 2010; Thakur, 1984b). In 2010, it was reported that a “substantial proportion” of patients were visiting private laboratories for diagnosis, which could explain underreporting documented for VL (Hasker et al., 2010). The quality of passive case reporting may have fluctuated greatly over the years, so a single correction factor cannot be applied across all years of data without further evidence. While the question remains if these are true epidemics, both surveys did result in the implementation of IRS programmes to control the disease.

When considering the current data set, given that the peaks originate from independent surveys, they may simply be artefacts highlighting the true burden of disease. As reported by Bora et al., case data recorded by GOB typically includes cases recorded through government health facilities only as private practitioners treating for VL fail to report cases they have treated to the health authorities (Bora et al., 1994). This was also the issue faced by Dye and Wolpert when modelling the Assam epidemics of 1875 and 1950; where poor quality data limited the findings for the cause behind the second epidemic (Dye & Wolpert, 1988). Given the discrepancies in data and sources, further evidence is required to confirm the presence of a VL disease cycle.

The anthroponotic nature of the disease in India, suggests records of PKDL could provide additional evidence to explain the peaks seen in 1977 and 1991 (Addy & Nandy, 1992). Open access case data sources used within this analysis, provided the number of deaths and cases only, restricting further exploration using this variable.

Historically VL vector monitoring within Bihar has been sporadic and limited by lack of spatial and temporal granularity (Dinesh et al., 2001; Malaviya et al., 2011; Picado et al., 2010; Tiwary et al., 2013). This has prevented the robust multi-annual evaluation of P. argentipes abundance peaks, physiology status and transmission effectiveness. For effective retrospective analysis to understand disease incidence, understanding the vector is crucial, furthermore climatic parameters can be further refined based on vector behaviour.

Time periods used within this study were identified as periods of high sand fly activity or adverse weather. The months included within the Sand Fly Peak time period (March - June and October - November) have been shown previously to coincide with months of high case numbers (Dinesh et al., 2001). This relationship needs further exploration, through spatial models, if monthly granular data were available, to further understand disease trends and outbreak prediction.

The available published data describes Bihar’s primary intervention against VL, IRS, as a dichotomous variable, limiting the understanding of its impact on the natural disease patterns. The VL incidence between IRS and non-IRS years is very similar when considering the overall dataset and interestingly increases (3.590) during years of IRS in comparison to years of non-IRS (3.000), suggesting that associations are potentially being masked by the lack of detail about the intervention. Other IRS programme indicators, such as spray quality, would also supplement understanding the success of IRS round, allowing for a robust statistical model to be developed (Coleman et al., 2015).

Since the start of the VL control efforts strategies adopted for case finding have changed, with the most concerted efforts ongoing considering the current elimination target. Accredited Social Health Activists (ASHAs) today play a key role in promoting healthcare activities such as immunisation, reproductive and child health and detection of VL cases (Das et al., 2016). Research has shown this incentive based approach for identifying VL patients led to an improved referral rate, translating into early detection of cases and a long-term drop in case numbers (Das et al., 2016).

The Indian VL elimination programme have taken significant steps to improve case data quality over recent years: increasing VL awareness, providing incentives to local health workers (ASHAs) and implementing district level tracking systems. For the VL elimination programme this means that to sustain the gains made a system such as reactive IRS as seen in other diseases must be established, or, the levels of VL will increase again as observed in Figure 1 (Coleman et al., 2008). Using the current WHO standards for initiating alerts in India, as in Brazil, after months when the incidence has been double the monthly average, is likely to be a sub-optimal criterion with limited infrastructure and funding to ensure rapid response: particularly after the elimination target has been achieved and the disease is no longer considered a public health concern (Coleman et al., 2008; Lewnard et al., 2014). A comprehensive understanding of the causative factors for disease trends is required to develop a suitable outbreak detection and response system within the VL Indian context.

Declarations

Data availability

Data underlying this study is available from Open Science Framework: http://doi.org/10.17605/OSF.IO/2T9DX (Deb et al., 2018)

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 21 Feb 2018
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
Gates Open Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Deb RM, Stanton MC, Foster GM et al. Visceral leishmaniasis cyclical trends in Bihar, India – implications for the elimination programme. [version 1; peer review: 1 approved, 2 approved with reservations]. Gates Open Res 2018, 2:10 (https://doi.org/10.12688/gatesopenres.12793.1)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 21 Feb 2018
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions

Are you a Gates-funded researcher?

If you are a previous or current Gates grant holder, sign up for information about developments, publishing and publications from Gates Open Research.

You must provide your first name
You must provide your last name
You must provide a valid email address
You must provide an institution.

Thank you!

We'll keep you updated on any major new updates to Gates Open Research

Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.