Spatial epidemiology of Japanese encephalitis virus and other infections of the central nervous system in Lao PDR (2003 – 2011): a retrospective analysis =========================================================================================================================================================== * Sayaphet Rattanavong * Audrey Dubot-Pérès * Mayfong Mayxay * Manivanh Vongsouvath * Sue J Lee * Julien Cappelle * Paul N. Newton * Daniel M. Parker ## ABSTRACT **Background** Central nervous system (CNS) infections are important contributors to morbidity and mortality and the causative agents for ∼50% patients are never identified. The causative agents of some CNS infections have distinct spatial and temporal patterns. **Methodology/Principal Findings** Here we present the results of a spatial epidemiological and ecological analysis of CNS infections in Lao PDR (2003 – 2011). The data came from hospitalizations for suspected CNS infection at Mahosot Hospital in Vientiane. Out of 1,065 patients, 450 were assigned a confirmed diagnosis. While many communities in Lao PDR are in rural and remote locations, most patients in these data came from villages along major roads. Japanese encephalitis virus ((JEV); n=94) and *Cryptococcus* spp. (n=70) were the most common infections. JEV infections peaked in the rainy season and JEV patients came from villages with higher surface flooding during the same month as admission. JEV infections were spatially dispersed throughout rural areas and were most common in children. *Cryptococcus* spp. infections clustered near Vientiane (an urban area) and among adults. **Conclusions/Significance** The spatial and temporal patterns identified in this analysis are related to complex environmental, social, and geographic factors. For example, JEV infected patients came from locations with environmental conditions (surface water) that are suitable to support larger mosquito vector populations. Most patients in these data came from villages that are near major roads; likely the result of geographic and financial access to healthcare and also indicating that CNS diseases are underestimated in the region (especially from more remote areas). As Lao PDR is undergoing major developmental and environmental changes, the space-time distributions of the causative agents of CNS infection will also likely change. There is a major need for increased diagnostic abilities; increased access to healthcare, especially for rural populations; and for increased surveillance throughout the nation. **AUTHOR SUMMARY** Infections of the central nervous system (CNS) are important with regard to public health. However many CNS infections are never diagnosed. In this analysis we investigated spatial and temporal patterns in hospitalized patients with suspected CNS infections in Lao PDR. We found that patients were most likely to come from villages located along major roads and highways. Patients from remote areas may have more difficulty reaching healthcare facilities. The most commonly diagnosed infection in these patients was Japanese encephalitis virus (JEV). Patients with this disease came from locations that were optimal for the mosquito vectors that spread JEV, rural areas with surface water and during the rainy season. Our analyses suggest that CNS infections should be a priority for public health workers in the region. Diagnostic capabilities should be increased throughout the nation; surveillance efforts should be broadened; and efforts should be increased toward providing easy access to healthcare for rural and remote populations. ## INTRODUCTION Numerous illnesses go undiagnosed and the causative agents of many infections are never identified. In regions where access to healthcare facilities is limited and where diagnostic capabilities are few, a smaller proportion of diseases are objectively diagnosed. Diseases with mild symptoms may more frequently go untreated, but in some areas even severe illnesses commonly go undiagnosed and untreated. Infections of the central nervous system (CNS) can be particularly severe, affecting the brain and/or spinal cord and/or the surrounding meninges, frequently resulting in death. Pathogens that invade and infect the CNS include viruses, bacteria, fungi, parasites, and prions. As with other infections, the causative agent(s) of many CNS infections are never determined (frequently <50%) [1–3]. Symptoms of CNS infection can range widely, even for single causative agents, leading to further difficulties with diagnosis. Diagnoses are frequently presumptive and non-specific (i.e. meningitis is often presumed to be caused by bacteria whereas encephalitis is presumed to be caused viruses [2,4–6]). Some causative agents are specific to regions (e.g. Japanese Encephalitis, Saint Louis Encephalitis, Rift Valley Fever Viruses) and exhibit seasonal fluctuations (e.g. vector borne infections), therefore geography and seasonality can facilitate presumptive diagnosis of CNS diseases [4,7]. In Southeast (SE) Asia, CNS infections are increasingly recognized as important contributors to morbidity and mortality [1,8,9]. However, detailed medical and epidemiological data are frequently lacking, especially for low income nations and from remote areas within middle-to-high income nations. Known important viral CNS infections in SE Asia include Japanese encephalitis, dengue, and rabies viruses [1]. Important bacterial CNS infections include *Streptococcus pneumoniae, Haemophilus influenzae, S. suis, Mycobacterium tuberculosis* and *Neisseria meningitidis. Orientia tsutsugamushi, Rickettsia typhi* and *Leptospira* spp. are increasingly recognized as important causes [9]. Detailed analyses that confirm the cause of CNS related infections or assess their spatial and temporal distribution in the region are rare [8–11,14]. We recently published the results of a study of the etiology and impact of CNS infections diagnosed among 1,065 patients at Mahosot Hospital, Vientiane, Lao PDR [14]. The goal of this secondary analysis was to investigate the spatial distribution(s) of CNS-related infections; to look for differential spatial distributions for major causative agents; and to explore potential geographic, demographic, and environmental correlates of these infections. ## DATA AND METHODS ### Data sources, processing and merging Data used in this research were compiled from four main sources (**Supporting Figure 1**). The epidemiological data come from an 8-year research project on CNS infections in Lao PDR from all patients who received diagnostic lumbar puncture (LP) at Mahosot Hospital in Vientiane, Lao PDR between January 2003 and August 2011 and consenting to participate [14]. All patients were admitted to the hospital because of suspected CNS infection and Mahosot Hospital is the only medical facility in Lao PDR capable of performing diagnostic LP and cerebral spinal fluid (CSF) analysis. Tests for a large panel of pathogens were performed at the Microbiology Laboratory following international standards (details provided in **Supporting Materials I** and [14]). Demographic (age, gender, ethnicity) and geographic (home village) characteristics of patients were recorded in the database. The epidemiological data were used to create two separate datasets: One aggregated at the village level (one row per location) and another was maintained at the individual level, with one row per individual. The official Lao PDR censuses from 2005 and 2015 were used to geocode villages (based on village name and administrative units listed in patient records) and to assign village population estimates to each village (taking a mean population size between 2005 and 2015). Village location and population sizes were then merged to both the individual- and village-level datasets. A subset of villages within the geographic region of the home villages of included patients was selected by overlaying a standard deviational ellipse with 3 standard deviations (calculations described in **Supporting Materials II**) around the patient home villages and then selecting all villages within that ellipse (**Supporting Figure 2**). These villages were then retained for village level comparisons between villages populated, and not populated, with patients admitted with CNS disease needing an LP. This subset of villages is hereafter referred to as the “study area”. Major road network data was taken from OpenStreetMaps ([http://www.openstreetmap.la](http://www.openstreetmap.la)), selecting “primary”, “secondary”, and all major connecting roads (downloaded in February 2017, **Supporting Figure 3**). Primary and secondary roads are the two largest road classifications for the nation. Primary roads link major towns and cities and secondary roads link mid-sized towns. Primary and secondary link roads are ramps or slip roads that connect other roads to primary or secondary roads. Together, these types of roads are hereafter referred to as “major roads”. Smaller roads were not included in this analysis as they are less likely to be accurately included in the OpenStreetMaps data. The Euclidian distance was then calculated from all villages in the census to the nearest point along a major road. These distances were merged to both the village- and the individual-level datasets. Environmental predictor variables for vegetation and surface water were derived from Moderate Resolution Imaging Spectroradiometer (MODIS) products (MOD13Q1/MYD13Q1 250 meter AQUA/TERRA 16 day composites). Since many infectious diseases, especially vector borne diseases, are strongly influenced by environmental factors we hypothesized that indicators of vegetation and surface flooding would correlate with some specific diagnoses. Three environmental indices were downloaded and considered in these analyses: a normalized flooding index (NFI) [15]; the normalized difference vegetative index (NDVI); and the enhanced vegetation index (EVI). NFI is indicative of surface water, NDVI is indicative of green surface vegetation, and EVI is an improved measure of green vegetation that is intended to account for dense forest canopies and atmospheric conditions that can lead to error in NDVI measurements. Data were downloaded for each of these environmental indices (EI) within each 16-day time period from February 2002 through December of 2011. The final analyses conducted in this research retained only the EVI and NFI for environmental predictors (summary statistics for NDVI are included) because NDVI and EVI were strongly correlated. The EIs are described in more detail in **Supporting Materials III**. The environmental raster data were then summarized and extracted based on varying buffer sizes (2km, 5km, and 10km) for each village in the individual- and village-level datasets. Permanent water bodies (such as the Mekong River and Nam Ngum Dam) were masked from the NFI calculations. For the village-level datasets, mean values of each environmental variable was calculated for the study period duration and used as an indicator of “average” vegetation or surface water characteristics of each village. For the individual-level dataset the values were extracted based on the admission date, using increasing durations of time prior to admission (within the same month, within the previous 2 months, within the previous year). The final datasets include the village-level data, that is a subsample of 98% of all villages with patients included in the study and the other census designated villages within the same region (the study area), and the individual level dataset that includes all patients included in the study. Variables used in this analysis and their descriptions are listed in **Table 1**. View this table: [Table 1:](http://medrxiv.org/content/early/2019/10/25/19005884/T1) Table 1: List of variables, their spatial and temporal scales, and transformations ### Exploratory spatial data analysis Summary statistics (median; Q1:Q3; mean) were calculated for the distances between villages and the nearest major road, and aggregated by whether or not the village was home to included patients and by specific diagnoses. Summary statistics (mean and 95% confidence intervals) were also calculated for all environmental variables, at each buffer size and temporal resolution, and for each of the major diagnoses. Tukey’s post hoc range test was used to assess statistically significant differences in environmental indices across diagnoses. Standard distance deviations (SDDs) and standard deviation ellipses (SDEs) were calculated (**details in Supporting Materials II**) and mapped to measure and visually analyze the central tendency and spatial distributions for all patient home villages and by each of the major single (mono-infection) diagnoses. ### Formal analyses Multivariable regressions were used to calculate model-adjusted odds ratios and confidence intervals. The regressions at the village level focused on study patient villages and the home villages of JEV diagnosed patients. A multivariable regression was also done at the individual level focusing on JEV infected patients. Other diagnoses were not included in these analyses because of small numbers of cases per village. Logistic generalized additive models (GAMs) were used for variable selection and specification (**detailed in Supporting Materials IV and in Table 1**) for the final models. The GAMs were used to examine different specifications of the continuous environmental, geographic, and demographic variables and for changes in model fit and strength of association across buffer sizes (i.e. 2km, 5km, or 10km buffers) and for different time durations for EI measurements (i.e. same month, 2 months prior, 12 months prior to hospital admission). The final model covariates were chosen based on a combination of *a priori* hypotheses, model fit (using the Akaike information criterion), and strength of association between the covariate and the model outcome variable. The final model for the individual-level analysis was a generalized logistic mixed model with a random effect for home village. The final model for the village-level analysis was a logistic regression. ### Software All maps were created using QGIS version 3.4.9. R cran version 3.5.2 was used for downloading, processing, and wrangling MODIS data (using the “MODIStsp”; “raster”; “rgdal”; and “maptools” packages) and for all regressions. The “mgcv” package was used for GAMs and the lme4 package was used for the generalized mixed models. Euclidian distances between villages and major roads were calculated using QGIS. ArcMap version 10.5.1 was used to calculate SDDs and SDEs. ### Ethics approval Ethical clearance for the CNS study was granted by the Oxford University Tropical Ethics Research Committee and by the Ethical Review Committee of the Faculty of Medical Sciences, National University of Laos. Between 2003 and 2006, both Oxford Tropical Ethics Committee and the Ethical Review Committee of the Faculty of Medical Sciences, National University of Laos, approved the use of oral witnessed consent. Oral consent was confirmed by a signed witness statement. Verbal consent (from 2003 – 2006) and written consent (from 2006 – 2011) were obtained from all recruited patients or immediate relatives. All ages were included provided that they or their guardian gave witnesses oral informed consent (2003 – 2006) and written informed consent since 2006. ## RESULTS ### Summary statistics A total of 1,065 patients were recruited with no LP contraindications and consented to have a diagnostic LP; 450 (42%) were assigned a final laboratory diagnosis. The most common of these were Japanese encephalitis virus ((JEV) 94 individuals); followed by *Cryptococcus* spp. with 70 individuals (9 were *C. gattii*)); scrub typhus (*Orientia tsutsugamushi*) 31; Dengue virus 27; *Leptospira* spp. 25; murine typhus (*Rickettsia typhi*) 24; *Streptococcus pneumoniae* in 22 and 20 with *Mycobacterium tuberculosis*. 124 patients died prior to discharge (out of 893 with recorded discharge type recorded). The majority (666, 63%) of patients were male, with the lowest sex bias in cryptococcal infections (40/70, 57% male) and the highest among dengue infections (22/27, 82% male) (**Table 2**). Age patterns were evident in JEV and cryptococcal infections, with median ages of 13 and 33 years, respectively (**Table 2**). Patients were linked to 582 different villages (multiple patients could come from the same village): 90 villages with JEV patients, 66 with cryptococcal patients, 31 with scrub typhus patients, 27 with dengue patients, 24 with leptospiral patients, and 24 with murine typhus patients. The majority (870, 82%) of patients came from within Vientiane Prefecture (678, 64%) or Vientiane Province (192, 18%). View this table: [Table 2:](http://medrxiv.org/content/early/2019/10/25/19005884/T2) Table 2: Age and gender of study patients. (Q1 and Q3 indicate the first and third quartiles, respectively). A total of 6,416 villages (of 10,520 recorded in 2005 [16]) were selected as the study area for further village level analyses (**Table 3**). Villages that were home to study patients were 11 times (0.7km versus 6.3km, from **Table 3**) closer to a major road when compared to other villages within the study area. Scrub typhus and JEV infected patient homes were further from major roads than other patients, but the difference was not statistically significant in univariate analyses. View this table: [Table 3:](http://medrxiv.org/content/early/2019/10/25/19005884/T3) Table 3: Distribution of distances (in km) to the nearest major road, by diagnosis type. Counts of villages are from within 3 standard deviational ellipses (SDEs) of all LP villages (referred to as the “study area” in text). In some cases, multiple patients came from the same village meaning that counts of villages will be smaller than counts of total patients. (Q1 and Q3 indicate the first and third quartiles, respectively). The home villages of JEV patients were more broadly dispersed (**Figure 1B**) than for patients with other etiologies (**Figure 1C**), evident from the circular, larger SDE and SDD. The distribution of these JEV patient home villages was also relatively isotropic, with the SDE and SDD being nearly identical. Conversely, patients with cryptococcal infections were clustered near Vientiane City and along the road leading North/Northwest from the urban center (**Figure 1C**). Scrub typhus and murine typhus infections were also both clustered around Vientiane City but showed perpendicular spatial distributions (**Supporting Figures 4D and 4G**) a pattern previously described from IgG seropositivity data from Vientiane City [17]). ### Characteristics of patient home villages Mean NFI values for the 2km radius tended to be higher than for either the 5km or 10km radius as surface flooding is heterogeneous and taking a mean across larger radii dilutes the measurement. The opposite pattern was observed for both vegetation indices. The 2km radius for both mean NDVI and mean EVI was usually smaller than at 5km or 10km radii (**Figures 2 and 3**). Study patient villages had higher mean NFI values than non-study patient villages (**Figure 2A**) (non-study patient villages are those in the same study area as study patient villages but were not home to a study patient). The home villages of study patients diagnosed with *dengue virus* and cryptococcal infections had high mean NFI values over the duration of the study period when compared to other major diagnoses (**Figure 2A**). Conversely, the home villages of study patients tended to have lower mean EVI values than non-study patient villages (**Figure 2C**). Villages from which patients who were diagnosed with JEV were an exception to this general pattern. JEV patient home villages had higher mean EVI when compared to home villages of patients with dengue and cryptococcal infections (**Figure 2C**). Home villages of patients with JEV diagnoses had lower mean NFI over the duration of the study period, but had higher NFI than other major diagnoses when looking specifically at the month of admission (especially when compared to cryptococcal infections and murine typhus (**Figure 3A**)). JEV patient home villages also had higher EVI during the month of admission than most other mono-infections, especially when compared to either cryptococcal infections or murine typhus (**Figure 3C**). Scrub typhus infections had higher EVI than cryptococcal infections when the measurement was taken at the 10km radius buffer (not detectable at smaller radii (**Figure 3C**)). Home villages of patients with dengue infections had higher NFI than murine typhus or *Cryptococcus* spp. patient home villages when considering the 2 months prior to admission (**Figure 3D**). *Cryptococcus* spp. patient home villages had particularly low EVI in the two months leading up to admission, especially when compared to JEV and *Leptospira* spp. patient home villages (**Figure 3F**). At one year prior to admission both *Cryptococcus* spp. patient home and *dengue virus* patient home villages had higher NFI than JEV patient home villages (**Figure 3G**). ### Logistic regressions for geographic, environmental and demographic predictors of CNS infections In agreement with univariate analyses, villages from which study patients came tended to be larger in population size (**Supporting Figure 4**), lower in elevation (**Supporting Figure 5**), and closer to a major road when compared to other villages within the study area (**Table 4**). They also had higher mean levels of surface flooding, with villages in the highest NFI quadrant having over two times the odds (AOR: 2.21; CI: 1.49 – 3.31) of being a home village for study patients when compared to neighboring villages in the study area (**Table 4 and Figure 2**). Villages from which JEV patients originated had few defining characteristics in the logistic regression, other than being larger in population size (AOR: 1.74; CI: 1.55 – 1.96) and at lower elevations (AOR: 0.69; CI: 0.46 – 0.97) than non-study patient villages (**Table 5**). View this table: [Table 4:](http://medrxiv.org/content/early/2019/10/25/19005884/T4) Table 4: Logistic regression and model adjusted odds ratios (AOR) for village level analysis of LP villages View this table: [Table 5:](http://medrxiv.org/content/early/2019/10/25/19005884/T5) Table 5: Logistic regression and model adjusted odds ratios (AOR) for village level analysis of JEV villages In the individual-level analysis (**Table 6**), age and season were the strongest predictors of JEV infections among all patients. Patients who were admitted between July and September had over seven times the odds (AOR: 7.40; CI: 1.45 – 37.67) of being diagnosed with JEV when compared to patients who were admitted between January and March (**Table 6**). JEV was most common in children aged 5 through 14 (AOR: 2.74; CI: 1.31 – 5.69; ages 0 – 4 as the comparison group). NFI during the month of admission (10km buffer used in the regression) was also a strong predictor of JEV infection. Individuals who came from villages in the highest quadrant of NFI measurements had approximately 3 times the odds being diagnosed with JEV (AOR: 3.06; CI: 1.04 – 8.96). While study patients came from villages with lower mean elevations, patients who were diagnosed with JEV came from higher elevation villages in comparison to the other patients (AOR: 1.36; CI: 1.11 – 1.66). EVI was a significant predictor in models that did not include distance to road, village population, and elevation (**Table 6 M1 and M2**). View this table: [Table 6:](http://medrxiv.org/content/early/2019/10/25/19005884/T6) Table 6: Mixed effects logistic regression and model adjusted odds ratios (AOR) for individual level analysis ## DISCUSSION Patients from this study were recruited based on symptomology and a medical procedure that is only available at a single location in the nation (diagnostic LP at Mahosot Hospital, Vientiane). The home villages of all included study patients, regardless of diagnosis, were approximately centered on Vientiane City and were closer to major roads than would be expected by chance alone. For many of the infections studied in this analysis, this association is likely more related to geographic and financial access to healthcare systems rather than exposure to environmental risk factors – especially for infections that are more associated with rural areas (e.g. JEV). This finding also suggests that the results here may be a vast underestimate of the true burden of CNS infections, with much of the Lao population not being in near proximity to a major road (**Table 3**). The causative agents of CNS infections differ in biology, ecology, and geography, and this is evident through the spatial distributions of the home villages of patients. The geographic, environmental, and demographic patterns exhibited by patients needing a diagnostic LP for suspected CNS infections, and for specific diagnoses, are the result of complex overlapping factors. A similar spatial pattern was described from an epidemiological analysis of CNS infections among children admitted to Ho Chi Minh City hospitals in Vietnam – with most patients coming from districts near the hospital [18]. While the majority of infections (55%) in the Vietnam study were presumed to be bacterial in origin, in this study from Lao PDR bacterial infections were identified in only 38% (170 out of 450 patients with diagnoses). JEV was the single largest identified cause of CNS infections in these data; it primarily affected children (median 13 years of age, **Table 2**), occurred predominantly during the rainy season (likely corresponding to peaks in mosquito vector populations), and in villages with recent high levels of surface water [19]. JEV is a vaccine-preventable disease, but the vaccine has historically been expensive and vaccine programs are frequently limited by access to remote communities. In 2013 the WHO approved a less expensive vaccine (produced in China by the Chengdu Institute of Biological Products) which has since been used in mass vaccination campaigns in Lao PDR in 2013 and 2015 [20]. The vaccine is now routinely given to all children less than 9 months of age but coverage may be low in some areas. The second largest contributor to CNS diseases were cryptococcal infections, which are opportunistic fungal infections with high mortality [21]. Of the 70 patients with cryptococcal infections, 12 died prior to discharge and another 8 likely died at home after leaving the hospital. Cryptococcal infections are generally acquired after inhalation of the yeast-like form of the fungus which has been associated with several ecological habitats (*Cryptococcus gattii* has been associated with over 50 species of trees; *Cryptococcus neoformans* has been associated with bird droppings but is also suspected to be associated with plants [22,23]). This disease has a long incubation period [24,25] and while *C. gattii* infections commonly occur among immunocompetent individuals, *C. neoformans* infections are frequently associated with HIV infections [26,27]. In these data, 20 of the patients diagnosed with cryptococcal infections also had HIV infections. While several studies have shown that cryptococcal species exist in specific ecological habitats and have inferred environmental exposure, the long incubation period and complex natural history likely obfuscate ecological correlations. There are several limitations to this research. Diagnostic LPs are only conducted in Mahosot Hospital in the national capital. Logistical and financial difficulties in accessing healthcare facilities, and especially for etiological diagnostic capabilities, likely leads to severe under-reporting of meningitis, encephalitis, or in the diagnosing the causative agent in these conditions when the patient does access care. All of these factors ultimately lead to small case counts for numerous different causative agents. The spatial patterns in points (villages) and ellipses exhibited in these data are likely influenced by the shape of the nation and it is possible that the point patterns and ellipses would differ if we had data from neighboring nations. Spatial and temporal patterns that differentiate different infections might be more obvious if the surveillance system instead focused on any symptomatic infections (rather than only suspected infections of the CNS). Some pathogens are neurotropic whereas others have tropism for other organs, while being capable of occasionally infecting the CNS. This may partially explain the higher case counts of JEV and why we were able to identify spatial, temporal, and environmental predictors for this causative agent. OpenStreetMaps data are volunteered data and may be prone to error. For this reason we focused on major roads, whose routes have changed very little over the last decades. For the regressions, the distances from all villages to the nearest major road was also rounded to the nearest 5km. Examination of satellite imagery in comparison with the major roads from OpenStreetMaps suggests that where error does exist, it is on a scale of +100 meters, meaning that measurements of distances, as used in this analysis, should not be strongly influenced. Some of these data now come from over a decade ago. Surveillance systems of this type (based on relatively vague symptomology), with a wide panel of possible contributing causative agents, and necessary intensive laboratory components are extremely labor and time intensive. Lao PDR is currently undergoing vast environmental, demographic, and economic changes. Road networks are increasing in range and density and several areas (i.e. Vientiane, Savannakhet) are undergoing expansive urbanization [28]. These environmental changes will most likely result in shifting patterns of infectious diseases. As the region undergoes urbanization (including both a decrease in urban landscape and movement of human populations to urban centers), pathogens that thrive in rural areas (e.g. JEV) may undergo reduced transmission, especially if vaccine campaigns are more capable of reaching rural populations. Conversely, infections that cluster in urban and peri-urban areas (such as dengue and murine typhus) may increase in frequency. Several environmental indices from remote sensing instruments have shown potential for predicting disease risk, differentiating disease types, or for other surveillance efforts in SE Asia [29,30] and globally [31–35]. This analysis, and others like it, illustrates the ability to differentiate some infections (namely JEV when compared to other diagnoses) through the use of freely available data (i.e. MODIS) and software (R and QGIS) and routinely collected healthcare data. Surveillance systems and potentially diagnostic algorithms [36] in developing settings could benefit from inclusion of such resources. A far-reaching surveillance system that is representative of the entire nation and includes likely CNS infections would be beneficial in order to assess the true burden of CNS infections – many of which would benefit from primary and secondary prevention through increased provision of vaccines, vector control, and early diagnosis and treatment. Given the inherent difficulties in accurately diagnosing and treating CNS infections, the predictors reported here and from other epidemiological studies for major contributors to CNS diseases (i.e. age, seasonality, location, and environmental characteristics) could be considered alongside clinical symptomology when presumptive diagnoses are being made. However, it will be important to consider current and ongoing demographic, environmental, and economic changes in Lao PDR. Finally, increasing population access to vaccines, diagnosis, and treatment would have clear benefits to overall population health. As with other parts of the developing world, a large fraction of the Lao population must travel long distances in order to reach primary healthcare centers. In 2005 73% of the Lao population was reported to live in rural areas, 21% without roads. By 2015 67% of the population were reported to live in rural villages with 8% in villages without roads [37]. For many communities, travel during the wet season remains difficult. Travel costs can also be prohibitive. Most of the CNS infections in this analysis occurred or developed symptoms during the wet season. Public health initiatives that help to decrease the distances between communities and the healthcare services that they need are warranted. ## Data Availability Data from this manuscript come from four different sources. 1.) Lao Bureau of Statistics Census data, available through: http://www.decide.la/en/ 2.) Remote sensing data available through: https://modis.gsfc.nasa.gov/data/ 3.) Roads and highway data available through: http://www.openstreetmap.la 4.) Patient hospital records which can be requested from the Mahidol Oxford Tropical Medicine Research Unit (following policies and guidelines): http://www.tropmedres.ac/data-sharing [http://www.decide.la/en/](http://www.decide.la/en/) [https://modis.gsfc.nasa.gov/data/](https://modis.gsfc.nasa.gov/data/) [http://www.openstreetmap.la](http://www.openstreetmap.la) [http://www.tropmedres.ac/data-sharing](http://www.tropmedres.ac/data-sharing) ## Funding This study was supported by the European Commission Innovate program (ComAcross project, grant no. DCI-ASIE/2013/315-047). The work of LOMWRU, PN and MM are funded by the Wellcome Trust of Great Britain, grant 106698/Z/14/Z. Travel to Lao PDR to support DMP for this research was obtained through a University of California Council on Research, Computing, and Libraries (CORCL) Faculty Grant. ## Supporting information legends **Supporting information file 1: Detailed explanation of variables and analyses** ## FIGURES Figure 1: Spatial distributions of the home villages of study patients, for A: all study patients, B: study patients with JEV infections, and C: with cryptococcal infections, D: scrub typhus infections, E: with *dengue virus* infections, F: with leptospiral infections, and G: with murine typhus infections. SDDs and SDEs are weighted by case numbers, with some patients coming from the same village. Maps were created using QGIS version 3.4.9. All layers were created by the authors of this manuscript. Figure 2: Environmental indices for villages with study patient homes for the duration of the study period (January 2003 through August 2011) for all study patient villages, non study patient villages in the study area, and for major diagnoses (JEV = Japanese Encephalitis virus; Crypto = cryptococcal infection; ST = scrub typhus; MT = murine typhus; dengue = Dengue virus; lepto = *Leptospira* spp. infection). Bar values are mean values and the error bars are 95% confidence intervals, using the t-distribution. NFI values here have a constant (0.25) added to them for visualization only. Figure 3: Environmental indices for study patients by major diagnosis and at different times leading up to the date of admission. JEV = Japanese Encephalitis virus; Crypto = cryptococcal infection; ST = scrub typhus; MT = murine typhus; dengue = *Dengue virus*; lepto = *Leptospira* spp. infection. Bar values are mean values and the error bars are 95% confidence intervals, using the t-distribution. NFI values here have a constant (0.25) added to them for visualization only. ## Acknowledgements We are very grateful to the patients and to Bounthaphany Bounxouei, past Director of Mahosot Hospital; to Bounnack Saysanasongkham, Director of Department of Health Care, Ministry of Health; to H.E. Bounkong Syhavong, Minister of Health, Lao PDR; to the staff of the Infectious Disease Center and Microbiology Laboratory for their help and support. ## Footnotes * Sayaphet Rattanavong: Sayaphet{at}tropmedres.ac, Audrey Dubot-Pérès: Audrey{at}tropmedres.ac, Mayfong Mayxay: Mayfong{at}tropmedres.ac, Manivanh Vongsouvath: Manivanh{at}tropmedres.ac, Sue J Lee: Sue{at}tropmedres.ac, Julien Cappelle: Julien.cappelle{at}cirad.fr, Paul N. Newton: Paul.Newton{at}tropmedres.ac * Received September 2, 2019. * Revision received October 15, 2019. * Accepted October 25, 2019. * © 2019, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## REFERENCES 1. 1.Tarantola A, Goutard F, Newton P, de Lamballerie X, Lortholary O, Cappelle J, et al. Estimating the Burden of Japanese Encephalitis Virus and Other Encephalitides in Countries of the Mekong Region. PLoS Neglected Tropical Diseases. 2014;8: 4. doi:10.1371/journal.pntd.0002533 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0002533&link_type=DOI) 2. 2.Glaser C a, Honarmand S, Anderson LJ, Schnurr DP, Forghani B, Cossen CK, et al. Beyond viruses: clinical profiles and etiologies associated with encephalitis. Clinical infectious diseases!]: an official publication of the Infectious Diseases Society of America. 2006;43: 1565–1577. doi:10.1086/509330 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/509330&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17109290&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000242126300007&link_type=ISI) 3. 3.Glaser CA, Gilliam S, Schnurr D, Forghani B, Honarmand S, Khetsuriani N, et al. In Search of Encephalitis Etiologies: Diagnostic Challenges in the California Encephalitis Project, 1998–2000. Clinical Infectious Diseases. 2003;36: 731–742. doi:10.1086/367841 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/367841&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12627357&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000181451200008&link_type=ISI) 4. 4.Solomon T, Michael BD, Smith PE, Sanderson F, Davies NWS, Hart IJ, et al. Management of suspected viral encephalitis in adults – Association of British Neurologists and British Infection Association National Guidelines. Journal of Infection. 2012;64: 347–373. doi:10.1016/j.jinf.2011.11.014 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2011.11.014&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22120595&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000301927300001&link_type=ISI) 5. 5.McGill F, Heyderman RS, Michael BD, Defres S, Beeching NJ, Borrow R, et al. The UK joint specialist societies guideline on the diagnosis and management of acute meningitis and meningococcal sepsis in immunocompetent adults. Journal of Infection. 2016;72: 405–438. doi:10.1016/j.jinf.2016.01.007 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2016.01.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26845731&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) 6. 6.Bharucha T, Vickers S, Ming D, Lee SJ, Dubot-Pérès A, de Lamballerie X, et al. Association between reported aetiology of central nervous system infections and the speciality of study investigators—a bias compartmental syndrome? Trans R Soc Trop Med Hyg. 2017;111: 579–583. doi:10.1093/trstmh/try008 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/trstmh/try008&link_type=DOI) 7. 7.Kennedy PGE. Viral Encephalitis: Causes, Differential Diagnosis, and Management. Journal of Neurology, Neurosurgery & Psychiatry. 2004;75: i10–i15. doi:10.1136/jnnp.2003.034280 [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czoxNDoiNzUvc3VwcGxfMS9pMTAiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAxOS8xMC8yNS8xOTAwNTg4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 8. 8.Olsen SJ, Campbell AP, Supawat K, Liamsuwan S, Chotpitayasunondh T, Laptikulthum S, et al. Infectious Causes of Encephalitis and Meningoencephalitis. 2015;21: 2003–2005. 9. 9.Dittrich S, Rattanavong S, Lee SJ, Panyanivong P, Craig SB, Tulsiani SM, et al. Orientia, rickettsia, and leptospira pathogens as causes of CNS infections in Laos: A prospective study. The Lancet Global Health. 2015;3: e104–e112. doi:10.1016/S2214-109X(14)70289-X [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2214-109X(14)70289-X&link_type=DOI) 10. 10.Mai NTH, Phu NH, Nhu LNT, Hong NTT, Hanh NHH, Nguyet LA, et al. Central Nervous System Infection Diagnosis by Next-Generation Sequencing: A Glimpse Into the Future? Open Forum Infect Dis. 2017;4. doi:10.1093/ofid/ofx046 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ofid/ofx046&link_type=DOI) 11. 11.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;17. doi:10.1186/s12879-017-2915-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-017-2915-6&link_type=DOI) 12. 12.Ho NT, Hoang VMT, L. NNT, Nguyen DT, Tran A, Kaki D, et al. A spatial and temporal analysis of paediatric central nervous system infections from 2005 to 2015 in Ho Chi Minh City, Vietnam. Epidemiology and Infection. 2017;145: 3307–3317. doi:10.1017/S095026881700228X [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S095026881700228X&link_type=DOI) 13. 13.Mawuntu AHP, Bernadus JBB, Dhenni R, Wiyatno A, Anggreani R, Feliana, et al. Detection of central nervous system viral infections in adults in Manado, North Sulawesi, Indonesia. PLOS ONE. 2018;13: e0207440. doi:10.1371/journal.pone.0207440 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0207440&link_type=DOI) 14. 14.Dubot-Pérès A, Mayxay M, Phetsouvanh R, Lee SJ, Rattanavong S, Vongsouvath M, et al. Management of Central Nervous System Infections, Vientiane, Laos, 2003–2011 - Volume 25, Number 5—May 2019 - Emerging Infectious Diseases journal - CDC. doi:10.3201/eid2505.180914 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2505.180914&link_type=DOI) 15. 15.Boschetti M, Nutini F, Manfron G, Brivio PA, Nelson A. Comparative analysis of normalised difference spectral indices derived from MODIS for detecting surface water in flooded rice cropping systems. PLoS ONE. 2014;9. doi:10.1371/journal.pone.0088741 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0088741&link_type=DOI) 16. 16.Lao Population and Housing Census 2015: Provisional Report [Internet]. Ministry of Planning and Investment Lao Statistics Bureau; 2015 Dec p. 64. Available: [https://lao.unfpa.org/sites/default/files/pub-pdf/Final%20report-editting-English1.pdf](https://lao.unfpa.org/sites/default/files/pub-pdf/Final%20report-editting-English1.pdf) 17. 17.Vallé J, Thaojaikong T, Moore CE, Phetsouvanh R, Richards AL, Souris M, et al. Contrasting spatial distribution and risk factors for past infection with scrub typhus and murine typhus in Vientiane city, Lao PDR. PLoS Neglected Tropical Diseases. 2010;4: 1–10. doi:10.1371/journal.pntd.0000909 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0000909&link_type=DOI) 18. 18.Ho NT, Hoang VMT, L. NNT, Nguyen DT, Tran A, Kaki D, et al. A spatial and temporal analysis of paediatric central nervous system infections from 2005 to 2015 in Ho Chi Minh City, Vietnam. Epidemiology and Infection. 2017;145: 3307–3317. doi:10.1017/S095026881700228X [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S095026881700228X&link_type=DOI) 19. 19.Campbell G, Hills S, Fischer M, Jacobson J, Hoke C, Hombach J, et al. Estimated global incidence of Japanese encephalitis: Bulletin of the World Health Organization. 2011;89: 766–774. doi:10.2471/BLT.10.085233 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2471/BLT.10.085233&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22084515&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295707700015&link_type=ISI) 20. 20.UNICEF. Japanese Encephalitis Vaccine Market and Supply Update. 2015; 21. 21.Slavin MA, Chakrabarti A. Opportunistic fungal infections in the Asia-Pacific region. Medical Mycology. 2012;50: 18–25. doi:10.3109/13693786.2011.602989 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/13693786.2011.602989&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21905945&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) 22. 22.Springer DJ, Chaturvedi V. Projecting global occurrence of Cryptococcus gattii. Emerging Infectious Diseases. 2010;16: 14–20. doi:10.3201/eid1601.090369 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1601.090369&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20031037&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) 23. 23.Ellis DH, Pfeiffer TJ. Ecology, life cycle, and infectious propagule of Cryptococcus neoformans. The Lancet. 1990;336: 923–925. doi:10.1016/0140-6736(90)92283-N [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0140-6736(90)92283-N&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1976940&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990EC45600014&link_type=ISI) 24. 24.Castrodale LJ, Gerlach RF, Preziosi DE, Frederickson P, Lockhart SR. Prolonged incubation period for cryptococcus gattii infection in cat, Alaska, USA. Emerging Infectious Diseases. 2013;19: 1034–1035. doi:10.3201/eid1906.130006 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1906.130006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23735429&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) 25. 25.Georgi A, Schneemann M, Tintelnot K, Calligaris-Maibach RC, Meyer S, Weber R, et al. Cryptococcus gattii meningoencephalitis in an immunocompetent person 13 months after exposure. Infection. 2009;37: 370–373. doi:10.1007/s15010-008-8211-z [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s15010-008-8211-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19390780&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269254600014&link_type=ISI) 26. 26.Wright P, Inverarity D. Human immunodeficiency virus (HIV) related cryptococcol meningitis in rural Central Thailand-Treatment difficulties and prevention strategies. SOUTHEAST ASIAN J TROP MED PUBLIC HEALTH. 2007;38: 4. 27. 27.Inverarity D, Bradshaw Q, Wright P, Grant A. The spectrum of HIV-related disease in rural Central Thailand. SOUTHEAST ASIAN J TROP MED PUBLIC HEALTH. 2002;33: 10. 28. 28.1. Steinberg F, 2. Hakim J Bank AD. Urban development in the Greater Mekong Subregion. Steinberg F, Hakim J, editors. Mandaluyong City, Philippines: Asian Development Bank; 2016. doi:10.1111/j.1467-9787.2010.00709.x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1467-9787.2010.00709.x&link_type=DOI) 29. 29.1. Schumann GJ-P Ledien J, Sorn S, Hem S, Huy R, Buchy P, Tarantola A, et al. Assessing the performance of remotely-sensed flooding indicators and their potential contribution to early warning for leptospirosis in Cambodia. Schumann GJ-P, editor. PLOS ONE. 2017;12: e0181044. doi:10.1371/journal.pone.0181044 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0181044&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28704461&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) 30. 30.Chadsuthi S, Chalvet-Monfray K, Wiratsudakul A, Suwancharoen D, Cappelle J. A remotely sensed flooding indicator associated with cattle and buffalo leptospirosis cases in Thailand 2011–2013. BMC Infectious Diseases. 2018;18. doi:10.1186/s12879-018-3537-3 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-018-3537-3&link_type=DOI) 31. 31.Rogers DJ, Randolph SE, Snow RW, Hay SI. Satellite imagery in the study and forecast of malaria. Nature. 2002;415: 710–715. doi:10.1038/415710a [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/415710a&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11832960&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173709100061&link_type=ISI) 32. 32.Lambin EF, Tran A, Vanwambeke SO, Linard C, Soti V. Pathogenic landscapes: Interactions between land, people, disease vectors, and their animal hosts. International Journal of Health Geographics. 2010;9: 54. doi:10.1186/1476-072X-9-54 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1476-072X-9-54&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20979609&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) 33. 33.Beck LR, Lobitz BM, Wood BL. Remote sensing and human health: new sensors and new opportunities. Emerg Infect Dis. 2000;6: 217–227. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid0603.000301&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10827111&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000087321300001&link_type=ISI) 34. 34.Ford TE, Colwell RR, Rose JB, Morse SS, Rogers DJ, Yates TL. Using Satellite Images of Environmental Changes to Predict Infectious Disease Outbreaks. Emerg Infect Dis. 2009;15: 1341–1346. doi:10.3201/eid/1509.081334 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid/1509.081334&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19788799&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269507500001&link_type=ISI) 35. 35.Lleo MM, Lafaye M, Guell A. Application of space technologies to the surveillance and modelling of waterborne diseases. Curr Opin Biotechnol. 2008;19: 307–312. doi:10.1016/j.copbio.2008.04.001 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.copbio.2008.04.001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18514503&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2019%2F10%2F25%2F19005884.atom) 36. 36.Pokharel S, White LJ, Aguas R, Celhay O, Pellé KG, Dittrich S. Algorithm in the diagnosis of febrile illness using pathogen-specific rapid diagnostic tests. Clin Infect Dis. doi:10.1093/cid/ciz665 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciz665&link_type=DOI) 37. 37.Results of the Population and Housing Census 2015 [Internet]. Lao Statistics Bureau; 2016 p. 282. Available: [https://lao.unfpa.org/sites/default/files/pub-pdf/PHC-ENG-FNAL-WEB\_0.pdf](https://lao.unfpa.org/sites/default/files/pub-pdf/PHC-ENG-FNAL-WEB_0.pdf)