Impact of cattle on joint dynamics and disease burden of Japanese encephalitis and leptospirosis

Japanese encephalitis (JE) is a mosquito-borne neglected tropical disease. JE is mostly found in rural areas where people usually keep cattle at home for their needs. Cattle in households reduce JE virus infections since they distract vectors and act as a dead-end host for the virus. However, the presence of cattle introduces risk of leptospirosis infections in humans. Leptospirosis is a bacterial disease that spreads through direct or indirect contact of urine of the infected cattle. Thus, cattle have both positive and negative impacts on human disease burden. This study uses a mathematical model to study the joint dynamics of these two diseases in the presence of cattle and to identify the net impact of cattle on the annual disease burden in JE-prevalent areas. Analysis indicates that the presence of cattle helps to reduce the overall disease burden in JE-prevalent areas. However, this reduction is dominated by the vector’s feeding pattern. To the best of our knowledge, this is the ﬁrst study to examine the joint dynamics of JE and leptospirosis.


Introduction
Japanese encephalitis viral disease (JE) was first documented in 1871 in Japan [1].Japanese encephalitis (JE) is the main cause of viral encephalitis in many countries of Asia and the western Pacific.24 countries in the WHO South-East Asia and Western Pacific regions have Japanese encephalitis virus (JEV) transmission risk, which includes more than 3 billion people [1,2].JEV is transmitted to humans through bites from infected mosquitoes of the Culex species (mainly Culex tritaeniorhynchus).The virus exists in a transmission cycle between mosquitoes, pigs and/or water birds and is transmitted to humans by infected mosquito bite.JE is predominantly found in rural and periurban settings [3,4].Each year there are nearly 68,000 clinical cases of JE globally, with approximately 13,600 to 20,400 deaths.JE primarily affects children.Most adults in endemic countries have natural immunity after childhood infection, but individuals of any age may be affected.Most people infected with JE do not have symptoms or have only mild symptoms.However, a small percentage of infected people develop inflammation of the brain (encephalitis), with symptoms including sudden onset of headache, high fever, disorientation, coma, tremors, and convulsions.Approximately 1 in 250 infections results in severe clinical illness.Severe disease is characterized by rapid onset of high fever, headache, neck stiffness, disorientation, coma, seizures, spastic paralysis, and ultimately death.[1,2].The case fatality rate for the disease can be as high as 30% among those with disease symptoms.In addition, 20-30% of those who survive suffer permanent neuropsychiatric sequelae [3].
Leptospirosis is a bacterial disease that affects humans and animals.Leptospirosis is considered to be the most widespread zoonotic disease in the world [5].It is caused by bacteria of the genus Leptospira [6].The bacteria are spread through the urine of infected animals, which can get into water or soil and can survive there for weeks to months.Humans can become infected through contact with urine (or other body fluids, except saliva) from infected animals; also, through contact with water, soil, or food contaminated with the urine of infected animals [7].In humans, leptospirosis may occur in two phases, where about 10% of infected people move to the 2nd phase [5].The first phase causes a wide range of symptoms, some of which may be mistaken for other diseases: fever, chills, headache, muscle aches, vomiting, or diarrhea.The second phase, if it occurs, is more severe; the person may have kidney damage, liver failure, meningitis, respiratory distress, and even death [6,8].It is estimated that more than 1 million cases occur worldwide annually, including almost 60,000 deaths [8].
In the last couple of decades, many ecological and field studies, along with a few mathematical ones, have been conducted to understand the dynamics of JEV transmission and to find some control measures to reduce JEV prevalence in humans.In 2001, a study used a probabilistic model of pathogen transmission to investigate various control measures for JEV transmission in humans [9].Outcomes of the study show that a combination of control measures of similar effect (strategies to reduce vector population, or strategies to reduce vector-human interactions) is more effective compared to the combination of control measures of different effects (strategy to reduce vector population and strategy to reduce human-vector interaction).In 2014, Khan et al. studied the dynamics of JEV transmission in a pig population in northwest Bangladesh [10].This study developed an SEIR model to understand transmission dynamics in pigs, and to estimate the potential impact of pig vaccination.Their results found that the prevalence of JE in pig populations can be reduced by up to 89% when 75% of susceptible pigs are vaccinated each year.The next year, in 2015, Lord et al. performed a study to rethink JEV transmission among hosts and vectors [11].They suggested using a mathematical model parameterized with data to quantify the relative roles of potential species in JEV transmission.
Very little research has been done, to the best of our knowledge, to understand the dynamics of leptospirosis in cattle and humans.In 2017, Chadsuthi et al. investigated the leptospirosis prevalence in livestock and humans in Thailand for 2010-2015 [12].They tested humans, buffaloes, cattle, pigs, and analyzed collected data.Their analysis found livestock more susceptible to leptospirosis infection compared to humans.Later, in 2018, another study was done to understand the spread of leptospirosis in lambs in New Zealand.Here, researchers used a simple SI model to predict conditions under which the disease would persist in the lamb population [13].Analysis of this study suggested that increasing the leptospira death rate in farms can reduce infection in livestock, and eventually in humans.
There are many countries in Asia where both JE and leptospirosis are prevalent.Taking this reality into account, this study is developed to investigate the prevalence of both diseases in humans.Cattle contribute to leptospirosis infections in humans, acting as a source of leptospira.However, the presence of cattle in a domestic or peridomestic setting can be considered as an additional host for JE vectors, where humans act as the primary host.Hence, the presence of cattle increases host richness in the setting, which might eventually reduce the combined burden of these two diseases for humans.Earlier studies showed that an additional host (host richness) in a setting helps to reduce human infections of vector-borne diseases if certain conditions are maintained [16,17,18,20].Johnson and Thieltges in 2010 identified that host diversity may help to reduce human infections depending on the relative abundance of additional host(s) relative to the focal host [16].In 2013, Miller and Huppert proved that species diversity in host populations can amplify or can dilute disease prevalence depending on vectors' preference of host [17].Recently in 2020, Zahid and Kribs showed that an additional host (dog, reservoir host) in a community can help in reducing the prevalence of visceral leishmaniasis in humans depending on the additional host's irritability to vector bites [18].In another study, the same authors established that the presence of an additional host (chicken, incompetent host) in a domestic setting might help to reduce the prevalence of Chagas disease in humans if the distance between the two host populations (humans and chickens) is maintained at a certain range [20].
The presence of cattle helps to reduce human-vector interactions by attracting mosquitoes towards them from humans.Cattle have no contribution to the transmission of JE [21] whereas they act as a source of leptospirosis infections in humans.Thus, in terms of infections, the presence of cattle has two opposite influences on human health -it helps to reduce JE prevalence in humans and contributes to leptospirosis risk in humans.The goal of this study is to understand the dynamics of both diseases in each population, and eventually to understand if the presence of cattle in a setting, where JE and leptospirosis both are prevalent, is helpful to reduce the combined burden of JE and leptospirosis in humans.To answer this question we compare two different peridomestic settings: a setting involving cattle with humans, pigs, and mosquitoes, and the other setting involving humans, pigs, and mosquitoes without cattle.To understand disease dynamics, and to compare the proposed two settings, we develop SIR models for cattle, pig and human populations, and an SI model for the mosquito population.This study estimates the total number of DALY s (Disability-Adjusted Life Year) generated from JEV and (in the presence of cattle) leptospirosis infections for both settings and compare them to identify the better scenario.To the best of our knowledge, this is the first work ever to consider JE and leptospirosis infections together.More broadly, this extends the well-studied question in disease ecology of the impact of an additional host to a multi-pathogen context.

Model development
The goal of this work is to understand the dynamics of Japanese encephalitis virus (JEV) and leptospirosis infections together and to estimate the impact of cattle presence on human infections of these diseases.Therefore, we select a location for this study where both diseases are prevalent.We choose Malkangiri, the southern district of Odisha (Orissa) State in India.Malkangiri has a history of JE outbreaks since 2009 where the cases are documented as Acute Encephalitis Syndrome (AES) [22].This AES can include infections other than JE.Studies showed that AES cases are not always cases of JE infections, they can be cases of the severe form (the 2nd phase) of leptospirosis infections [23,24,25].This fact, along with the high prevalence of leptospirosis in cattle in neighbor districts [26,27], illustrates the presence of human leptospirosis cases in our study location.
We use an SIR model to understand the dynamics of JEV among pig, mosquito, and human populations; also, to understand the dynamics of leptospirosis between cattle and human populations.JEV mainly exists in a transmission cycle between mosquitoes and pigs, and humans get the infection of JE from bites of infected mosquitoes (I m ) [1, 2].However, infected humans do not infect feeding mosquitoes due to the lack of sufficient viremia.So, humans are dead-end hosts for JEV [1,2,3,36].Also, cattle have no role in the maintenance of JEV in nature [21].Thus, humans and cattle have no contribution to JEV transmission.We consider total number of pig bites as λ r = bmNm krNr+kcNc+N h k r N r and total number of human bites as λ h = bmNm krNr+kcNc+N h N h , where N m , where N r , N c , and N h represent sizes of mosquito, pig, cattle and human populations respectively, b m N m represents total mosquito bites.Here, k r and k c are defined based on mosquitoes' feeding preference among hosts (see Table 2).In our model development, we consider the vertical transmission of JEV in pigs, and mosquitoes.All these ideas are illustrated in Figures 1, 2, and 3.
Leptospirosis spreads through the urine of infected cattle or pigs and is transmitted through the contact of contaminated (by the urine of infected animal) surfaces or water [7].Here, we ignore leptospirosis infections through pigs, because people usually keep pigs in a separate place from houses.Another important reason is its population size (N r ), which is very small compared to the cattle population size (N c ).Therefore, in our model cattle are the only source of Leptospira, which causes the disease.Hence, cattle and humans are getting infections just due to the presence of cattle (Figure 2 and Figure 3).Here, we consider that Leptospira can be transmitted in cattle through vertical transmission also.
Susceptible humans get leptospirosis infections and move to SL1 h from SS h .Most leptospirosis infected humans recover from the first phase, while the remaining move to the second phase of the infection (SL2 h ).Finally, patients from the second phase move to SR h upon their recovery.The SR h compartment represents people who already recovered from leptospirosis infection; however, they are still susceptible to JEV.A portion of susceptible humans get JE infections by infected mosquito bites before leptospirosis infection and move to ES h .Some of these infected humans get leptospirosis infection before their recovery from JE and move to EL1 H .The remaining population of ES h move directly to the JE recovered class (RS h ).Similar to SL1 h , most people from EL1 h move to EL2 h while others move to ER h .Recovered people from JEV infections (RS h ) get leptospirosis infections and move to RL1 h .Then, similar to SL1 h and El1 h , most of them move to RR h after recovery from leptospirosis, while others enter the second phase of infections (RL2 h ).Finally, people recover from RL2 h and move to RR h , which contains people who recovered from infections of both diseases.
In our model development, we consider constant populations for cattle, pigs, and mosquitoes.For disease transmission, we choose standard incidence for JEV, because infected mosquitoes (I m ) are free to bite any individual among our host populations.On the other hand, we use mass-action incidence for Leptospira transmission, because the new infections depend on the availability of infected cattle (L c ).Also, we assume our cattle, pig and mosquito populations to be constant.

Parameter estimation
The World Bank recorded the life expectancy of people at birth in India in 2017 as 69.165 years [28].Taking the reciprocal of this value we get the natural death rate for humans as µ h = 3.96×10 −5 /day.As per the 2014 census by the local health department, Malkangiri district had a population of 641,385 in 2014 with 109,483 households [29].However, according to the World Bank data rural population in India in 2014 and 2018 was 876,035,725 and 892,321,651 respectively [30].We use the ratio of these two populations to estimate the population of Malkangiri district in 2018, which gives us the population as 653,309 (N h ) with 111,518 households.In our literature review, we did not find any value for the human recruitment rate.Hence, we use one of our equations (9th in the system) in a disease-free state which gives the equation Λ h = µ h N h .Then we use this equation, and our already estimated values of µ h and N h , to get Λ h = 25.87 human/day.We did not find any documented data for the size of cattle and pig populations (N c and N r respectively) in Malkangiri district.However, we found one documented data source which men- .Now we need to know the feeding pattern of JEV vectors (value of k c and k r ) to estimate λ r and λ h .We found a study that estimated the feeding preference for JE vectors in the southern part of India as 46.4% on cattle, 4.8% on pigs, and 1.5% on humans [31].Using these values we have k c = 46.41.5 human/cattle and k r = 4.8 1.5 human/pig.
To estimate λ r and λ h , we also need to have the size of JEV carrier mosquito populations (N m ), along with their biting rate (b m ).A field study in rural villages of Western Yunnan Province of China was done in 2013 to estimate the abundance of mosquitoes in an Asian rural setting [32].Researchers collected mosquitoes from two households in each of four studied villages in each month for a 12 month-study.The total collection of mosquitoes was 85,307.It is equivalent to 85,307 (2×4)×12 =888 mosquitoes/household. So, for our study we have a total of (111,518×888)=99,027,984 mosquitoes.Another study carried out in some villages of Malkangiri district showed that 69.1% of available mosquitoes are vectors for JEV [22].Hence, the total number of JEV vectors we have is N m =99,027,984×69.1%=68,428,337.We found documented data for the biting rate of JEV vectors which give b m =0.25 bite/mosquito-day (mean of 0.2 and 0.3) [33].Anyway, not all mosquito bites belong to these three hosts.Also, blood-meal analysis showed that some mosquitoes bite more than one host [31].Hence, we estimated that about 70% of mosquito bites are distributed among cattle, pig and human population.Thus, λ r = 0.25×68,428,337×70% 4.8 1.5 ×37,616+ 46.4   1.5 ×302,911+653,309 × ( 4.8 1.5 × 37, 616)=142,102 bite (total)/day and λ h = 0.25×68,428,337×70% 4.8 1.5 ×37,616+ 46.4   1.5 ×302,911+653,309 × 653, 309=771,250 bite (total)/day.In 2000, the Indian Council of Medical Research carried out a study which found the transmission and infection rate for JEV vectors [34].The study showed that vectors have 82% infection probabi-  [36].Taking the reciprocal of this value we get γ r = 1 4 /day.For JE vectors' mortality rate, we found a study that estimated their life expectancy as 59.8 days [37], and this gives us γ m = 1 59.8 /day.Humans' mean recovery period during the 2012 JE outbreak in Odisha was 7 days [35], which gives us γ e = 1 7 /day.For the estimation of µ r , we found that domestic pigs have an average lifespan of 6-10 years, but, that can be shorter due to certain problems [38].Thus, instead of the median (8 years) of the range, we choose 7 years as the life expectancy for pig population which eventually estimated µ r = 1 7×365 /day = 3.91 × 10 −4 /day.
A review study from 2014 mentioned the mean durations for the acute phase (1st phase, in our model) and immune phase (2nd phase, in our model) as 10 days and 14 days respectively [39].Taking reciprocals of these values we get γ l1 = 1 10 /day and γ l2 = 1 14 /day.Another study in coastal south India estimated that 51 out of 202 (25.24%) leptospirosis patients move to the 2nd phase [40].Hence, we use this result and the ratio r l γ l1 +r l +µ h = 51  202 to estimate the transfer rate from phase 1 to phase 2, and found r l = 0.033/day.We found a couple of documented data for cattle recovery period of leptospirosis which are inconsistent in values [41,42,43].Hence, based on those available values we choose that cattle can shed leptospira in their urine for 60 weeks, which gives γ lc = 1 60×7 /day = 2.38 × 10 −3 /day.In 2016, a different disease in cattle in Odisha state was studied which had the highest cattle age group of 6.5-7.5 years; however, they did not mention the highest or average age [44].So, we take 8 years as the lifespan for cattle in Odisha, which gives us µ c = 1 8×365 /day = 3.42 × 10 −4 /day.
We did not find enough evidence for vertical transmission of leptospirosis in cattle (p c ) and of JEV in pigs (p r ), and therefore we assume p c = 0 and p r = 0.However, we found 4% effective vertical transmission for JEV vectors, which gave p m = 0.04 [45].In 2016, in Malkangiri district 37 died out of 175 JE patients [29].Using this ratio and the relation αe αe+γe+µ h = 37  175 , we estimated α e = 0.0383/day.The case fatality rate for leptospirosis infection is 7.69% [50].Also, we already know that about 25.24% leptospirosis cases move to the 2nd phase [40].Hence, the case fatality rate among patients of 2nd phase is 7.69  25.24 × 100, which implies α l = 3.13 × 10 −2 /day.
During our literature review, we did not find any value for infection rates of leptospirosis in cattle and in humans (β lc and β lh respectively).Hence, we follow the method used in [20,19] to estimate β lc and β lh which gives where y c and y h are prevalence of leptospirosis in cattle and in humans respectively.We did not find any documented data for y c in Malkangiri.However, we found one study from 2013 which estimated leptospirosis prevalence as 20.7% in cattle in a neighbour district of Malkangiri [27].We choose y c = 20.7% 2 = 11.35%since the sample collections of this study were not random, rather were mostly from the villages with a history of abortions and other disorders.Also, we did not find documented data for leptospirosis prevalence in humans in Malkangiri.However, we found human and cattle prevalence of leptospirosis in another prevalent location of leptospirosis, northeast Thailand.We estimated leptospirosis prevalence in cattle for this area to be 34+2+58 130+183+238 = 94 551 = 17% using the data from a study published in 2016 [47].A different study estimated 12.5 annual cases of human leptospirosis per 100,000 people [48].We used the ratio of estimated leptospirosis prevalence in cattle for our study area (11.35%) and of northeast Thailand (17%), and the annual cases of human leptospirosis in northeast Thailand (12.5 per 100,000) to calculate annual cases of human leptospirosis for our study area, and found 12.5×11.3517 = 8.35 cases per 100,000 people, which is equivalent to 0.00835%.Also, we have a weighted recovery period as (10+14)×25.24%+10×74.76%100% = 13.53 days, which gives us 365  13.53 = 27 as the number of generations in a year for infected people.Finally, we got y h = 0.00835% 27 = 3.09 × 10 −6 .Using these prevalence values and our expressions for leptospirosis infection rates, we got β lc = 9.66 × 10 −8 /day-cattle and β lh = 9.11 × 10 −12 /day-cattle.

Analysis
We begin with the general model where we consider the presence of cattle which includes both JE and leptospirosis diseases.Later, we consider the scenario without cattle; this special scenario has JEV infection only.Finally, we estimate total disease burdens for these two scenarios and compare them to understand the impact of cattle on disease burden in JE prevalent areas.The goal of this study is actually to understand how the presence of domestic animals, here cattle, helps to reduce disease burden.
First, we analyse the cattle system, which is decoupled from the other two subsystems -the pig-mosquito system, and the human system.We set all the equations of this system equal to zero and solve them for state variables to find equilibria values.After performing some basic algebra we get the disease free equilibrium (DF E) as S c = N c , L c = 0, R c = 0, and the endemic equilibrium (EE) as Then we use the next generation method [49] to obtain the basic reproduction number (BRN ) which gives R 0l = β lc Nc+µcpc γ lc +µc .We find that the EE exists under the condition R 0l > 1.Here the EE is globally asymptotic stable if and only if R 0l > 1, otherwise the DF E is globally stable (see appendix 1 for details).
Before we begin our discussion on the pig-mosquito subsystem, here we consider a simplifying assumption that the two diseases do not significantly affect the magnitude of the human population, so that we can replace N h (t) by N h (0) in the denominator of λ r and λ h , in order to simplify the analysis.This assumption is supported by both our literature review and the later numerical analysis.As per the study [29], JEV infection causes 37  641,385 × 100, 000 = 5.77 annual deaths per 100,000 people.Also, leptospirosis infection causes 12.5 × 7.69 = 0.96 annual death per 100,000 people [48,50].These two estimations give us a total of 44 (= 653,309 100,000 × (5.77 + 0.96)) annual deaths due to both diseases, which is 0.0067% of the total human population.Under this assumption, mosquito's biting rates to pigs, and humans ( λr Nr and λ h N h respectively) become constants, and the vector-reservoir (pig-mosquito) subsystem is decoupled from the human subsystem.The pig-mosquito subsystem is already decoupled from the cattle subsystem.Now, following the same approaches, like the cattle system, for the pig-mosquito system, we get the DF E as S r = N r , I r = 0, R r = 0, S m = N m , I m = 0, and the EE as ( We get the BRN for the pig-mosquito system as Here the EE exists and is locally asymptotically stable (LAS) if and only if R 0e > 1, otherwise the DF E is LAS.Next, we develop a strong Lyapunov function to show that the DF E is globally stable when R 0e < 1 (see appendix 2 for details).Then, we analyse the stability of the EE numerically which indicates that the EE here is globally stable if and only if R 0e > 1.
The cattle subsystem and the pig-mosquito subsystem are decoupled from the human subsystem respectively by the model and the assumption that the magnitude of the human population is insignificantly affected by the effect of JE and leptospirosis.As these two subsystems go to equilibrium, the limiting system for system 1 is given by the human subsystem with L c and I m replaced by their equilibrium values (please see system (3)).Here N m is constant which makes the term I * m Nm to be constant and λ h0 N h (= bmNm krNr+kcNc+N h0 =constant) is also constant.Thus, the resulting system (system (3)) is linear and has only one equilibrium (see equation (9) in Appendix 3) where the equilibrium values are functions of L * c and I * m .Hence, the behavior of the limiting system follows the behavior of our cattle and pig-mosquito subsystems.Now, a theorem by Horst R. Thieme [51] leads us to the fact that the behavior of the entire system is asymptotic to the behavior of the limiting system, which is linear.Eventually, the entire system is governed by the behavior of the cattle and the pig-mosquito subsystems.
Next, we extend our analysis to the entire system, which includes the cattle, the pig-mosquito, and the human subsystems.We obtain the invasion reproduction number (IRN ) for both diseases following the approach by Mitchell & Kribs [52].Each of these IRN s is found to be the same as the corresponding BRN s.Thus, both of the IRN s are independent of the presence of the other disease.This result is easily comprehensible since neither of these diseases' transmissions is affected by the other.Our preceding analyses found two different values for both L * c and I * m ; hence, their all possible combinations produces four different equilibria for the entire system.These four equilibria are -DF E, endemic in leptospirosis and free of JE, endemic in JE and free of leptospirosis, and endemic in both diseases respectively.So, for the entire system, we have four different scenarios -(i) R 0l < 1 and R 0e < 1 (ii) R 0l > 1 and R 0e < 1, (iii) R 0l < 1 and R 0e > 1, and (iv) R 0l > 1 and R 0e > 1. where Now we consider the scenario without cattle which eliminates the possibility of leptospirosis infections in humans.Here the BRN for JE (R 0e ) is different compared to the BRN (R 0e ) from from the general model since the λ r in R 0e is a function of cattle population size.The absence of cattle reduces the human system to a three-dimensional system, which gives the EE for the human system as where λ h0 = bmNm krNr+N h0 N h (here N h0 = N h (0)).This special case has DF E, and the only endemic equilibrium (EE) when R 0e > 1.This EE is completely governed by I * m (see equation ( 2) for I * m ).
To better understand the impact of cattle, some numerical analyses are done based on our parameter estimations.Using our estimated parameter values, we find R 0l = 10.35 and R 0e = 1.008 respectively in the presence of cattle.However, our estimations in the absence of cattle found R 0e = 12.97.Based on our estimated parameter values, in the presence of cattle we estimate 72 and 228 annual cases for leptospirosis and JE respectively, whereas the annual deaths are 6 and 48 respectively.In the absence of cattle, we estimate 9,407 and 1,988 number of annual cases and annual deaths due to JEV infections.Next, we perform numerical analysis to estimate total cases and the BRN of JE and leptospirosis as the average number of cattle per household varies and plot the related graphs (Figure 4).The Figure 4(a) shows that the JEV infection decreases as the number of cattle increases; however, the human leptospirosis incidence increases with the cattle; the Figure 4(b) illustrates that BRN for JE exponentially decreases, and for leptospirosis linearly increases as the number of average cattle increases.So, the presence of cattle has positive and negative impacts on human health.
To quantify the impact of cattle on human health, we calculate the total disease burden for our two settings -with cattle and without cattle.In this calculation we follow an approach similar to the approach of the study [53] where they followed WHO guidelines [54].To calculate the total disease burden they estimated total DALY (Disability-Adjusted Life Year) where one DALY means the loss of one year of healthy life.The DALY is composed of Y LD (Years Lost due to Disability) resulting from infections and Y LL (Years of Life Lost ) caused by disease-induced premature deaths which are calculated using formulas Y LD = I × DW × L 1 and Y LL = D × L 2 , where I and D represent the total number of infections (cases) and total number of deaths respectively.Here L 1 = average duration (in years) of illness which is the reciprocal of the recovery rate for survival cases and of the death rate for non-survival cases, and L 2 = standard life expectancy at age of death (in years) = average life expectancy of the community -average age at premature death and the term DW is the disability weight for diseases which ranges from 0 (perfect health) to 1(death).It can be thought of as the proportional reduction in perfect health due to any adverse health condition caused by infections.So, the total burden of disease can be represented as the sum of Y LD and Y LL, which gives In this work, we study two different diseases, JE and leptospirosis.Hence, we will have different values of I, D, L 1 , L 2 and of DW for JE, and leptospirosis.Here the quantities L 1 s are the reciprocals of the γ l1 , γ l2 for the first and second phases of leptospirosis respectively, and the reciprocal of the γ e for JE.To estimate L 2 s, we need to know the average ages of infections which are 40.48 years [40] and 5 years [35] for leptospirosis and JEV infections respectively.The average age of infections and the average age of deaths are the same for leptospirosis; however, they are different for JE (5 years and 3 years respectively [35]).The values of DW differs between leptospirosis and JE; it also differs between phases of leptospirosis.Furthermore, the DW has different values for survival and non-survival.All these different values of DW are taken from the study [55] based on the severity (mild or severe) of illness and on the onward complications upon recovery.Therefore, using the above formula for DALY we calculate total DALY s separately for JE and leptospirosis infections and add them to have total DALY cause from diseases combined.
Based on our parameter estimations, the above formula and values together calculate 4,396 DALY s and 157,642 DALY s for the setting with cattle and without cattle respectively.These results clearly show that the presence of cattle has a huge positive impact on the disease burdens in JE-prevalent areas.Next, we calculate total DALY s varying the average number of cattle per household to understand how the variation in the average number of cattle per household changes annual disease burden.The results of these calculations are portrayed in Figure 5, which evidently illustrates that the presence of cattle is unconditionally helpful in terms of the disease burden.However, there is an optimum value for the average number of cattle per household to ensure the minimum disease burden which is 2.75.Maintaining cattle beyond this threshold value the annual disease burden increases slowly which is generated from the infections of leptospirosis.On the other hand, keeping fewer cattle per household than the optimum number causes a sharp rise in the annual disease burden due to a significant increase in JE infections.

Discussion
To the best of our knowledge, this is the first work to understand the joint dynamics of JE and leptospirosis.Our qualitative analysis exhibits the classic threshold behavior of the system.The quantitative analysis shows that the presence of cattle reduces total disease burden even though it introduces leptospirosis.It happens because the number of JE cases is reduced by the cattle presence at a higher rate than the rate at which humans get leptospirosis infections due to the presence of cattle.Also, the case fatality rate (CFR) of JE infections is higher than the CFR of leptospirosis infections.The disability weight is higher for JE infections than for the leptospirosis infections, and the average age of JEV infections is very low compared to the average age of leptospirosis infections.
Even though the presence of cattle is always helpful, there is an optimum average number of cattle per household which minimizes the annual disease burden.
Our analysis found that the presence of cattle helps to reduce the annual disease burden in JE-prevalent areas because of the vector's (mosquitoes') feeding preference.Mosquitoes prefer to bite cattle at a much higher proportion compared to the proportion to bite humans or pigs.Thus, cattle reduce mosquito-human interactions and eventually reduce human cases of JEV infections.This result is consistent with the result of an ecological study by Miller and Huppert [17] where they identified that vectors' feeding preference decides whether or not host richness is helpful to dilute disease prevalence.Also, our results found that the JE prevalence in humans decreases as the cattle population increases.This finding is consistent with the finding of another ecological study on vector-borne diseases by Johnson and Thieltges [16] where they found that the strength of the dilution effects depends on the relative abundance of the dilution host to the focal host.
This study didn't consider seasonality in its model development.Considering seasonality will provide a more accurate picture of short-term dynamics.We didn't include seasonality because the concern of our study is primarily with the overall disease burden, not with the detailed disease dynamics.In our analysis, we assumed constant biting rates to pigs and humans which should not have any significant impact on our results since the mosquito population is not limited by the available host populations, rather by available breeding sites.We also assumed a homogeneous distribution among households and those cattle close enough to households to affect the infection risks for JEV and leptospirosis.Large cattle farms farther away from homes than pigs are, are not likely to affect those risks (except to employees).However, studies have shown that heterogeneity does affect disease risk [56,57].In our case, heterogeneity is the heterogeneous distribution of household cattle per household.

Appendix 1. Stability of the endemic equilibrium of the cattle system
Proposition 1: The endemic equilibrium for the cattle subsystem is locally stable if and only if R 0l > 1. Proof: We find the Jacobian matrix of the cattle system at the EE which gives the characteristic equation Here, λ = −µ c is negative.Next, we use Routh-Hurwitz Criteria to check if the real part of remaining eigenvalues are negative.As per Routh-Hurwitz Criteria, the coefficients of λ and the constant term need to be positive to ensure the negativity of real parts of eigenvalues.Here it is evident that the coefficient of λ is always positive and the constant term is positive when R 0l > 1.
Hence, the EE for the cattle subsystem is locally stable if and only if R 0l > 1.
Proposition 2: The endemic equilibrium for the cattle subsystem is globally stable when R 0l > 1, otherwise the DF E is globally stable.Proof: The first two equations of the cattle subsystem are It is evident from the system (4) that S c > 0 if S c = 0 and L c = 0 if L c = 0, which means the system (4) is well posed.Now, the system (4) gives So, no limit cycle exists.
Hence, by Poincaré-Bendixson Theorem all solutions approach the EE when R 0l > 1 which means the EE is globally stable if and only if R 0l > 1, otherwise the DF E is globally stable.

Appendix 2. Stability of the endemic equilibrium of the pig-mosquito system
Proposition 3: The endemic equilibrium for the pig-mosquito subsystem is locally stable if and only if R 0e > 1. Proof: We can reduce our system to a system of three equations since pig and mosquito populations are constant.We Replace S r = N r − (I r + R r ) and S m = N m − I m in the 3rd and 5th equation of the system (1) and get the reduced pig-mosquito system as Next, as per Routh-Hurwitz Criteria, we verify A 1 , A 2 > 0 and A 1 A 2 > A 3 to check if the real parts of all eigenvalues are negative.Here, A 1 is clearly positive.Also, our analysis found A 2 > 0 and A 1 A 2 > A 3 .
Hence, the endemic equilibrium for the pig-mosquito subsystem is locally stable if and only if R 0e > 1.
Proposition 4: The DF E for the pig-mosquito subsystem is globally stable when R 0e < 1. Proof: We define our Lyapunov function as Here we use a alternative threshold quantity R to understand if the function V satisfies all the conditions of a strong Lyapunov function which is defined as where R > 1 if and only if R 0e > 1.Now, the alternative threshold value R, equation (7) and equation (8) easily verify that V (0, 0, 0) = 0 and V (I r , R r , I m ) > 0 for (I r , R r , I m ) = (0, 0, 0) when R < 1.Also, V (0, 0, 0) = 0 and V (I r , R r , I m ) < 0 for (I r , R r , I m ) = (0, 0, 0) when R < 1.

Figure 1 :
Figure 1: Flow diagram for pig and mosquito population

Figure 4 :
Figure 4: The effect of cattle on infections and BRNs

Figure 5 :
Figure 5: Relation between avg cattle/household and annual disease burden (in DALY s) or L c ≥ N c which shows all solutions are bounded.Now we define β(S c , L c ) = 1 ScLc to use Dulac's Criterion to check if the EE approaches a limit cycle.Here we assume F (S c , L c ) = S c and G(S c , L c ) = L c which give

Table 1 :
Model variables with definition Variable Definition S c , S p , S m Susceptible cattle, pigs, mosquitoes L c , I r , I m Infected cattle, pigs, mosquitoes R c , R r Recovered cattle, pigs SS h Humans susceptible to both JE and leptospirosis ES h Humans infected with JE, but susceptible to leptospirosis RS h Humans recovered from JE and susceptible to leptospirosis SL1 h Humans susceptible to JE, but infected with 1st phase of leptospirosis EL1 h Humans infected with JE and 1st phase of leptospirosis RL1 h Humans recovered from JE and infected with 1st phase of leptospirosis SL2 h Humans susceptible to JE, but infected with 2nd phase of leptospirosis EL2 h Humans infected with JE and 2nd phase of leptospirosis RL2 h Humans recovered from JE and infected with 2nd phase of leptospirosis SR h Humans susceptible to JE, but recovered from leptospirosis ER h Humans infected with JE, but susceptible to leptospirosis RR h Humans recovered from both JE and leptospirosis tioned that the numbers of humans, cattle and pigs in Korkunda block of Malkangiri district in 2014 were 143,867, 80,583 (75,772+4,811), and 10,007 respectively in 29,667 households [22].We use these numbers to estimate cattle and pigs per household, and then use our estimated number of households (111,518) to calculate numbers of total cattle and total pigs in the district in 2018, which gives N c = 80,583 29,667 × 111, 518 = 302, 911 and N r = 10,007 29,667 × 111, 518 = 37, 616