# Modelling Japanese encephalitis virus transmission dynamics and human exposure in a Cambodian rural multi-host system

Héléna Ladreyt, Véronique Chevalier, Benoit Durand

► To cite this version:

Héléna Ladreyt, Véronique Chevalier, Benoit Durand. Modelling Japanese encephalitis virus transmission dynamics and human exposure in a Cambodian rural multi-host system. PLoS Neglected Tropical Diseases, 2022, 16 (7), pp.e0010572. 10.1371/JOURNAL.PNTD.0010572. hal-03745907

HAL Id: hal-03745907

https://hal.inrae.fr/hal-03745907

Submitted on 4 Aug 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Modelling Japanese encephalitis virus transmission dynamics and human exposure in a Cambodian rural multi-host system

Héléna Ladreyt¹, Véronique Chevalier¹,²,³,⁴, Benoit Durand⁵*

1 International Center of Research in Agriculture for Development (CIRAD), UMR ASTRE, Montpellier, France, 2 CIRAD, UMR ASTRE, Antananarivo, Madagascar, 3 Epidemiology and Public Health Unit, Institut Pasteur du Cambodge, Phnom Penh, Cambodia, 4 Epidemiology and Clinical Research Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar, 5 Epidemiology Unit, Laboratory for Animal Health, French Agency for Food, Environmental and Occupational Health and Safety (ANSES), University Paris-Est, Maisons-Alfort, France

* benoit.durand@anses.fr

## Abstract

Japanese encephalitis (JE) is a vector-borne zoonosis and the leading cause of human viral encephalitis in Asia. Its transmission cycle is usually described as involving wild birds as reservoirs and pigs as amplifying hosts. JE is endemic in Cambodia, where it circulates in areas with low pig densities (<70 pigs per km²), and could be maintained in a multi-host system composed of pigs, but also poultry as competent hosts, and dogs, cattle and humans as non-competent hosts. We used a mathematical model representing Japanese encephalitis virus (JEV) transmission in a traditional Cambodian village that we calibrated with field data collected in 3 districts of Kandal province, Cambodia. First, R₀ calculations allowed us to assess the capacity of the epidemiological system to be invaded by JEV and sustain virus transmission in villages in the 3 districts, and we predicted human exposure at the epidemiological equilibrium, based on simulations. Changes in spatial density of livestock, in agricultural practices, and epizootics (e.g., African swine fever), can profoundly alter the composition of host communities, which could affect JEV transmission and its impact on human health. In a second step, we then used the model to analyse how host community composition affected R₀ and the predicted human exposure. Lastly, we evaluated the potential use of dog JE seroprevalence as an indicator of human exposure to JEV. In the modeled villages, the calculated R₀ ranged from 1.07 to 1.38. Once the equilibrium reached, predicted annual probability of human exposure ranged from 9% to 47%, and predicted average age at infection was low, between 2 and 11 years old, highlighting the risk of severe forms of JEV infection and the need to intensify child immunization. According to the model, increasing the proportion of competent hosts induced a decrease in age at infection. The simulations also showed that JEV could invade a multi-host system with no pigs, reinforcing the assumption of poultry acting as reservoirs. Finally, the annual human exposure probability appeared linearly correlated with dog seroprevalence, suggesting that in our specific study area, dog seroprevalence would be a good proxy for human exposure.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Ladreyt H, Chevalier V, Durand B (2022) Modelling Japanese encephalitis virus transmission dynamics and human exposure in a Cambodian rural multi-host system. PLoS Negl Trop Dis 16(7): e0010572. https://doi.org/10.1371/journal.pntd.0010572

**Editor:** Elvina Viennet, Australian Red Cross Lifefood, AUSTRALIA

**Received:** May 12, 2021

**Accepted:** June 8, 2022

**Published:** July 11, 2022

**Copyright:** © 2022 Ladreyt et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** This study was funded by the French Ministry of Agriculture (HL), the International Centre for Research in Agriculture and Development (Cirad), and the French Agency for Food, Environmental and Occupational Health and Safety (Anses). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
**Competing interests:** The authors have declared that no competing interests exist.

## Author summary

Japanese encephalitis virus (JEV) is endemic in Cambodia and remains the most common cause of acute viral encephalitis, particularly in children and adolescents. The traditionally described cycle of JEV, involving wild birds as reservoirs, pigs as amplifying hosts and *Culex* mosquitoes as vectors is questioned, with increasing evidence of a more complex multi-host system involved in areas where densities of pigs are low. In Cambodia, the infection could be maintained in a multi-host system consisting of pigs and poultry as competent hosts, and dogs, cattle and humans as non-competent hosts. We defined a compartmental dynamic model of JEV transmission in a multi-host system representing a rural Cambodian village, to predict human exposure to JEV in the studied area, and to analyse how host community composition may affect human exposure and R₀ value. Our theoretical approach showed that variations of the composition of the multi-host system may have an impact on human exposure to JEV, and thus on the disease burden in humans, especially in young children. Besides children vaccination in JEV endemic areas, a proper evaluation of the impact on human health is needed to target prevention actions and reduce JEV burden in Cambodia.

## Introduction

Japanese encephalitis virus (JEV) is the leading cause of human acute viral encephalitis in Asia. Even if the global burden of JEV on human health is difficult to evaluate, this vaccine-preventable vector-borne zoonosis was estimated to have caused around 100,000 cases and more than 25,000 deaths in 2015 [1], with an observed case fatality rate up to 30%. Severe neurological sequelae may persist in 30% to more than 50% of survivors [2-6]. Approximately three quarters of cases concern children and JE remains a substantial public health issue even in areas where human vaccination programs are implemented [7]. Furthermore, JEV exposure figures are likely to be underestimated due to under-detection and under-diagnosis of acute encephalitis, particularly in developing countries, as well as cross-reactivity of serological tests with other flaviviruses, especially dengue virus [8,9].

It is commonly accepted that JEV is transmitted from Ardeid birds (wild reservoir hosts) or domestic pigs (amplifying hosts) to humans through the bites of *Culex*, and probably some *Aedes* mosquitoes [10,11]. However, JEV circulates in areas with low densities of domestic pigs or Ardeid birds, such as in Singapore, where its circulation was detected in sentinel poultry years after pig farming was phased out [12]. Thus, the epidemiology of JE may differ from one region to another [13,14]. Domestic chickens and ducks have been shown to be exposed to JEV in several regions in Asia [8,15-18]. A serosurvey conducted in Kandal province, Cambodia, showed that JEV force of infection, i.e. the instantaneous probability to become infected, exerted on ducks was comparable to that exerted on pigs [16]. Chicks and ducklings also develop significant viremia after being bitten by infected mosquitoes [19-22] and may be able to infect susceptible mosquitoes feeding on them [23,24], suggesting that young poultry may be competent hosts, i.e. able to transmit JEV to vectors [25]. However, no study has yet fully investigated the potential of poultry as competent hosts for JEV under natural conditions.

*Culex spp.* mosquitoes, the main vectors of JEV, are opportunistic and can feed on various host species, depending on local abundance of hosts and vector species-specific intrinsic factors [26-28]. These hosts can be competent, such as pigs, waterfowl and poultry [8,12,29,30], or not, such as cattle, dogs, humans and horses, which may “dilute” mosquito bites [31-34].
The basic reproduction number ($R_0$) is the expected number of secondary cases generated by a primary case in an entirely susceptible population [35]. $R_0$ is the indicator commonly used to measure whether or not a pathogen can invade a population. $R_0$ greater than 1 is a necessary but not sufficient condition for self-sustained transmission, which also requires a renewal of the susceptible and competent host pool. Variations in the composition of the multi-host system may thus affect $R_0$ value.

Japanese encephalitis is endemic in Cambodia, and remains the most common cause of acute encephalitis, particularly in children and adolescents [33,36]. Cambodia is a predominantly rural country, although the density of domestic pigs is low (around 20 pigs per km² on average) compared to other countries where JEV is circulating, such as Japan, China or Vietnam where densities can reach 700 pigs per km² [13,37–39]. In Kandal province, a rural area of Cambodia surrounding the capital Phnom Penh, the majority of the livestock is raised in backyards, with close proximity between pigs, chickens, ducks, cattle and humans. In this area, a recent survey showed a yearlong presence of JEV vectors feeding on pigs, chickens, humans, dogs and cattle [26,40]. Although JEV host densities are low (353 chickens and 66 pigs per km² in Kandal province [37,41]), this multi-host system could allow for year-round circulation of JEV.

The overall exposure of the human population to JEV is probably underestimated and the clinical incidence of JE in Cambodia remains difficult to assess [1,7], although this would help to promote vaccination campaigns. Measuring the exposure of dogs to JEV could be an alternative solution. Indeed, exposure of dogs to JEV has been shown in Singapore and Japan (with JEV neutralizing antibodies detected in 17% to 39.6% of dogs sampled) [32,42]. In Kandal province, the dog-to-human ratio has been estimated to 1:3.8 in 2017 [43]. A recent serological survey carried out in the same area showed that 35% of sampled dogs had JEV neutralizing antibodies [16], and blood meal analysis of *Culex spp.* trapped in this area confirmed that JEV vectors can feed on dogs [26]. As dogs live in close proximity to humans, their exposure to JEV could be a good proxy for human exposure.

The objectives of this study were (i) to evaluate $R_0$ and predict, based on simulations, human exposure to JEV in villages of Kandal, a rural province of Cambodia, (ii) to analyse how, in such villages, host community composition affects $R_0$ value and human exposure, and (iii) to assess the potential use of dogs as sentinels of human exposure to JEV.

# Materials and methods

## Study area, host populations and serological data

The study focuses on three rural districts of Kandal province numbered 1–3 below (D1: Khsach Kandal, D2: Kien Svay and D3: Kaoh Thum), located northeast, east and southeast of Phnom Penh respectively. In these 3 districts, a JEV serological survey was carried out in 2018 to assess the relative exposure to JEV of pigs, ducks, chickens, and dogs [16]. During this survey, the number of animals and humans per household was noted for each backyard where animals were sampled. To calculate the average population size per animal species and per village in each of the three districts, we used the number of households per village provided by the communal offices of 9 of the 20 villages visited. According to these village authorities, on average, 20% of the households raised pigs and 90% raised chickens. The percentage of households breeding ducks (30%) was computed based on our field observations, by calculating a ratio of households breeding ducks to households breeding chickens (Table 1). In the three districts of the study area, 112 pigs, 185 chickens, 128 ducks and 188 dogs from 20 villages were blood-sampled between March and December 2018. Samples were analysed by hemagglutination inhibition and virus neutralization assays. Depending on the district, neutralizing
<table><caption>Table 1. Average host population sizes in a traditional village of the three studied districts.</caption><thead><tr><th>District</th><th>Pigs</th><th>Ducks</th><th>Chicken</th><th>Dogs</th><th>Cattle</th><th>Humans</th></tr></thead><tbody><tr><td>D1</td><td>516.4</td><td>2020</td><td>7330</td><td>232</td><td>53.3</td><td>1630</td></tr><tr><td>D2</td><td>3644</td><td>509</td><td>14766</td><td>1460</td><td>39.2</td><td>4062</td></tr><tr><td>D3</td><td>3275</td><td>683</td><td>10312</td><td>1203</td><td>156</td><td>3880</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0010572.t001

antibodies were found in 0 to 42% of the pigs, 0 to 33% of the ducks, 0 to 0.02% of the chickens, and 35% of the dogs [16].

## Model description

**Structure of the model.** We developed a model of JEV transmission between competent hosts (pigs, ducks, and chickens) and vectors (*Culex spp.*), and from vectors to non-competent hosts (cattle, humans, and dogs) at the village scale. This model was built to allow computing $R_0$ and to analyse the epidemiological equilibrium situation. For this reason, and since previous studies conducted in the study area had shown high abundance of JEV vectors [40,44,45] and an intense circulation of the virus (seroconversion before 4 months of age in sentinel pigs, on average) [44,45], we chose to implement the model as a deterministic system, using ordinary differential equations. A detailed model description is given in S1 File. The epidemiological system represented a rural village of Kandal province where swine, ducks, chickens, cattle, dogs and humans are living. Because of different life expectancy, swine were split into two compartments: fattening pigs and sows. In the study area, these hosts are exposed to bites from *Culex tritaeniorhynchus*, *Culex vishnui* and *Culex gelidus*, the most abundant vectors of JEV in this region, and the feeding preferences of these vectors are similar [26]. For this reason, and due to the small surface area of villages in the study area (3–4 km²) compared to the flight range of *Culex spp.* [46], these different vector species were represented in the model by a single population of *Culex spp.* We assumed that direct pig-to-pig transmission and disease-induced death were negligible [47].

In our study area, pigs are raised in open pens, and poultry are often free roaming in the gardens during the day and locked in pens at night. Mosquito nets are rarely used in animal pens. Based on this livestock organization, we assumed that all hosts, including humans, were homogeneously exposed to the vector population. For this reason, and because direct transmission was neglected, it was not necessary to structure the host population by households.

We assumed that there was little or no exchange of infected hosts or vectors between the modelled and neighboring villages. Movements of live animals between villages are indeed extremely rare [37,41]. Only pigs may sometimes be transported to a slaughterhouse in a nearby village where they are slaughtered within two days. We assumed constant host population sizes, with births offsetting deaths or slaughters. The model incorporated seasonal variations in vector population size according to sinusoidal dynamics. A parameter denoted $\psi$ (with $0 \le \psi \le 1$) represented the amplitude of these variations, relative to the yearly average value $N_v$: vector population size could then vary between $N_v(1-\psi)$ in the middle of the dry season, and $N_v(1+\psi)$ in the middle of the rainy season, with $\psi = 0$ denoting a constant size of vector population throughout the year (see S1 File).

Fig 1 shows the flowchart of the model. Hosts were categorized into 5 infection states: protected by their maternal antibodies (*M*), susceptible (*S*), latent (*E*), infectious (*I*) and recovered (*R*). Hosts born from recovered females were assumed to be protected by maternal antibodies before becoming susceptible. After being bitten by an infectious mosquito, hosts had a species-specific probability (*p*) of entering the latent state (*E*). They then entered the infectious state
<figure><img src="image_4.png" /><figcaption>**Fig 1. Schematic flowchart of the JEV model.** The boxes represent the health states of hosts (h) and vectors (v): under maternal immunity (M), susceptible (S), latent (E), infected (I) and recovered (R). Solid-line arrows depict the flow into and out of compartments, determined by the associated parameters: $\mu =$ mortality rate, $\delta = 1$/duration of maternal immunity, $\Phi = 1$/duration of incubation period, $\gamma = 1$/duration of viremic period. $N_v$ is the total population of vectors. Dashed arrows depict the forces of infection (foi) exerted by an infected vector on a susceptible host (foi$_h$) and by infectious competent hosts on a susceptible vector (foi$_v$). The force of infection exerted on a vector by a non-competent host is null as the probability of being infected (q) in this case is zero.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0010572.g001

after an incubation period ($1/\Phi_h$). After the infectious (viremic) period ($1/\gamma_h$), hosts entered the recovered state (R) where they stayed until death ($1/\mu_h$). Although JEV vertical transmission has been documented in *Culex spp.*, this phenomenon appears to be limited [48,49]. Moreover, the vector population and intensity of JEV transmission in the study area are such that we assumed this vertical transmission of JEV to be negligible: all vectors emerged in the susceptible state. Susceptible vectors had a given probability ($q$) of being infected and entering the latent state when biting an infectious host, $q$ being null for non-competent hosts. After the extrinsic incubation period (of duration $1/\Phi_v$), vectors became infectious and remained in the infectious state until death.

**Forces of infection.** The force of infection in hosts, $foi_i(t)$, is the instantaneous rate at time $t$ that a susceptible host of species $i$ becomes infected. It represents the transition rate between the susceptible state (S) and the latent state (E) and was defined in our model by:

$$foi_i(t) = bx_i p_i \frac{I_v(t)}{N_i} \qquad (1)$$

where:

* $b$ was the biting rate of the vector (number of bites per vector per day), or the inverse of the duration of its gonotrophic cycle,

* $p_i$ was the transmission probability per bite from the vector to a host of species $i$,

* $I_v(t)$ was the number of infected vectors at time $t$, and $N_i$ the number of hosts of species $i$.

* $x_i$ was the probability that a vector will choose a host of species $i$ for its blood meal, calculated by multiplying the preference $\pi_i$ for host $i$ by the number of hosts of species $i$, divided by the sum of all preferences times host population sizes:

$$x_i = \frac{\pi_i N_i}{\sum_j \pi_j N_j} \qquad (2)$$
$x_i p_i$ was then the host-specific per bite transmission rate from the vector to host $i$, defined as the fraction of successful transmission events from one infected vector to a susceptible host of species $i$, per bite [50].

The global force of infection on vectors at time $t$ was the sum of the forces exerted by the different host species at time $t$:

$$f_{oi_v}(t) = b \sum_i x_i q_i \frac{I_i(t)}{N_i} \qquad (3)$$

Where $q_i$ was the probability that a susceptible mosquito became infected after a blood meal taken on a viremic host of species $i$, the value of this parameter being > 0 for the competent host species (pigs, sows, ducks and chickens) and null for the non-competent host species (humans, dogs, cattle).

**Parameterization.** The definitions and values of the model parameters are provided in Table 2, and details of how fixed parameters were obtained are given in S1 File

We estimated three of the model parameters which were unavailable in the literature, using seroprevalence data collected in the same area in pigs, poultry and dogs: the feeding preferences of *Culex spp.* for ducks ($\pi_d$) and dogs ($\pi_{dog}$), relative to pigs (used as a reference), as well as the number of vectors involved in JEV transmission. This latter parameter was separately estimated for each district ($N_{v1}$, $N_{v2}$, $N_{v3}$) to control for differences between districts that would not be captured by the model, such as inter-district variations in vector breeding sites densities or in the accessibility of vectors to hosts. We used the Nelder-Mead optimization algorithm to minimize the negative log-likelihood of our serological data (number of positive among the tested animals in each host group and district), with respect to $N_{v1}$, $N_{v2}$, $N_{v3}$, $\pi_d$ and $\pi_d$ (see S1 File for the definition of the likelihood function).

We used the “optim” function of R (version 4.0.2) [73] to minimize the log-likelihood. The variance-covariance matrix, obtained by inverting the Hessian matrix, was used to compute the confidence intervals of the estimated parameters.

## Model exploitation

**Evaluation of R₀ and human exposure in the study area.** We first used the model to characterize the circulation of JEV in villages in each of the three study districts, simulating human exposure and evaluating the basic reproduction number ($R_0$). The details of its computation are provided in S1 File. We then used three indicators computed on the 30th simulated year to predict human exposure in the studied districts, assuming the epidemiological system was in a steady state:

* the annual probability of human exposure, which is the probability of receiving at least one infective bite over 1 year. This indicator was given by: $1 - \exp(-\sum_{t=tmax-365}^{tmax} f_{oi_h}(t))$, with $tmax$ the last day of simulation (10,958th day, see S1 File) and $f_{oi_h}(t)$ the value of the force of infection in humans at time $t$;

* the average age at infection, in years. This indicator was calculated as $1 / \sum_{t=tmax-365}^{tmax} f_{oi_h}(t)$, assuming that people living in the modeled village are equally exposed to JEV since birth (this indicator was calculated only in situations where $R_0$ was above 1);

* the annual incidence rate of human infections, which represents the proportion of persons infected in one year, at the epidemiological equilibrium. This indicator was given by
<table><caption>Table 2. Parameters of the deterministic model.</caption><thead><tr><th>Parameter</th><th>Definition</th><th>Value</th><th>References</th></tr></thead><tbody><tr><td colspan="4">Hosts</td></tr><tr><td>1/μ<sub>p</sub>, 1/μ<sub>c</sub></td><td>Average lifespan of pigs &amp; chickens</td><td>6m<sup>*</sup></td><td>Field data</td></tr><tr><td>1/μ<sub>s</sub></td><td>Average lifespan of sows</td><td>3y</td><td>Field data</td></tr><tr><td>1/μ<sub>d</sub></td><td>Average lifespan of ducks</td><td>2y</td><td>Field data</td></tr><tr><td>1/μ<sub>h</sub></td><td>Average lifespan of humans</td><td>70y</td><td>[51]</td></tr><tr><td>1/μ<sub>b</sub></td><td>Average lifespan of cattle</td><td>7y</td><td>Field data</td></tr><tr><td>1/μ<sub>dog</sub></td><td>Average lifespan of dogs</td><td>5y</td><td>[43]</td></tr><tr><td>1/δ<sub>p</sub></td><td>Duration of maternal immunity in pigs &amp; sows</td><td>2.5m</td><td>[44,45]</td></tr><tr><td>1/δ<sub>d</sub>, 1/δ<sub>c</sub></td><td>Duration of maternal immunity in ducks &amp; chickens</td><td>1m</td><td>[52,53]</td></tr><tr><td>1/δ<sub>h</sub></td><td>Duration of maternal immunity in humans</td><td>5m</td><td>[54,55]</td></tr><tr><td>1/δ<sub>b</sub>, 1/δ<sub>dog</sub></td><td>Duration of maternal immunity in cattle &amp; dogs</td><td>3m</td><td>[56,57]</td></tr><tr><td>1/Φ<sub>p</sub></td><td>Incubation period in pigs &amp; sows</td><td>2d</td><td>[58–60]</td></tr><tr><td>1/Φ<sub>d</sub></td><td>Incubation period in ducks</td><td>2d</td><td>[19–22]</td></tr><tr><td>1/Φ<sub>c</sub></td><td>Incubation period in chickens</td><td>1.5d</td><td>[19,61]</td></tr><tr><td>1/Φ<sub>h</sub></td><td>Incubation period in humans</td><td>10d</td><td>[38]</td></tr><tr><td>1/Φ<sub>b</sub>, 1/Φ<sub>dog</sub></td><td>Incubation period in cattle &amp; dogs</td><td>4d</td><td>[62,63]</td></tr><tr><td>1/γ<sub>p</sub></td><td>Viremic period in pigs &amp; sows</td><td>1.5d</td><td>[58,59]</td></tr><tr><td>1/γ<sub>d</sub>, 1/γ<sub>c</sub></td><td>Viremic period in ducks &amp; chickens</td><td>3d</td><td>[19–21,61]</td></tr><tr><td>1/γ<sub>h</sub>, 1/γ<sub>b</sub>, 1/γ<sub>dog</sub></td><td>Recovery period in humans &amp; cattle &amp; dogs</td><td>5d</td><td>[62–65]</td></tr><tr><td colspan="4">Vectors</td></tr><tr><td>1/μ<sub>v</sub></td><td>Average lifespan of Cx. spp.</td><td>25d</td><td>[66]</td></tr><tr><td>b</td><td>Cx. spp. biting rate</td><td>0.25</td><td>[67,68]</td></tr><tr><td>1/Φ<sub>v</sub></td><td>Extrinsic incubation period in Cx. spp.</td><td>10d</td><td>[21,69]</td></tr><tr><td>ψ</td><td>Seasonal variations of vector population size</td><td>0</td><td>[40,44]</td></tr><tr><td colspan="4">Vector/Host interactions</td></tr><tr><td>π<sub>p</sub></td><td>Feeding preference of Cx. spp. for pigs &amp; sows</td><td>1<sup>**</sup></td><td>[26]</td></tr><tr><td>π<sub>d</sub></td><td>Feeding preference of Cx. spp. for ducks<sup>**</sup></td><td></td><td>Estimated</td></tr><tr><td>π<sub>c</sub></td><td>Feeding preference of Cx. spp. for chickens<sup>**</sup></td><td>0.09</td><td>[26]</td></tr><tr><td>π<sub>h</sub></td><td>Feeding preference of Cx. spp. for humans<sup>**</sup></td><td>0.5</td><td>[26]</td></tr><tr><td>π<sub>b</sub></td><td>Feeding preference of Cx. spp. for cattle<sup>**</sup></td><td>1.7</td><td>[26]</td></tr><tr><td>π<sub>dog</sub></td><td>Feeding preference of Cx. spp. for dogs<sup>**</sup></td><td></td><td>Estimated</td></tr><tr><td>p</td><td>Vector (Cx. spp.) to host transmission probability</td><td>0.5</td><td>[70,20]</td></tr><tr><td>q<sub>p</sub>, q<sub>d</sub>, q<sub>c</sub></td><td>Competent host to vector (Cx. spp.) transmission probability</td><td>0.5</td><td>[21,69–72]</td></tr><tr><td>q<sub>b</sub>, q<sub>dog</sub>, q<sub>h</sub></td><td>Non-competent host to vector (Cx. spp.) transmission probability</td><td>0</td><td>[62,63]</td></tr></tbody></table>

Subscripts: p = pigs, s = sows, d = ducks, c = chickens, b = cattle, h = human.

* y = year; m = month; d = day.

** Feeding preferences relative to pig, used as a reference.

https://doi.org/10.1371/journal.pntd.0010572.t002

$\frac{\sum_{t=t_{max}-365}^{t_{max}} f_{oh}(t)S_h(t)}{N_h}$, with $S_h(t)$ the number of susceptible people at day t, and $N_h$ the total number or persons in the village.

The 95% confidence intervals (CIs) of $R_0$ and of the three above indicators were calculated by drawing 1000 joint values of $N_{v1}$, $N_{v2}$, $N_{v3}$, $\pi_d$, and $\pi_{dog}$, from a multidimensional Gaussian distribution centered on the estimated parameter values, and whose dispersion parameters were given by the variance-covariance matrix (produced by parameter estimation procedure). We ran the model using these 1000 joint values, computed $R_0$ and human exposure indicators:
CI bounds were the 2.5% and 97.5% percentiles of the resulting distribution. The model was implemented using R (version 4.0.2) [73].

**Influence of host community composition on R₀ and human exposure.** In a second step, we used the model to analyse how the composition of the host community affected the ability of the epidemiological system to be invaded by JEV and influenced human exposure. To assess the effect of the host community composition variation on the estimated indicators, the size of the vector population had to remain constant in order to distinguish the effect of a change in host community composition from that of a change in the number of vectors. Second, variations in host community composition had to be realistic: since most animals are raised for self-consumption in the study area, the host community had to vary while maintaining a roughly constant level of food resources for villagers. For *Culicoides* midges, vector population size has been shown to depend on the total available body surface area (BSA) of hosts [74]. Assuming the same for *Culex spp.*, varying the host community composition while keeping their total BSA equal to a reference value satisfied the first constraint. It also satisfied the second constraint, as BSA is linked to body mass by an allometric relationship. We thus varied the host community composition by changing the proportion of each host group in a constant total BSA. We analysed the effects on R₀ and the three indicators of human exposure defined above (calculated from the simulations at epidemiological equilibrium), considering 3 types of variations of the host community composition:

* variation 1: relative share of competent hosts BSA, versus non-competent hosts BSA (percentage of competent hosts BSA ranging from 5% to 95% of the whole system BSA). As non-competent hosts “dilute” mosquito bites, we expected an impact on the proportion of infectious vectors;

* variation 2: relative share of pigs, chickens and ducks BSA among competent hosts BSA was also expected to influence JEV transmission dynamics, as each competent host has species-specific infection parameters (e.g. duration of viraemia), and vector feeding preferences differ from one species to another. For this variation, we fixed the percentage of chickens among poultry BSA at 5% (scenario A), 55% (scenario B) and 95% (scenario C) and made vary the proportion of pigs among competent hosts BSA.

* variation 3: relative share of cattle BSA among cattle-and-pigs BSA together (percentage of cattle BSA ranging from 0% to 100% of the cattle-and-pigs BSA). Cattle, being non-competent hosts on which *Culex spp.* feed, can dilute infecting bites and decrease R₀. This variation was intended to explore this dilution effect on human exposure, the human population size being constant (contrary to variation 1).

For each of the three variations, the total reference BSA (held constant) was calculated based on the average number of hosts per category (Table 1) in villages of the 3 districts (using as weights the number of villages visited during the field survey): 3420 people, 2446 pigs, 285 sows, 941 ducks, 10639 chickens, 107 cattle and 1044 dogs. The BSA of a typical individual was calculated using existing allometric formulas [75–77]: 1.53 m² for a pig, 3.47 m² for a sow, 0.15 m² for a duck, 0.13 m² for a chicken, 1.81 m² for a human, 4.29 m² for a cattle, and 0.61 m² for a dog.

The constant size of the vector population was the average of the three district-specific estimates of N<sub>v1</sub>, N<sub>v2</sub> and N<sub>v3</sub>.

**Use of dogs as sentinels of human exposure.** Finally, we used the model to analyse the relationship between JE seroprevalence in dogs and the annual probability of exposure in humans. The relationship between dog seroprevalence and annual probability of human exposure may be influenced by (i) the proportion of each host species in the multi-host system
(including the human/dog ratio), (ii) the average vector population size ($N_v$), and (iii) the life-
span of dogs (linked to the mortality rate $\mu_{dog}$) as the probability of being seropositive increases
with age in the case of an endemic disease and when antibodies remain all life long after a sin-
gle infection as it is the case with Flaviruses. Therefore, we ran 10,000 simulations for which
the composition of the host community, $N_v$, and the lifespan of dogs varied randomly, and $R_0$
was >1. For each simulation, the proportion of BSA of each host species was randomly drawn
between 0.01 and 0.99 while ensuring that the overall BSA of the system remained constant. $N_v$
was randomly drawn between 1,000 and 200,000 vectors, and the lifespan of the dogs was ran-
domly drawn between 2 and 11 years, with their average observed lifespan being 5 years [43].

We then plotted the annual probability of human exposure to JEV as a function of JEV seroprevalence in dogs and visually analysed the relationship between these two indicators.

Sensitivity analysis

To rank the model parameters according to their influence on R₀, we first performed a semi-
quantitative sensitivity analysis, using the Morris method [78]. Details are available in S2 File.

We then focused the sensitivity analysis on specific parameters, (i) for which hypotheses
had been made when parameterizing the model, *i.e.* the seasonal variations of vector popula-
tion size *N*<sub>v</sub> (*ψ*), (ii) for which data available in literature were highly variable, *i.e.* the transmis-
sion probability from viremic hosts to vectors (*q*<sub>p</sub>, *q*<sub>d</sub>, *q*<sub>c</sub>), or (iii) for which the available data
may not represent the field condition, *i.e.* the feeding preferences of *Culex* spp. for humans
(*π*<sub>h</sub>). Following a “one at a time” plan and for the traditional villages of the three studied dis-
tricts (Table 1) as well as for the three types of variations of the host community composition,
we successively considered alternative values of these parameters. For each of these values, we
re-estimated N<sub>v1</sub>, N<sub>v2</sub>, N<sub>v3</sub>, π<sub>duck</sub>, and π<sub>dog</sub>, computed R<sub>0</sub> and the indicators of human exposure
when the epidemiological system has reached a steady state. Two alternative values were con-
sidered for q<sub>p</sub>, q<sub>d</sub>, and q<sub>c</sub>: 90% lower and 90% higher. A decrease from 10% to 50% was consid-
ered for π<sub>h</sub>. Finally, and to be consistent with field observations [40,44], we considered a
seasonal variation of the vector population size (ψ) of +/-20% relative to the average annual
value.

Results

Parameter estimation and model fit

The estimated average size of the vector population was 17,789 mosquitoes for a traditional vil-
lage of district D1 (N<sub>v1</sub>, 95% CI: 12,304–25,723), 52,353 mosquitoes for district D2 (N<sub>v2</sub>, 95%
CI: 44,276–61,905), and 59,536 mosquitoes for district D3 (N<sub>v3</sub>, 95% CI: 49,993–70,902) (see
Table 1 for average host population sizes in these villages). The estimated feeding preference of
*Culex spp.* (relative to pigs that were used as reference) was 0.43 for ducks (π<sub>duck</sub>, 95% CI: 0.22–
0.85), and 0.12 for dogs (π<sub>dog</sub>, 95% CI: 0.08–0.21). Model fit results based on observed and pre-
dicted JEV ser prevalences are reported in Table 3.

Model exploitation

**Evaluation of R₀ and human exposure in the study area.** R₀ was above 1 in the 3 districts (Table 4). The estimated annual probability of human exposure to JEV ranged from 9% to 47%, and the estimated incidence rate of human JEV infections ranged from 1.37% to 1.46%, corresponding to 24 to 56 people infected per year per village (Table 4, population sizes are given in Table 1). The estimated average age at infection varied between 2 years old and 11 years old.
<table><caption>Table 3. Observed and predicted JEV seroprevalence per species in traditional villages of the three studied districts.</caption><thead><tr><th>District</th><th>Species</th><th>Observed seroprevalence [16]</th><th>Predicted seroprevalence</th></tr></thead><tbody><tr><td rowspan="3">D1</td><td>Pigs</td><td>0<sup>a</sup>/5<sup>b</sup></td><td>0.05<sup>c</sup> (0.0001–0.16)<sup>d</sup></td></tr><tr><td>Ducks</td><td>1/81</td><td>0.01 (0.00–0.07)</td><td>0.01 (0.00–0.05)</td></tr><tr><td>Chickens</td><td>1/82</td><td>0.01 (0.00–0.07)</td><td>0.004 (0.00–0.01)</td></tr><tr><td rowspan="3">D2</td><td>Pigs</td><td>15/59</td><td>0.25 (0.15–0.38)</td><td>0.25 (0.15–0.36)</td></tr><tr><td>Ducks</td><td>0/5</td><td>0.00 (0.00–0.52)</td><td>0.07 (0.03–0.14)</td></tr><tr><td>Chickens</td><td>1/46</td><td>0.02 (0.00–0.12)</td><td>0.01 (0.007–0.02)</td></tr><tr><td rowspan="4">D3</td><td>Pigs</td><td>20/48</td><td>0.42 (0.28–0.57)</td><td>0.40 (0.27–0.52)</td></tr><tr><td>Ducks</td><td>14/42</td><td>0.33 (0.20–0.50)</td><td>0.28 (0.18–0.40)</td></tr><tr><td>Chickens</td><td>0/57</td><td>0.00 (0.00–0.06)</td><td>0.02 (0.02–0.04)</td></tr><tr><td>Dogs</td><td>65/188</td><td>0.35 (0.28–0.42)</td><td>0.32 (0.26–0.38)</td></tr></tbody></table>

<sup>a</sup>Number of seropositive animals.

<sup>b</sup>Number of tested animals.

<sup>c</sup>Seroprevalence rate.

<sup>d</sup>95% confidence interval

https://doi.org/10.1371/journal.pntd.0010572.t003

## Influence of host community composition on R₀ and human exposure. Variation 1: Influence of the relative share of competent hosts BSA versus non-competent hosts BSA.

In an average village of Kandal province, competent hosts represented 46.1% of the total BSA of the village (Fig 2, dashed line). As *Culex spp.* had lower feeding preferences for competent hosts (pigs, poultry) than for cattle, which are non-competent hosts, a sufficient quantity of competent hosts was required to allow JEV to invade the epidemiological system. R₀ became greater than 1 when the percentage of competent hosts BSA reached 15% of the whole system’s BSA. R₀ then increased with the proportion of competent hosts and stabilized at about 1.35. When R₀ was >1, the annual probability of human exposure to JEV at the equilibrium varied from 0.024 to 0.45 and the annual incidence rate of human JEV infections varied from 0.012 to 0.014.

We compared two contrasted host community compositions, highlighted in Fig 2 (dotted lines), the first with 15% of competent hosts BSA, and the second with 85% of competent hosts BSA. In the first host community composition (left dotted line), an individual had an annual probability of exposure to infection of 0.02. The predicted annual incidence rate of human infections was 0.0122, corresponding to 72.4 infections among 5939 people. In the second host community composition, the annual incidence rate of human infections was similar (0.0139), although the predicted number of infections was lower than in the first situation (13 infections) due to the lower population size (952 people). Conversely, the annual exposure probability was much higher in the second than in the first host community composition, i.e. 0.44.

<table><caption>Table 4. Estimated values of R₀ and of human exposure indicators in traditional villages of the three studied districts.</caption><thead><tr><th>District</th><th>R<sub>0</sub> (95% CI)</th><th>Annual probability of human exposure to JEV (95% CI)</th><th>Annual incidence rate of human JEV infections (95% CI)</th><th>Annual incidence of human JEV infections (95% CI)</th><th>Average age at infection in humans (years) (95% CI)</th></tr></thead><tbody><tr><td>D1</td><td>1.07 (0.996–1.20)</td><td>0.09 (0.0002–0.3)</td><td>0.0146 (0.0002–0.015)</td><td>23.7 (0.3–25.1)</td><td>10.6 (3.1–>100)</td></tr><tr><td>D2</td><td>1.25 (1.16–1.37)</td><td>0.34 (0.21–0.47)</td><td>0.0137 (0.0135–0.0139)</td><td>55.8 (54.9–56.4)</td><td>2.9 (2.0–4.6)</td></tr><tr><td>D3</td><td>1.38 (1.29–1.53)</td><td>0.47 (0.34–0.62)</td><td>0.0139 (0.0137–0.0140)</td><td>53.9 (53.3–54.3)</td><td>2.0 (1.5–2.8)</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0010572.t004
<figure><img src="image_5.png" /><figcaption>Fig 2. R₀, annual probability of human exposure to JEV (blue solid line), proportion of susceptible humans (blue dashed line), annual incidence rate of human JEV infections (red line), and average age at infection, at the epidemiological equilibrium, according to the proportion of competent hosts in the whole system (BSA). The vertical dashed line corresponds to the host community composition observed in an average village of Kandal province. The two vertical dotted lines correspond to two contrasted host community compositions: one with 15% and the second with 85% of competent hosts BSA. Total BSA of hosts and the size of vector population (48,663 mosquitoes, 95% CI: 40,347–58,863) remain constant.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0010572.g002
Interestingly, the annual incidence rate of human infections was quite stable although the probability of exposure increased with the proportion of competent hosts in the system. Assuming homogeneous exposure of all age groups, the predicted average age at infection was 40 years old in the first host community composition (left dotted line) where the probability of exposure was low. In the second situation (right dotted line), where the probability of exposure was high, the predicted average age at infection was 2 years old. In the first host community composition, the annual incidence of infections was therefore distributed among all age groups, whereas in the second host community composition, infections were concentrated on the youngest age groups.

**Variation 2: Influence of the relative share of pigs, chickens and ducks BSA among competent hosts BSA.** In an average village of the study area, pigs BSA represented 35% of the whole village's BSA and 76% of the competent hosts BSA. Chickens BSA represented 91% of the poultry (chickens and ducks) BSA (Fig 3).

<figure><img src="image_6.png" /><figcaption>Fig 3. $R_0$, annual probability of human exposure to JEV (blue solid line), proportion of susceptible humans (dashed blue line), annual incidence rate of human JEV infection (red line), and average age at infection, at the epidemiological equilibrium, according to the proportion of pigs in competent hosts (BSA), for three fixed percentages of chickens among poultry (BSA): 5% (scenario A), 55% (scenario B) and 95% (scenario C). Total BSA of hosts and the size of vector population (48,663 mosquitoes, 95% CI: 40,347–58,863) remain constant.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0010572.g003
Due to the strong feeding preference of vectors for pigs and sows (lower than for cattle but twice higher than for ducks and 10 times higher than for chicken, Table 2), R₀ was mainly influenced by the proportion of pigs BSA, and it increased almost linearly with this proportion, regardless of the percentage of chickens among poultry BSA. In scenario A (percentage of chicken in poultry BSA set to 5%), R₀ was >1 when there was more than 30% of pigs in competent hosts BSA. In scenario B (percentage of chicken in poultry BSA set to 55%), R₀ was >1 with only 10% of pigs in competent hosts BSA. At first, R₀ slightly increased also with the percentage of chickens in poultry BSA because their lifespan (6 months) is shorter than that of ducks (2 years): a short lifespan ensures a faster turnover of susceptible animals in the system. However, when the percentage of chicken among competent hosts was very high (95%, scenario C), R₀ became >1 with at least 35% of pigs in competent hosts BSA. Indeed, the mosquitoes' feeding preference for chickens is the lowest among the six studied species ($π_c$ = 0.09). When there were almost only chickens among the competent hosts, vectors bit the non-competent hosts more frequently, resulting in a decrease of R₀.

Starting from 25% of chickens in poultry BSA (S1 Fig), and with 55% of chickens in poultry BSA (Fig 3, scenario B), and no pigs, R₀ was close to 1 and its 95% confidence interval encompassed 1.

When R₀ was >1, the annual probability of human exposure to JEV at the equilibrium varied from 0.004 to 0.58 and the annual incidence rate of human JEV infections varied from 0.003 to 0.015. In scenario B, 20% of pigs in competent hosts BSA led to a predicted annual probability of exposure of 0.05 and a predicted annual incidence rate of 0.0155 (corresponding to 53.2 infections per year per village). In the same scenario but with 80% of pigs in competent hosts BSA, the annual exposure probability reached 0.40, and the annual incidence rate was similar: 0.0138 (corresponding to 47.2 infections per year per village).

Again, the annual incidence rate of infection varied slightly whereas the annual probability of exposure increased with the proportion of pigs in the competent hosts BSA. As the latter increased, the average age at infection decreased and infections were more likely concentrated in the younger age groups. In scenario B, the average age at infection was over 70 years old (the graphs have been truncated above the lifespan of humans), over 20 years old and less than 5 years old with respectively 10%, 20% and more than 50% of pigs in competent hosts BSA.

**Variation 3: Influence of the relative share of cattle BSA among cattle-and-pigs BSA.** In an average village of Kandal province, cattle BSA represented 8.9% of the total cattle-and-pigs BSA. Due to the strong feeding preference of vectors for cattle (the highest of all competent and non-competent hosts, Table 2), R₀ and human exposure decreased as the percentage of cattle increased (Fig 4). From 65% of cattle among the total cattle-and-pigs BSA, R₀ fell below 1. The annual incidence rate was maximal, i.e. 0.015, when there was 60% of cattle among cattle-and-pigs BSA (corresponding to 23% of cattle and 15% of pigs in the whole system's BSA). This incidence rate corresponded to 51.3 human infections per year in the village (3420 people). For this system composition, the annual probability of human exposure to JEV was 0.06 and the average age at infection was approximately 18 years old.

**Use of dogs as sentinels of human exposure to JEV.** The relationship between the predicted annual probability of human exposure and the predicted seroprevalence in dogs appeared roughly linear (Fig 5), with the slope of the curves influenced only by the average lifespan of dogs (S2 Fig). It is worth noting that, for the dog average lifespan estimated in the study area (5 years [43]), the slope of the curve was close to 1, the annual probability of human exposure to infective bites being then approximately equal to the seroprevalence in dog. This result suggests that, in the studied area, JE seroprevalence in dogs would be a good proxy for human exposure to JEV infection.
<figure><img src="image_7.png" /><figcaption>Fig 4. R₀, annual probability of human exposure to JEV (blue solid line), proportion of susceptible humans (blue dashed line), annual incidence rate of human JEV infections (red), and average age at infection, at the epidemiological equilibrium, according to the proportion of cattle among cattle and pigs (BSA). The vertical dashed line corresponds to the host community composition observed in an average village of Kandal province. Total BSA of hosts and the size of vector population (48,663 mosquitoes, 95% CI: 40,347–58,863) remain constant.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0010572.g004
<figure><img src="image_8.png" /><figcaption>Fig 5. Annual probability of human exposure to JEV according to dog JEV seroprevalence.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0010572.g005

## Sensitivity analysis

Details of Morris sensitivity analysis results are available in S2 File. Among the 51 model parameters, 9 had a significant and linear influence on R₀. As expected, these were mostly vector and pig related parameters: the vector biting rate *b*, the vector mortality rate *μ*_v, the average vector population size *N*_v, the pig recovery rate *γ*_p, the probability of JEV transmission from pigs to vectors *q*_p, the probability of JEV transmission from vectors to pigs *p*_p, and the feeding preference of vectors for pigs *π*_p.

Details of sensitivity analysis of model parameter estimates, R₀ and human exposure indicators to host-to-vector transmission probability and to *Culex spp.* feeding preference for humans, in traditional villages of the three studied districts, are given in S2 File. Large variations in the value of host-to-vector transmission probability had almost no impact on R₀ and on human exposure indicators in the host community composition observed in the field (Table A in S2 File). A 50% decrease in the vector feeding preference for humans, based on the assumption that humans may be less exposed to mosquito bites than animals living outdoors, did not influence parameter estimates, but resulted in a doubling of the average age of infection in district D1 (Table B in S2 File). Finally, incorporating seasonal variations of vector population size did not significantly influence R₀ and human exposure indicators (Table C in S2 File).

The indicators (R₀, annual exposure probability, annual incidence rate and average age at infection of humans) were then calculated for each variation in host community composition with the new input parameter values and parameter estimates (S2 File). In variation 1 (relative share of competent hosts versus non-competent hosts BSA), varying input parameters had almost no impact on output indicators, except an expected decrease of human exposure
probability when $\pi_h$ was halved. Overall, the main qualitative results obtained previously, i.e. an $R_0$ increasing with the proportion of competent hosts and a distribution of annual incidence of infections concentrated on the youngest age groups when there is a majority of competent hosts in the system, were not changed (Fig B in S2 File). In variation 2 (relative share of pigs, chickens and ducks among competent hosts), the 90% decrease in $q_p$ led to a profound change in the dynamics of the system since, on the one hand, the estimated averaged $N_v$ was much larger than in the initial situation, and on the other hand, the system converged towards a transmission dynamic where the virus circulated mainly between birds. When modifying the other input parameters ($q_p, q_c, \pi_h, \psi$), key qualitative results were not changed compared to those obtained from simulations with the default parameters (Fig C in S2 File). In variation 3 (relative share of cattle among cattle and pigs), varying input parameters had almost no impact on output indicators and key qualitative results presented above remained similar with those obtained from simulations performed with the default parameters (Fig D in S2 File).

## Discussion

Our model, calibrated on field-collected serological and demographic data first showed an intense circulation of JEV in Kandal province, with $R_0$ values ranging from 1.07 and 1.38 and an annual probability of human exposure from 9 to 47% depending on the district. The average age at infection was always low, i.e. between 2 and 11 years old suggesting an important clinical impact on children health in absence of vaccination. Secondly, the simulation results confirm previous experimental results [19–24]: poultry could serve as a reservoir and JEV could invade a system without pigs. Lastly, dogs might be a good proxy for human JEV exposure in the study area.

Estimation of the model parameters was based on the assumption that the system in which samples were collected to calculate seroprevalences was at endemic equilibrium state, which may not be accurate. However, serological studies performed in the same region suggest little seasonal and inter-annual variation in JEV circulation in pigs [44,45] and ducks [8,16], and allowed us to make this assumption.

Available data on the clinical incidence of JEV in Cambodia are scarce and come mainly from hospital-based studies, sometimes in areas where access to care structures and clinical case identification capacities are limited [33,36,79,80]. The community-level clinical JEV incidence is thus probably underestimated. Under these conditions, Tarantola et al. (2014) estimated that 1/250 to 1/500 JEV infections resulted in symptomatic cases in Cambodia [81]. Based on these ratios and the annual incidence of human JEV infections predicted by our model (between 23.7 and 56 infections per year in a village of 1630 and 4062 people respectively), one can expect 0.05 (if the rate is 1/500) to 0.09 (if the rate is 1/250) clinical cases per year per village in district D1, and 0.11 to 0.22 clinical cases per year per village in districts D2 and D3. Although this is only a rough estimate, it would correspond to a clinical incidence of 2.8 to 5.8/100,000 JEV cases per year. In comparison, Mao et al. (2020), and Tian et al. (2015), estimated the average annual incidence of JEV cases in Yunnan and Changsha provinces, China, to be 0.16/100,000 (in 2017, after 10 years of vaccination program) and 0.15/100,000 respectively, by collecting hospital data for 6 to 10 years [82,83]. After 19 years of follow-up, Montini et al. (2020) estimated an incidence of 0.16/100,000 cases per year in Malaysia [84]. In Bhutan, 0.3/100,000 and 0.8/100,000 cases per year were estimated in 2020 for adults and children respectively, based on data from 5 sentinel hospitals [85]. Finally, Campbell et al. (2011) extrapolated hospital-based data of 12 southeast Asia countries, and estimated a global incidence of 1.8/100,000 and 5.4/100,000 cases per year for adults and children respectively, in all the 24 countries where JEV circulates [7]. Although they only concern the province of Kandal
in Cambodia, our estimates are consistent with the latter figures. However, these exposure fig-
ures are to be set against the JE vaccination data in Cambodia since the consequence of this
exposure will depend on the vaccination coverage of the population. The implementation of
vaccination in Cambodia is recent, and the reported number of individuals vaccinated went
from 0 before 2010 to more than 500,000 in 2015 [1]. Even if Quan et al, 2020 estimated a con-
sequent reduction of JE cases through vaccination [1], JEV was still reported in 2017, after the
main vaccination campaigns in 2016, as the primary cause of acute meningitis-encephalitis in
children, with 35% of 1160 patients confirmed or highly probable to have JEV infection [36].

The proportion of severe cases occurring in Kandal province remains unknown. However, the predicted average age at infection in the modeled villages was very low, (ranging from 2 to 11 years old). Since young children are more likely to develop severe forms after JEV infection, as well as severe sequelae [2,33,82,84,85], these results and the current knowledge we have about JEV circulation in Cambodia would justify to intensify child immunization.

According to our results, variations of host community composition would influence the average age at infection: in peri-urban areas (with few competent hosts in the system (Fig 2) or few pigs among competent hosts (Fig 3)), infections might be concentrated in adults, whereas in rural areas (with more competent hosts in the system (Fig 2) or more pigs among competent hosts (Fig 3)), infections might be concentrated in the younger age groups. In the three studied districts, the modelled villages belonged to the latter category where the predicted average age of infection was low while the exposure probability was high.

Results of simulations confirmed the major role played by pigs in JEV circulation. However, it appeared that in a system without pigs, the upper bound of R₀ confidence interval remained >1 when there was more than 25% of chickens in poultry BSA (S1 Fig). This suggests that, in a pig-free system, chickens and ducks could be sufficient for the virus to invade the epidemiological system. As shown by the sensitivity analysis, our results appeared robust with respect to the values of uncertain parameters (qₚ, q₉, qₑ, πₕ) or for which specific assumptions had been made during model parameterization.

As Lord et al. (2015) and Bae et al. (2018) have already pointed it out, the epidemiology of JEV needs to be rethought, as JEV circulation could be maintained, or not, within complex and area-dependent epidemiological systems [13,14].

Epidemiological systems are bound to change, for political, economic, cultural or sanitary reasons. Traditional pig farming has been widespread in Cambodia since before the 1950s. The Khmer Rouge period constituted a rupture in the country's history, and pig farming resumed only from the 1980s [37,41]. Since then, traditional pig farming has been facing sanitary crises such as Classical Swine fever and recently African Swine fever (ASF) [86], which could drastically reduce swine density. Moreover, the industrialization of the sector and the competition with intensive production, notably from Thailand, is driving down the price of meat and gradually discouraging small farmers for whom traditional pig farming is no longer profitable (National Animal Health and Production Research Institute, personal communication). On the other hand, landscape management, through drainage or reduction of rice fields due to industrialization, as well as climate change and climatic hazards could also transform the epidemiological system by directly affecting host and vector populations, and thus the transmission pattern of JEV. The changes in host community composition we have simulated are thus likely to occur and the current sanitary context related to the circulation of ASF in Asia is a concrete example of these changes. In Vietnam, more than 21% of the total pig herd was decimated in 2020 due to the ASF outbreak. The reduction of the total number of pig herds has led to a rapid growth of cattle production (+5% in one year) [87]. Our model confirmed that R₀ might decrease and the level of virus circulation at endemic equilibrium be
lower in systems where the proportion of cattle increases. This zooprophylactic-like effect is explained by the fact that cows are non-competent hosts on which JEV vectors feed [26,28].

In our study area, it has been shown that dogs, which are numerous and live in close proximity to humans, are widely exposed to JEV infection [16]. In an endemic area where dogs live an average of 5 years, the annual probability of human exposure to JEV was similar to the value of seroprevalence in dogs (Fig 5). The simulations indicated that the relationship between these two indicators may be generalizable, as it was approximately linear regardless of the size and composition of the host community and the size of the vector population (Fig 5). In the particular context of Kandal province, seroprevalence in dogs may be a good proxy for human exposure, and a tool for estimating the impact of JEV on public health. If this result was confirmed in varied epidemiological contexts, the practical use of dogs as sentinels for human JEV exposure would depend on the epidemiological situation of the region. In endemic contexts as in Cambodia, a verified correlation between dog seroprevalence and human exposure probability would help quantifying people exposure, by implementing serological surveys in dogs. In epidemic contexts where JEV circulation is seasonal as in Thailand or Vietnam [39,88], detecting JEV circulation in sentinel dogs would help implementing prevention or information measures ahead of the expected waves of exposure. In disease-free areas, sentinel dogs could be used as an early-surveillance system of JEV emergence in risky areas. Exposure data from sentinel dogs can also be used to target vaccination to areas where expected human exposure is greatest or where access to JEV vaccine or resources to implement vaccination are limited. Even if in-depth surveys under various environmental conditions should be further implemented to infer the potential use of dogs as JEV sentinels in the future, this complements the results of studies suggesting that dogs could be used as sentinels for other flaviviruses such as WNV [89,90].

Our theoretical approach showed that variations of the composition of the multi-host system identified in Cambodia may have an impact on the ability of the epidemiological system to sustain JEV transmission, on the human exposure to JEV, and thus on the disease burden in humans, especially in young children. Besides children vaccination in JEV endemic areas, a proper evaluation of the impact on human health is needed as well as further investigation on the potential use of dog as sentinels of human exposure, to target prevention actions and reduce JEV burden in Cambodia.

## Supporting information

S1 File. Equations of the deterministic compartmental model (I), Model parametrization and likelihood function (II) and Next Generation Matrix method for R₀ computation (III).

(DOCX)

S2 File. Morris sensitivity analysis of R₀ to parameter values (I), sensitivity analysis of model parameter estimates, R₀ and human exposure indicators to host-to-vector transmission probability (*q*<sub>p</sub>, *q*<sub>d</sub>, *q*<sub>c</sub>), to *Culex spp.* feeding preference for humans (*π*<sub>h</sub>), and to +/-20% seasonal variations of vector population size (*ψ* = 0.2) in traditional villages of the three studied districts (II) and sensitivity analysis of R₀ and human exposure indicators to *q*<sub>p</sub>, *q*<sub>d</sub>, *q*<sub>c</sub>, *π*<sub>h</sub>, and *ψ*, for the 3 variations of host community composition (III).

(DOCX)

S1 Fig. Variations of R₀ according to the proportion of chickens among poultry (BSA), and for five percentages of pigs among competent hosts (BSA). The black point corresponds to
the R₀ for 100% of pigs in competent hosts BSA.
(TIF)

S2 Fig. Annual probability of human exposure to JEV according to dog JEV seroprevalence, for varying human-dog ratios and vector population sizes.
(TIF)

**Acknowledgments**

The authors would like to thank Dr Didier Fontenille (Head of KIM RIVE-Défi RIVOC, MIVEGEC Unit, IRD), for supporting and enabling this work, Dr Philippe Dussart (Director of Institut Pasteur de Madagascar), Dr Veasna Duong (Head of the Virology Unit of Institut Pasteur du Cambodge), Dr Sebastien Boyer (Head of the Medical and Veterinary Entomology Unit of Institut Pasteur du Cambodge) for their support and collaboration, Dr Heidi Auerswald (senior researcher in virology in Institut Pasteur du Cambodge) for her help in interpreting serological test results and the Cambodian National Animal Health and Production Research Institute for their valuable help in collecting and analysing the field data.

**Author Contributions**

**Conceptualization:** Héléna Ladreyt, Véronique Chevalier, Benoit Durand.

**Data curation:** Héléna Ladreyt, Véronique Chevalier, Benoit Durand.

**Formal analysis:** Héléna Ladreyt, Véronique Chevalier, Benoit Durand.

**Funding acquisition:** Héléna Ladreyt, Véronique Chevalier, Benoit Durand.

**Investigation:** Héléna Ladreyt.

**Methodology:** Héléna Ladreyt, Véronique Chevalier, Benoit Durand.

**Project administration:** Véronique Chevalier, Benoit Durand.

**Resources:** Véronique Chevalier, Benoit Durand.

**Software:** Héléna Ladreyt, Benoit Durand.

**Supervision:** Véronique Chevalier, Benoit Durand.

**Validation:** Véronique Chevalier, Benoit Durand.

**Visualization:** Héléna Ladreyt.

**Writing – original draft:** Héléna Ladreyt.

**Writing – review & editing:** Héléna Ladreyt, Véronique Chevalier, Benoit Durand.

**References**

1. Quan TM, Thao TTN, Duy NM, Nhat TM, Clapham H. Estimates of the global burden of Japanese encephalitis and the impact of vaccination from 2000–2015. Franco E, editor. eLife. 2020 May 26; 9: e51027. https://doi.org/10.7554/eLife.51027 PMID: 32450946

2. Mayxay M, Douangdala P, Vilayhong C, Phommasone K, Chansamouth V, Vongsouvath M, et al. Outcome of Japanese Encephalitis Virus (JEV) Infection in Pediatric and Adult Patients at Mahosot Hospi-tal, Vientiane, Lao PDR. Am J Trop Med Hyg. 2020 Dec 21; 104(2):567–575. https://doi.org/10.4269/ajtmh.20-0581 PMID: 33350379; PMCID: PMC7866311.

3. Fischer M, Hills S, Staples E, Johnson B, Yaich M, Solomon T. Japanese Encephalitis Prevention and Control: Advances, Challenges, and New Initiatives. Emerging Infections 8. 2008 Jan 1;93–124.
4. Ooi MH, Lewthwaite P, Lai BF, Mohan A, Clear D, Lim L, et al. The Epidemiology, Clinical Features, and Long-Term Prognosis of Japanese Encephalitis in Central Sarawak, Malaysia, 1997–2005. Clinical Infectious Diseases. 2008 Aug 15; 47(4):458–68. https://doi.org/10.1086/590008 PMID: 18616397

5. Ding D, Hong Z, Zhao S jun, Clemens JD, Zhou B, Wang B, et al. Long-Term Disability from Acute Childhood Japanese Encephalitis in Shanghai, China. Am J Trop Med Hyg. 2007 Sep 1; 77(3):528–33. PMID: 17827373

6. Yin Z, Wang X, Li L, Li H, Zhang X, Li J, et al. Neurological Sequelae of Hospitalized Japanese Encephalitis Cases in Gansu Province, China. Am J Trop Med Hyg. 2015 Jun 3; 92(6):1125–9. https://doi.org/10.4269/ajtmh.14-0148 PMID: 25918213

7. Campbell GL, Hills SL, Fischer M, Jacobson JA, Hoke CH, Hombach JM, et al. Estimated global incidence of Japanese encephalitis: a systematic review. Bull World Health Organ. 2011 Oct 1; 89 (10):766–74, 774A-774E. https://doi.org/10.2471/BLT.10.085233 PMID: 22084515

8. Auerswald H, Ruget AS, Ladreyt H, In S, Mao S, Sorn S, et al. Serological evidence for Japanese encephalitis and West Nile virus infections in domestic birds in Cambodia. Frontiers in Veterinary Science. 2020; 7:9.

9. Maeki T, Tajima S, Ikeda M, Kato F, Taniguchi S, Nakayama E, et al. Analysis of cross-reactivity between flaviviruses with sera of patients with Japanese encephalitis showed the importance of neutralization tests for the diagnosis of Japanese encephalitis. J Infect Chemother. 2019 Oct; 25(10):786–790. https://doi.org/10.1016/j.jiac.2019.04.003 Epub 2019 May 16. PMID: 31105002.

10. de Wispelaere M, Desprès P, Choumet V. European Aedes albopictus and Culex pipiens Are Competent Vectors for Japanese Encephalitis Virus. PLoS Negl Trop Dis. 2017 Jan 13; 11(1). https://doi.org/10.1371/journal.pntd.0005294 PMID: 28085881

11. Faizah AN, Kobayashi D, Amoa-Bosompem M, Higa Y, Tsuda Y, Itokawa K, et al. Evaluating the competence of the primary vector, Culex tritaeniorhynchus, and the invasive mosquito species, Aedes japonicus japonicus, in transmitting three Japanese encephalitis virus genotypes. PLoS Negl Trop Dis. 2020 Dec 28; 14(12):e0008986. https://doi.org/10.1371/journal.pntd.0008986 PMID: 33370301

12. Yap G, Lim XF, Chan S, How CB, Humaidi M, Yeo G, et al. Serological evidence of continued Japanese encephalitis virus transmission in Singapore nearly three decades after end of pig farming. Parasites & Vectors. 2019 May 17; 12(1):244. https://doi.org/10.1186/s13071-019-3501-0 PMID: 31101069

13. Lord JS, Gurley ES, Pulliam JRC. Rethinking Japanese Encephalitis Virus Transmission: A Framework for Implicating Host and Vector Species. PLoS Negl Trop Dis. 2015 Dec 10; 9(12):e0004074. https://doi.org/10.1371/journal.pntd.0004074 PMID: 26657648

14. Bae W, Kim JH, Kim J, Lee J, Hwang ES. Changes of epidemiological characteristics of Japanese encephalitis viral infection and birds as a potential viral transmitter in Korea. Journal of Korean Medical Science. 2018; 33(9). https://doi.org/10.3346/jkms.2018.33.e70 PMID: 29441740

15. Adi AAAM, Astawa N, Ayu Asri Damayanti P, Made Kardena I, Gusti Made Krisna Erawan I, Suardana IW, et al. Seroepidemiological Evidence for the Presence of Japanese Encephalitis Virus Infection in Ducks, Chickens, and Pigs, Bali-Indonesia. Bali Medical Journal. 2016 Nov 14; 5:189.

16. Ladreyt H, Auerswald H, Tum S, Ken S, Heng L, In S, et al. Comparison of Japanese Encephalitis Force of Infection in Pigs, Poultry and Dogs in Cambodian Villages. Pathogens. 2020 Sep; 9(9):719. https://doi.org/10.3390/pathogens9090719 PMID: 32882890

17. Kalaiyarasu S, Mishra N, Khetan RK, Singh VP. Serological evidence of widespread West Nile virus and Japanese encephalitis virus infection in native domestic ducks (Anas platyrhynchos var domesticus) in Kuttanad region, Kerala, India. Comparative Immunology, Microbiology and Infectious Diseases. 2016 Oct 1; 48:61–8. https://doi.org/10.1016/j.cimid.2016.08.002 PMID: 27638121

18. Pant GR, Lunt RA, Rootes CL, Daniels PW. Serological evidence for Japanese encephalitis and West Nile viruses in domestic animals of Nepal. Comparative Immunology, Microbiology and Infectious Diseases. 2006; 29(2–3):166–75. https://doi.org/10.1016/j.cimid.2006.03.003 PMID: 16697904

19. Cleton NB, Bosco-Lauth A, Page MJ, Bowen RA. Age-Related Susceptibility to Japanese Encephalitis Virus in Domestic Ducklings and Chicks. Am J Trop Med Hyg. 2014 Feb 5; 90(2):242–6. https://doi.org/10.4269/ajtmh.13-0161 PMID: 24394476

20. Di D, Li C, Zhang J, Hameed M, Wang X, Xia Q, et al. Experimental Infection of Newly Hatched Domestic Ducklings via Japanese Encephalitis Virus-Infected Mosquitoes. Pathogens. 2020 May; 9(5):371. https://doi.org/10.3390/pathogens9050371 PMID: 32408553

21. Karna AK, Bowen RA. Experimental Evaluation of the Role of Ecologically-Relevant Hosts and Vectors in Japanese Encephalitis Virus Genotype Displacement. Viruses. 2019 Jan; 11(1):32. https://doi.org/10.3390/v11010032 PMID: 30621345

22. Xiao C, Wang X, Cui G, Pang L, Xu J, Li C, et al. Possible pathogenicity of Japanese encephalitis virus in newly hatched domestic ducklings. Vet Microbiol. 2018 Dec; 227:8–11. https://doi.org/10.1016/j.vetmic.2018.10.016 PMID: 30473356
23. Dhanda V, Banerjee K, Deshmukh PK, Ilkal MA. Experimental viraemia and transmission of Japanese encephalitis virus by mosquitoes in domestic ducks. Indian Journal of Medical Research. 1977; 66 (6):881-8. PMID: 205503

24. Gresser I, Hardy JL, Hu SMK, Scherer WF. Factors Influencing Transmission of Japanese B Encephalitis Virus by a Colonized Strain of Culex Tritaeniorhynchus Giles, from Infected Pigs and Chicks to Susceptible Pigs and Birds. Am J Trop Med Hyg. 1958 Jul 1; 7(4):365-73. https://doi.org/10.4269/ajtmh.1958.7.365 PMID: 13559585

25. Stewart Merrill TE, Johnson PTJ. Towards a mechanistic understanding of competence: a missing link in diversity-disease research. Parasitology. 2020 Sep; 147(11):1159-70. https://doi.org/10.1017/S0031182020000943 PMID: 32517830

26. Boyer S, Durand B, Yean S, Brengues C, Maquart PO, Fontenille D, et al. Host-Feeding Preference and Diel Activity of Mosquito Vectors of the Japanese Encephalitis Virus in Rural Cambodia. Pathogens. 2021 Mar; 10(3):376. https://doi.org/10.3390/pathogens10030376 PMID: 33800999

27. Yeo G, Chan S, How CB, Humaidi M, Lim XF, Mailepessov D, et al. Molecular Analysis of the Blood-meals of Culex spp. Mosquitoes at Natural Habitats in Singapore to Investigate the Potential Risk of Japanese Encephalitis Virus and West Nile Virus Transmission. Vector-Borne and Zoonotic Diseases. 2020 Jun 16; 20(9):703-14. https://doi.org/10.1089/vbz.2019.2576 PMID: 32931404

28. Mwandawiro C, Boots M, Tuno N, Suwonkerd W, Tsuda Y, Takagi M. Heterogeneity in the host preference of Japanese encephalitis vectors in Chiang Mai, northern Thailand. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2000; 94(3):238-42. https://doi.org/10.1016/S0035-9203(00)90303-1 PMID: 10974986

29. Ladreyt H, Durand B, Dussart P, Chevalier V. How Central Is the Domestic Pig in the Epidemiological Cycle of Japanese Encephalitis Virus? A Review of Scientific Evidence and Implications for Disease Control. Viruses. 2019 Oct 15; 11(10). https://doi.org/10.3390/v11100949 PMID: 31618959

30. Jamgaonkar AV, Yergolkar PN, Geevarghese G, Joshi GD, Joshi MV, Mishra AC. Serological evidence for Japanese encephalitis virus and West Nile virus infections in water frequenting and terrestrial wild birds in Kolar District, Karnataka State, India. A retrospective study. Acta Virol. 2003; 47(3):185-8. PMID: 14658848

31. Kako N, Suzuki S, Sugie N, Kato T, Yanase T, Yamakawa M, et al. Japanese encephalitis in a 114-month-old cow: pathological investigation of the affected cow and genetic characterization of Japanese encephalitis virus isolate. BMC Vet Res. 2014 Mar 11; 10:63. https://doi.org/10.1186/1746-6148-10-63 PMID: 24618225

32. Shimoda H, Ohno Y, Mochizuki M, Iwata H, Okuda M, Maeda K. Dogs as Sentinels for Human Infection with Japanese Encephalitis Virus. Emerg Infect Dis. 2010 Jul; 16(7):1137-9. https://doi.org/10.3201/eid1607.091757 PMID: 20587189

33. Touch S, Hills S, Sokhal B, Samnang C, Sovann L, Khieu V, et al. Epidemiology and burden of disease from Japanese encephalitis in Cambodia: results from two years of sentinel surveillance. Trop Med Int Health. 2009 Nov; 14(11):1365-73. https://doi.org/10.1111/j.1365-3156.2009.02380.x PMID: 19747185

34. Gulati BR, Singha H, Singh BK, Virmani N, Kumar S, Singh RK. Isolation and genetic characterization of Japanese encephalitis virus from equines in India. J Vet Sci. 2012 Jun; 13(2):111-8. https://doi.org/10.4142/jvs.2012.13.2.111 PMID: 22705732

35. Keeling MJ, Rohani P. Modeling Infectious Diseases in Humans and Animals. Princeton University Press; 2007.

36. Horwood PF, Duong V, Laurent D, Mey C, Sothy H, Santy K, et al. Aetiology of acute meningoencephalitis in Cambodian children, 2010–2013. Emerg Microbes Infect. 2017 May 24; 6(5):e35. https://doi.org/10.1038/emi.2017.15 PMID: 28536430

37. Knips V. Review of the livestock sector in the Mekong countries. FAO, AGAL; 2004.

38. Zhao S, Lou Y, Chiu APY, He D. Modelling the skip-and-resurgence of Japanese encephalitis epidemics in Hong Kong. Journal of Theoretical Biology. 2018; 454:1–10. https://doi.org/10.1016/j.jtbi.2018.05.017 PMID: 29792875

39. Ruget AS, Beck C, Gabassi A, Trevennec K, Lecollinet S, Chevalier V, et al. Japanese encephalitis circulation pattern in swine of northern Vietnam and consequences for swine's vaccination recommendations. Transboundary and Emerging Diseases. 2018; 65(6):1485-92. https://doi.org/10.1111/tbed.12885 PMID: 29740970

40. Boyer S, Peng B, Pang S, Chevalier V, Duong V, Gorman C, et al. Dynamics and diversity of mosquito vectors of Japanese encephalitis virus in Kandal province, Cambodia. Journal of Asia-Pacific Entomology. 2020 Dec 1; 23(4):1048-54.

41. Huynh TTT, Aarnink AJA, Drucker A, Verstegen MWA. Pig Production in Cambodia, Laos, Philippines, and Vietnam: A Review. 2002; 3(1):22.
42. Ting SHL, Tan HC, Wong WK, Ng ML, Chan SH, Ooi EE. Seroepidemiology of neutralizing antibodies to Japanese encephalitis virus in Singapore: Continued transmission despite abolishment of pig farming? Acta Tropica. 2004; 92(3):187-91. https://doi.org/10.1016/j.actatropica.2004.04.010 PMID: 15533286

43. Chevalier V, Davun H, Sorn S, Ly P, Pov V, Ly S. Large scale dog population demography, dog management and bite risk factors analysis: A crucial step towards rabies control in Cambodia. PLoS ONE. 2021 Jul 8; 16(7):e0254192. https://doi.org/10.1371/journal.pone.0254192 PMID: 34237103

44. Cappelle J, Duong V, Pring L, Kong L, Yakovleff M, Prasetyo DB, et al. Intensive Circulation of Japanese Encephalitis Virus in Peri-urban Sentinel Pigs near Phnom Penh, Cambodia. PLoS Negl Trop Dis. 2016 Dec 7; 10(12):e0005149. https://doi.org/10.1371/journal.pntd.0005149 PMID: 27926937

45. Di Francesco J, Choeung R, Peng B, Pring L, Pang S, Duboz R, et al. Comparison of the dynamics of Japanese encephalitis virus circulation in sentinel pigs between a rural and a peri-urban setting in Cambodia. PLoS Negl Trop Dis. 2018; 12(8). https://doi.org/10.1371/journal.pntd.0006644 PMID: 30138381

46. Clements A. Biology of Mosquitoes, Volume 2 CABI Publishing; 1999.

47. Diallo AOI, Chevalier V, Cappelle J, Duong V, Fontenille D, Duboz R. How much does direct transmission between pigs contribute to Japanese Encephalitis virus circulation? A modelling approach in Cambodia. PLoS ONE. 2018; 13(8):e0201209. https://doi.org/10.1371/journal.pone.0201209 PMID: 30114197

48. Kumari R, Kumar K, Rawat A, Singh G, Yadav NK, Chauhan LS. First indigenous transmission of Japanese Encephalitis in urban areas of National Capital Territory of Delhi, India. Trop Med Int Health. 2013 Jun; 18(6):743-9. https://doi.org/10.1111/tmi.12104 PMID: 23682856

49. Rosen L, Lien JC, Shroyer DA, Baker RH, Lu LC. Experimental vertical transmission of Japanese encephalitis virus by Culex tritaeniorhynchus and other mosquitoes. Am J Trop Med Hyg. 1989 May; 40 (5):548-56. https://doi.org/10.4269/ajtmh.1989.40.548 PMID: 2567124

50. Fischer EAJ, Boender GJ, Nodelijk G, de Koeijer AA, van Roermund HJW. The transmission potential of Rift Valley fever virus among livestock in the Netherlands: a modelling study. Vet Res. 2013 Jul 22; 44:58. https://doi.org/10.1186/1297-9716-44-58 PMID: 23876054

51. Life expectancy at birth, total (years)—Cambodia | Data. 2019. [cited 2021 Mar 31]. Available from: https://data.worldbank.org/indicator/SP.DYN.LE00.IN?locations=KH

52. Nemeth NM, Bowen RA. Dynamics of passive immunity to West Nile virus in domestic chickens (Gallus gallus domesticus). Am J Trop Med Hyg. 2007 Feb; 76(2):310-7. PMID: 17297041

53. Buescher EL, Scherer WF, Rosenberg MZ, Kutner LJ, McClure HE. Immunologic Studies of Japanese Encephalitis Virus in Japan: IV. Maternal Antibody in Birds. The Journal of Immunology. 1959 Dec 1; 83 (6):614-9. PMID: 13805723

54. Adimy M, Mancera PFA, Rodrigues DS, Santos FLP, Ferreira CP. Maternal Passive Immunity and Dengue Hemorrhagic Fever in Infants. Bull Math Biol. 2020 Jan 28; 82(2):24. https://doi.org/10.1007/s11538-020-00699-x PMID: 31993759

55. Capeding RZ, Brion JD, Caponpon MM, Gibbons RV, Jarman RG, Yoon IK, et al. The Incidence, Characteristics, and Presentation of Dengue Virus Infections during Infancy. Am J Trop Med Hyg. 2010 Feb 5; 82(2):330-6. https://doi.org/10.4269/ajtmh.2010.09-0542 PMID: 20134013

56. Horimoto M, Sakai T. Maternally derived antibodies to Japanese encephalitis virus in cattle. Kansenshogaku Zasshi. 1990 Sep; 64(9):1205-8. https://doi.org/10.11150/kansenshogakuzasshi1970.64. 1205 PMID: 2243192

57. Chastant S, Mila H. Passive immune transfer in puppies. Anim Reprod Sci. 2019 Aug; 207:162-70. https://doi.org/10.1016/j.anireprosci.2019.06.012 PMID: 31255495

58. Williams DT, Daniels PW, Lunt RA, Wang LF, Newberry KM, Mackenzie JS. Experimental infections of pigs with Japanese encephalitis virus and closely related Australian flaviviruses. Am J Trop Med Hyg. 2001; 65(4):379-87. https://doi.org/10.4269/ajtmh.2001.65.379 PMID: 11693888

59. Ricklin ME, García-Nicolás O, Brechbühl D, Python S, Zumkehr B, Posthaus H, et al. Japanese encephalitis virus tropism in experimentally infected pigs. Veterinary Research. 2016; 47(1). https://doi.org/10.1186/s13567-016-0319-z PMID: 26911997

60. Park SL, Huang YJS, Lyons AC, Ayers VB, Hettenbach SM, McVey DS, et al. North American domestic pigs are susceptible to experimental infection with Japanese encephalitis virus. Sci Rep. 2018 May 21; 8(1):7951. https://doi.org/10.1038/s41598-018-26208-8 PMID: 29784969

61. Nemeth N, Bosco-Lauth A, Oesterle P, Kohler D, Bowen R. North American Birds as Potential Amplifying Hosts of Japanese Encephalitis Virus. Am J Trop Med Hyg. 2012 Oct 3; 87(4):760-7. https://doi.org/10.4269/ajtmh.2012.12-0141 PMID: 22927494
62. Ilkal MA, Dhanda V, Rao BU, George S, Mishra AC, Prasanna Y, et al. Absence of viraemia in cattle after experimental infection with Japanese encephalitis virus. Trans R Soc Trop Med Hyg. 1988; 82 (4):628–31. https://doi.org/10.1016/0035-9203(88)90536-6 PMID: 2855679

63. Shimoda H, Tamaru S, Morimoto M, Hayashi T, Shimojima M, Maeda K. Experimental infection of Japanese encephalitis virus in dogs. J Vet Med Sci. 2011 Sep; 73(9):1241–2. https://doi.org/10.1292/jvms.11-0142 PMID: 21566398

64. Edelman R, Schneider RJ, Vejjajiva A, Pornpibul R, Voodhikul P. Persistence of Virus-Specific IgM and Clinical Recovery after Japanese Encephalitis. Am J Trop Med Hyg. 1976 Sep 1; 25(5):733–8. https://doi.org/10.4269/ajtmh.1976.25.733 PMID: 183556

65. Prince HE, Tobler LH, Yeh C, Gefter N, Custer B, Busch MP. Persistence of West Nile virus-specific antibodies in viremic blood donors. Clin Vaccine Immunol. 2007 Sep; 14(9):1228–30. https://doi.org/10.1128/CVI.00233-07 PMID: 17652525

66. Subra R. Biology and control of Culex pipiens quinquefasciatus Say, 1823 (Diptera, Culicidae) with special reference to Africa. Int J Trop Insect Sci. 1981 Dec; 1(04):319–38.

67. Birley MH, Rajagopalan PK. Estimation of the survival and biting rates of Culex quinquefasciatus (Diptera: Culicidae). J Med Entomol. 1981 Jun; 18(3):181–6. https://doi.org/10.1093/jmedent/18.3.181 PMID: 7328603

68. Madder DJ, Surgeoner GA, Helson BV. Number of generations, egg production, and developmental time of Culex pipiens and Culex restauns (Diptera: Culicidae) in southern Ontario. J Med Entomol. 1983 May 26; 20(3):275–87. https://doi.org/10.1093/jmedent/20.3.275 PMID: 6876091

69. Chapman GE, Sherlock K, Hesson JC, Blagrove MSC, Lycett GJ, Archer D, et al. Laboratory transmission potential of British mosquitoes for equine arboviruses. Parasit Vectors. 2020 Aug 12; 13(1):413. https://doi.org/10.1186/s13071-020-04285-x PMID: 32787904

70. van den Hurk AF, Nisbet DJ, Hall RA, Kay BH, MacKenzie JS, Ritchie SA. Vector competence of Australian mosquitoes (Diptera: Culicidae) for Japanese encephalitis virus. J Med Entomol. 2003 Jan; 40 (1):82–90. https://doi.org/10.1603/0022-2585-40.1.82 PMID: 12597658

71. Gajanana A, Rajendran R, Samuel PP, Thenmozhi V, Tsai TF, Kimura-Kuroda J, et al. Japanese encephalitis in south Arcot district, Tamil Nadu, India: a three-year longitudinal study of vector abundance and infection frequency. J Med Entomol. 1997 Nov; 34(6):651–9. https://doi.org/10.1093/jmedent/34.6.651 PMID: 9439119

72. Gould DJ, Barnett HC, Suyemoto W. Transmission of Japanese encephalitis virus by Culex gelidus theobald. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1962; 56(5). https://doi.org/10.1016/0035-9203(62)90018-4 PMID: 13949551

73. The R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing; 2020.

74. Viennet E, Garros C, Gardès L, Rakotoarivony I, Allène X, Lancelot R, et al. Host preferences of Palaearctic Culicoides biting midges: implications for transmission of orbiviruses. Medical and Veterinary Entomology. 2013; 27(3):255–66. https://doi.org/10.1111/j.1365-2915.2012.01042.x PMID: 22985009

75. Wang J, Hihara E. A unified formula for calculating body surface area of humans and animals. Eur J Appl Physiol. 2004 Jun; 92(1–2):13–7. https://doi.org/10.1007/s00421-004-1074-9 PMID: 15015002

76. Seamans TW, Hamershock DW, Bernhardt GE. Determination of body density for twelve bird species. Ibis. 1995; 137(3):424–8.

77. Haycock GB, Schwartz GJ, Wisotsky DH. Geometric method for measuring body surface area: a height-weight formula validated in infants, children, and adults. J Pediatr. 1978 Jul; 93(1):62–6. https://doi.org/10.1016/s0022-3476(78)80601-5 PMID: 650346

78. Morris MD. Factorial Sampling Plans for Preliminary Computational Experiments. null. 1991 May 1; 33 (2):161–74.

79. Chheng K, Carter MJ, Emary K, Chanpheaktra N, Moore CE, Stoesser N, et al. A prospective study of the causes of febrile illness requiring hospitalization in children in Cambodia. PLoS ONE. 2013; 8(4): e60634. https://doi.org/10.1371/journal.pone.0060634 PMID: 23593267

80. Turner P, Suy K, Tan LV, Sar P, Miliya T, Hong NTT, et al. The aetiologies of central nervous system infections in hospitalised Cambodian children. BMC Infect Dis. 2017 29; 17(1):806. https://doi.org/10.1186/s12879-017-2915-6 PMID: 29284418

81. Tarantola A, Goutard F, Newton P, Lamballerie X de, Lortholary O, Cappelle J, et al. Estimating the Burden of Japanese Encephalitis Virus and Other Encephalitides in Countries of the Mekong Region. PLoS Negl Trop Dis. 2014 Jan 30; 8(1):e2533. https://doi.org/10.1371/journal.pntd.0002533 PMID: 24498443
82. Mao X, Zhou H. The spatiotemporal distribution of Japanese Encephalitis cases in Yunnan Province, China, from 2007 to 2017. PLoS ONE. 2020 Apr 14; 15(4):e0231661. https://doi.org/10.1371/journal.pone.0231661 PMID: 32287313

83. Tian HY, Bi P, Cazelles B, Zhou S, Huang SQ, Yang J, et al. How environmental conditions impact mosquito ecology and Japanese encephalitis: An eco-epidemiological approach. Environment International. 2015; 79:17-24. https://doi.org/10.1016/j.envint.2015.03.002 PMID: 25771078

84. Montini Maluda MC, Jelip J, Ibrahim MY, Suleiman M, Jeffree MS, Binti Aziz AF, et al. Nineteen Years of Japanese Encephalitis Surveillance in Sabah, Malaysian Borneo. Am J Trop Med Hyg. 2020 Aug; 103(2):864-868. https://doi.org/10.4269/ajtmh.19-0928 Epub 2020 Jun 4. PMID: 32524958; PMCID: PMC7410460.

85. Wangchuk S, Tamang TD, Darnal JB, Pelden S, Lhazeen K, Mynak ML, et al. Japanese Encephalitis Virus as Cause of Acute Encephalitis, Bhutan. Emerg Infect Dis. 2020 Sep; 26(9):2239-42. https://doi.org/10.3201/eid2609.200620 PMID: 32818416

86. Mighell E, Ward MP. African Swine Fever spread across Asia, 2018-2019. Transbound Emerg Dis. 2021 Sep; 68(5):2722-2732. https://doi.org/10.1111/tbed.14039 Epub 2021 Apr 27. PMID: 33599077.

87. Woonwong Y, Do Tien D, Thanawongnuwech R. The Future of the Pig Industry After the Introduction of African Swine Fever into Asia. Animal Frontiers. 2020 Oct 30; 10(4):30-7. https://doi.org/10.1093/af/vfaa037 PMID: 33150009

88. Gingrich JB, Nisalak A, Latendresse JR, Pomsdhit J, Paisansilp S, Hoke CH, et al. A longitudinal study of Japanese encephalitis in suburban Bangkok, Thailand. The Southeast Asian journal of tropical medicine and public health. 1987; 18(4):558-66. PMID: 2896391

89. Durand B, Haskouri H, Lowenski S, Vachieri N, Beck C, Lecollinet S. Seroprevalence of West Nile and Usutu viruses in military working horses and dogs, Morocco, 2012: dog as an alternative WNV sentinel species? Epidemiology & Infection. 2016 Jul; 144(9):1857-64. https://doi.org/10.1017/S095026881600011X PMID: 26838515

90. Pham-Thanh L, Nguyen-Tien T, Magnusson U, Bui-Nghia V, Bui-Ngoc A, Le-Thanh D, et al. Dogs as Sentinels for Flavivirus Exposure in Urban, Peri-Urban and Rural Hanoi, Vietnam. Viruses. 2021 Mar 19; 13(3). https://doi.org/10.3390/v13030507 PMID: 33808552