Bayesian spatial modeling of childhood overweight and obesity prevalence in Costa Rica.

BACKGROUND: Childhood overweight and obesity levels are rising and becoming a concern globally. In Costa Rica, the prevalence of these conditions has reached alarming values. Spatial analyses can identify risk factors and geographical patterns to develop tailored and effective public health actions in this context. METHODS: A Bayesian spatial mixed model was built to understand the geographic patterns of childhood overweight and obesity prevalence in Costa Rica and their association with some socioeconomic factors. Data was obtained from the 2016 Weight and Size Census (6 - 12 years old children) and 2011 National Census. RESULTS: Average years of schooling increase the levels of overweight and obesity until reaching an approximate value of 8 years, then they start to decrease. Moreover, for every 10-point increment in the percentage of homes with difficulties to cover their basic needs and in the percentage of population under 14 years old, there is a decrease of 7.7 and 14.0 points, respectively, in the odds of obesity. Spatial patterns show higher values of prevalence in the center area of the country, touristic destinations, head of province districts and in the borders with Panama. CONCLUSIONS: Especially for childhood obesity, the average years of schooling is a non-linear factor, describing a U-inverted curve. Lower percentages of households in poverty and population under 14 years old are slightly associated with higher levels of obesity. Districts with high commercial and touristic activity present higher prevalence risk.

Mortality risk attributable to high and low ambient temperature in Pune city, India: A time series analysis from 2004 to 2012.

BACKGROUND: Exposure to high and low ambient temperatures is associated with morbidity and mortality across the globe. Most of these studies assessing the effects of non-optimum temperatures on health and have been conducted in the developed world, whereas in India, the limited evidence on ambient temperature and health risks and has focused mostly on the effects of heat waves. Here we quantify short term association between all temperatures and mortality in urban Pune, India. METHODS: We applied a time series regression model to derive temperature-mortality associations based on daily mean temperature and all-cause mortality records of Pune city from year January 2004 to December 2012. We estimated high and low temperature-mortality relationships by using standard time series quasi-Poisson regression in conjunction with a distributed lag non-linear model (DLNM). We calculated temperature attributable mortality fractions for total heat and total cold. FINDINGS: The analysis provides estimates of the total mortality burden attributable to ambient temperature. Overall, 6∙5% [95%CI 1.76-11∙43] of deaths registered in the observational period were attributed to non-optimal temperatures, cold effect was greater 5.72% [95%CI 0∙70-10∙06] than heat 0∙84% [0∙35-1∙34]. The gender stratified analysis revealed that the highest burden among men both for heat and cold. CONCLUSION: Non-optimal temperatures are associated with a substantial mortality burden. Our findings could benefit national, and local communities in developing preparedness and prevention strategies to reduce weather-related impacts immediately due to climate change.

The effect of community-driven larval source management and house improvement on malaria transmission when added to the standard malaria control strategies in Malawi: a cluster-randomized controlled trial.

BACKGROUND: Current standard interventions are not universally sufficient for malaria elimination. The effects of community-based house improvement (HI) and larval source management (LSM) as supplementary interventions to the Malawi National Malaria Control Programme (NMCP) interventions were assessed in the context of an intensive community engagement programme. METHODS: The study was a two-by-two factorial, cluster-randomized controlled trial in Malawi. Village clusters were randomly assigned to four arms: a control arm; HI; LSM; and HI + LSM. Malawi NMCP interventions and community engagement were used in all arms. Household-level, cross-sectional surveys were conducted on a rolling, 2-monthly basis to measure parasitological and entomological outcomes over 3 years, beginning with one baseline year. The primary outcome was the entomological inoculation rate (EIR). Secondary outcomes included mosquito density, Plasmodium falciparum prevalence, and haemoglobin levels. All outcomes were assessed based on intention to treat, and comparisons between trial arms were conducted at both cluster and household level. RESULTS: Eighteen clusters derived from 53 villages with 4558 households and 20,013 people were randomly assigned to the four trial arms. The mean nightly EIR fell from 0.010 infectious bites per person (95% CI 0.006-0.015) in the baseline year to 0.001 (0.000, 0.003) in the last year of the trial. Over the full trial period, the EIR did not differ between the four trial arms (p = 0.33). Similar results were observed for the other outcomes: mosquito density and P. falciparum prevalence decreased over 3 years of sampling, while haemoglobin levels increased; and there were minimal differences between the trial arms during the trial period. CONCLUSIONS: In the context of high insecticide-treated bed net use, neither community-based HI, LSM, nor HI + LSM contributed to further reductions in malaria transmission or prevalence beyond the reductions observed over two years across all four trial arms. This was the first trial, as far as the authors are aware, to test the potential complementary impact of LSM and/or HI beyond levels achieved by standard interventions. The unexpectedly low EIR values following intervention implementation indicated a promising reduction in malaria transmission for the area, but also limited the usefulness of this outcome for measuring differences in malaria transmission among the trial arms. Trial registration PACTR, PACTR201604001501493, Registered 3 March 2016, https://pactr.samrc.ac.za/ .

COVID-19 pandemic and associated lockdown as a "Global Human Confinement Experiment" to investigate biodiversity conservation.

Efforts to curtail the spread of the novel coronavirus (SARS-CoV2) have led to the unprecedented concurrent confinement of nearly two-thirds of the global population. The large human lockdown and its eventual relaxation can be viewed as a Global Human Confinement Experiment. This experiment is a unique opportunity to identify positive and negative effects of human presence and mobility on a range of natural systems, including wildlife, and protected areas, and to study processes regulating biodiversity and ecosystems. We encourage ecologists, environmental scientists, and resource managers to contribute their observations to efforts aiming to build comprehensive global understanding based on multiple data streams, including anecdotal observations, systematic assessments and quantitative monitoring. We argue that the collective power of combining diverse data will transcend the limited value of the individual data sets and produce unexpected insights. We can also consider the confinement experiment as a "stress test" to evaluate the strengths and weaknesses in the adequacy of existing networks to detect human impacts on natural systems. Doing so will provide evidence for the value of the conservation strategies that are presently in place, and create future networks, observatories and policies that are more adept in protecting biological diversity across the world.

Extended excess hazard models for spatially dependent survival data.

Relative survival represents the preferred framework for the analysis of population cancer survival data. The aim is to model the survival probability associated with cancer in the absence of information about the cause of death. Recent data linkage developments have allowed for incorporating the place of residence into the population cancer databases; however, modeling this spatial information has received little attention in the relative survival setting. We propose a flexible parametric class of spatial excess hazard models (along with inference tools), named "Relative Survival Spatial General Hazard," that allows for the inclusion of fixed and spatial effects in both time-level and hazard-level components. We illustrate the performance of the proposed model using an extensive simulation study, and provide guidelines about the interplay of sample size, censoring, and model misspecification. We present a case study using real data from colon cancer patients in England. This case study illustrates how a spatial model can be used to identify geographical areas with low cancer survival, as well as how to summarize such a model through marginal survival quantities and spatial effects.

Spatial analysis of modern contraceptive use among women who need it in Ethiopia: Using geo-referenced data from performance monitoring for action.

INTRODUCTION: The challenge of achieving maternal and neonatal health-related goals in developing countries is significantly impacted by high fertility rates, which are partly attributed to limited access to family planning and access to the healthcare systems. The most widely used indicator to monitor family planning coverage is the proportion of women in reproductive age using contraception (CPR). However, this metric does not accurately reflect the true family planning coverage, as it fails to account for the diverse needs of women in reproductive age. Not all women in this category require contraception, including those who are pregnant, wish to become pregnant, sexually inactive, or infertile. To effectively address the contraceptive needs of those who require it, this study aims to estimate family planning coverage among this specific group. Further, we aimed to explore the geographical variation and factors influencing contraceptive uptake of contraceptive use among those who need. METHOD: We used data from the Performance Monitoring for Action Ethiopia (PMA Ethiopia) survey of women of reproductive age and the service delivery point (SDP) survey conducted in 2019. A total of 4,390 women who need contraception were considered as the analytical sample. To account for the study design, sampling weights were considered to compute the coverage of modern contraceptive use disaggregated by socio-demographic factors. Bayesian geostatistical modeling was employed to identify potential factors associated with the uptake of modern contraception and produce spatial prediction to unsampled locations. RESULT: The overall weighted prevalence of modern contraception use among women who need it was 44.2% (with 95% CI: 42.4%-45.9%). Across regions of Ethiopia, contraceptive use coverage varies from nearly 0% in Somali region to 52.3% in Addis Ababa. The average nearest distance from a woman's home to the nearest SDP was high in the Afar and Somali regions. The spatial mapping shows that contraceptive coverage was lower in the eastern part of the country. At zonal administrative level, relatively high (above 55%) proportion of modern contraception use coverage were observed in Adama Liyu Zone, Ilu Ababor, Misrak Shewa, and Kefa zone and the coverage were null in majority of Afar and Somali region zones. Among modern contraceptive users, use of the injectable dominated the method-mix. The modeling result reveals that, living closer to a SDP, having discussions about family planning with the partner, following a Christian religion, no pregnancy intention, being ever pregnant and being young increases the likelihood of using modern contraceptive methods. CONCLUSION: Areas with low contraceptive coverage and lower access to contraception because of distance should be prioritized by the government and other supporting agencies. Women who discussed family planning with their partner were more likely to use modern contraceptives unlike those without such discussion. Thus, to improve the coverage of contraceptive use, it is very important to encourage/advocate women to have discussions with their partner and establish movable health systems for the nomadic community.

Spatio-temporal modeling of infectious diseases by integrating compartment and point process models.

Infectious disease modeling plays an important role in understanding disease spreading dynamics and can be used for prevention and control. The well-known SIR (Susceptible, Infected, and Recovered) compartment model and spatial and spatio-temporal statistical models are common choices for studying problems of this kind. This paper proposes a spatio-temporal modeling framework to characterize infectious disease dynamics by integrating the SIR compartment and log-Gaussian Cox process (LGCP) models. The method's performance is assessed via simulation using a combination of real and synthetic data for a region in São Paulo, Brazil. We also apply our modeling approach to analyze COVID-19 dynamics in Cali, Colombia. The results show that our modified LGCP model, which takes advantage of information obtained from the previous SIR modeling step, leads to a better forecasting performance than equivalent models that do not do that. Finally, the proposed method also allows the incorporation of age-stratified contact information, which provides valuable decision-making insights. Supplementary Information: The online version contains supplementary material available at 10.1007/s00477-022-02354-4.

Joint spatial modeling of the risks of co-circulating mosquito-borne diseases in Ceará, Brazil.

Mosquito-borne diseases such as dengue and chikungunya have been co-circulating in the Americas, causing great damage to the population. In 2021, for instance, almost 1.5 million cases were reported on the continent, being Brazil the responsible for most of them. Even though they are transmitted by the same mosquito, it remains unclear whether there exists a relationship between both diseases. In this paper, we model the geographic distributions of dengue and chikungunya over the years 2016 to 2021 in the Brazilian state of Ceará. We use a Bayesian hierarchical spatial model for the joint analysis of two arboviruses that includes spatial covariates as well as specific and shared spatial effects that take into account the potential autocorrelation between the two diseases. Our findings allow us to identify areas with high risk of one or both diseases. Only 7% of the areas present high relative risk for both diseases, which suggests a competition between viruses. This study advances the understanding of the geographic patterns and the identification of risk factors of dengue and chikungunya being able to help health decision-making.

Spatial modeling and socioeconomic inequities of COVID-19 in the urban area of the city of Cali, Colombia.

COVID-19 has spread worldwide with a high variability in cases and mortality between populations. This research aims to assess socioeconomic inequities of COVID-19 in the city of Cali, Colombia, during the first and second peaks of the pandemic in this city. An ecological study by neighborhoods was carried out, were COVID-19 cases were analyzed using a Bayesian hierarchical spatial model that includes potential risk factors such as the index of unsatisfied basic needs and socioeconomic variables as well as random effects to account for residual variation. Maps showing the geographic patterns of the estimated relative risks as well as exceedance probabilities were created. The results indicate that in the first wave, the neighborhoods with the greatest unsatisfied basic needs and low socioeconomic strata, were more likely to report positive cases for COVID-19. For the second wave, the disease begins to spread through different neighborhoods of the city and middle socioeconomic strata presents the highest risk followed by the lower strata. These findings indicate the importance of measuring social determinants in the study of the distribution of cases due to COVID-19 for its inclusion in the interventions and measures implemented to contain contagions and reduce impacts on the most vulnerable populations.

The Changing Trend of Paediatric Emergency Department Visits in Malaysia Following the COVID-19 Pandemic.

Background The coronavirus disease 2019 (COVID-19) pandemic has impacted the emergency department (ED) due to the surge in medical demand and changes in the characteristics of paediatric visits. Additionally, the trend for paediatric ED visits has decreased globally, secondary to implementing lockdowns to stop the spread of COVID-19. We aim to study the trend and characteristics of paediatric ED visits following Malaysia's primary timeline of the COVID-19 pandemic. Methods and materials A five-year time series observational study of paediatric ED patients from two tertiary hospitals in Malaysia was conducted from March 17, 2017 (week 11 2017) to March 17, 2022 (week 12 2022). Aggregated weekly data were analysed using R statistical software version 4.2.2 (R Foundation for Statistical Computing, Vienna, Austria) against significant events during the COVID-19 pandemic to detect influential changepoints in the trend. The data collected were the number of ED visits, triage severity, visit outcomes and ED discharge diagnosis. Results Overall, 175,737 paediatric ED visits were recorded with a median age of three years and predominantly males (56.8%). A 57.57% (p<0.00) reduction in the average weekly ED visits was observed during the Movement Control Order (MCO) period. Despite the increase in the proportion of urgent (odds ratio (OR): 1.23, p<0.00) and emergent or life-threatening (OR: 1.79, p<0.00) cases, the proportion of admissions decreased. Whilst the changepoints during the MCO indicated a rise in respiratory, fever or other infectious diseases, or gastrointestinal conditions, diagnosis of complications originating from the perinatal period declined from July 19, 2021 (week 29 2021). Conclusion The incongruent change in disease severity and hospital admission reflects the potential effects of the healthcare system reform and socioeconomic impact as the pandemic evolves. Future studies on parental motivation to seek emergency medical attention may provide insight into the timing and choice of healthcare service utilisation.

rspatialdata: a collection of data sources and tutorials on downloading and visualising spatial data using R.

Spatial and spatio-temporal data are used in a wide range of fields including environmental, health and social disciplines. Several packages in the statistical software R have been recently developed as clients for various databases to meet the growing demands for easily accessible and reliable spatial data. While documentation on how to use many of these packages exist, there is an increasing need for a one stop repository for tutorials on this information. In this paper, we present  rspatialdata a website that provides a collection of data sources and tutorials on downloading and visualising spatial data using R. The website includes a wide range of datasets including administrative boundaries of countries, Open Street Map data, population, temperature, vegetation, air pollution, and malaria data. The goal of the website is to equip researchers and communities with the tools to engage in spatial data analysis and visualisation so that they can address important local issues, such as estimating air pollution, quantifying disease burdens, and evaluating and monitoring the United Nation's sustainable development goals.

Modeling infectious disease dynamics: Integrating contact tracing-based stochastic compartment and spatio-temporal risk models.

Major infectious diseases such as COVID-19 have a significant impact on population lives and put enormous pressure on healthcare systems globally. Strong interventions, such as lockdowns and social distancing measures, imposed to prevent these diseases from spreading, may also negatively impact society, leading to jobs losses, mental health problems, and increased inequalities, making crucial the prioritization of riskier areas when applying these protocols. The modeling of mobility data derived from contact-tracing data can be used to forecast infectious trajectories and help design strategies for prevention and control. In this work, we propose a new spatial-stochastic model that allows us to characterize the temporally varying spatial risk better than existing methods. We demonstrate the use of the proposed model by simulating an epidemic in the city of Valencia, Spain, and comparing it with a contact tracing-based stochastic compartment reference model. The results show that, by accounting for the spatial risk values in the model, the peak of infected individuals, as well as the overall number of infected cases, are reduced. Therefore, adding a spatial risk component into compartment models may give finer control over the epidemic dynamics, which might help the people in charge to make better decisions.

Local mortality impacts due to future air pollution under climate change scenarios.

The health impacts of global climate change mitigation will affect local populations differently. However, most co-benefits analyses have been done at a global level, with relatively few studies providing local level results. We aimed to quantify the local health impacts due to fine particles (PM2.5) under the governance arrangements embedded in the Shared Socioeconomic Pathways (SSPs1-5) under two greenhouse gas concentration scenarios (Representative Concentration Pathways (RCPs) 2.6 and 8.5) in local populations of Mozambique, India, and Spain. We simulated the SSP-RCP scenarios using the Global Change Analysis Model, which was linked to the TM5-FASST model to estimate PM2.5 levels. PM2.5 levels were calibrated with local measurements. We used comparative risk assessment methods to estimate attributable premature deaths due to PM2.5 linking local population and mortality data with PM2.5-mortality relationships from the literature, and incorporating population projections under the SSPs. PM2.5 attributable burdens in 2050 differed across SSP-RCP scenarios, and sensitivity of results across scenarios varied across populations. Future attributable mortality burden of PM2.5 was highly sensitive to assumptions about how populations will change according to SSP. SSPs reflecting high challenges for adaptation (SSPs 3 and 4) consistently resulted in the highest PM2.5 attributable burdens mid-century. Our analysis of local PM2.5 attributable premature deaths under SSP-RCP scenarios in three local populations highlights the importance of both socioeconomic development and climate policy in reducing the health burden from air pollution. Sensitivity of future PM2.5 mortality burden to SSPs was particularly evident in low- and middle- income country settings due either to high air pollution levels or dynamic populations.

SARS-CoV-2 genomes from Saudi Arabia implicate nucleocapsid mutations in host response and increased viral load.

Monitoring SARS-CoV-2 spread and evolution through genome sequencing is essential in handling the COVID-19 pandemic. Here, we sequenced 892 SARS-CoV-2 genomes collected from patients in Saudi Arabia from March to August 2020. We show that two consecutive mutations (R203K/G204R) in the nucleocapsid (N) protein are associated with higher viral loads in COVID-19 patients. Our comparative biochemical analysis reveals that the mutant N protein displays enhanced viral RNA binding and differential interaction with key host proteins. We found increased interaction of GSK3A kinase simultaneously with hyper-phosphorylation of the adjacent serine site (S206) in the mutant N protein. Furthermore, the host cell transcriptome analysis suggests that the mutant N protein produces dysregulated interferon response genes. Here, we provide crucial information in linking the R203K/G204R mutations in the N protein to modulations of host-virus interactions and underline the potential of the nucleocapsid protein as a drug target during infection.

Detection of spatial disease clusters with LISA functions.

Detection of disease clusters is an important tool in epidemiology that can help to identify risk factors associated with the disease and in understanding its etiology. In this article we propose a method for the detection of spatial clusters where the locations of a set of cases and a set of controls are available. The method is based on local indicators of spatial association functions (LISA functions), particularly on the development of a local version of the product density, which is a second-order characteristic of spatial point processes. The behavior of the method is evaluated and compared with Kulldorff's spatial scan statistic by means of a simulation study. It is shown that the LISA method yields high sensitivity and specificity when it is used to detect simulated clusters of different sizes and shapes. It also performs better than the spatial scan statistic when they are used to detect clusters of irregular shape; however, it presents relatively high type I error in situations where the number of cases is high. Both methods are applied for detecting spatial clusters of kidney disease in the city of Valencia, Spain, in the year 2008.

Bayesian spatial modelling of geostatistical data using INLA and SPDE methods: A case study predicting malaria risk in Mozambique.

Bayesian spatial models are widely used to analyse data that arise in scientific disciplines such as health, ecology, and the environment. Traditionally, Markov chain Monte Carlo (MCMC) methods have been used to fit these type of models. However, these are highly computationally intensive methods that present a wide range of issues in terms of convergence and can become infeasible in big data problems. The integrated nested Laplace approximation (INLA) method is a computational less-intensive alternative to MCMC that allows us to perform approximate Bayesian inference in latent Gaussian models such as generalised linear mixed models and spatial and spatio-temporal models. This approach can be used in combination with the stochastic partial differential equation (SPDE) approach to analyse geostatistical data that have been collected at particular sites to predict the spatial process underlying the data as well as to assess the effect of covariates and model other sources of variability. Here we demonstrate how to fit a Bayesian spatial model using the INLA and SPDE approaches applied to freely available data of malaria prevalence and risk factors in Mozambique. We show how to fit and interpret the model to predict malaria risk and assess the effect of covariates using the R-INLA package, and provide the R code necessary to reproduce the results or to use it in other spatial applications.

The emergence and transmission of COVID-19 in European countries, 2019-2020: a comprehensive review of timelines, cases and containment.

When it emerged in late 2019, COVID-19 was carried via travelers to Germany, France and Italy, where freedom of movement accelerated its transmission throughout Europe. However, effective non-pharmaceutical interventions introduced by European governments led to containment of the rapid increase in cases within European nations. Electronic searches were performed to obtain the number of confirmed cases, incident rates and non-pharmaceutical government measures for each European country. The spread and impact of non-pharmaceutical interventions throughout Europe were assessed and visualized. Specifically, heatmaps were used to represent the number of confirmed cases and incident rates for each of the countries over time. In addition, maps were created showing the number of confirmed cases and incident rates in Europe on three different dates (15 March, 15 April and 15 May 2020), which allowed us to assess the geographic and temporal patterns of the disease.

Identifying Plasmodium falciparum transmission patterns through parasite prevalence and entomological inoculation rate.

Background: Monitoring malaria transmission is a critical component of efforts to achieve targets for elimination and eradication. Two commonly monitored metrics of transmission intensity are parasite prevalence (PR) and the entomological inoculation rate (EIR). Comparing the spatial and temporal variations in the PR and EIR of a given geographical region and modelling the relationship between the two metrics may provide a fuller picture of the malaria epidemiology of the region to inform control activities. Methods: Using geostatistical methods, we compare the spatial and temporal patterns of Plasmodium falciparum EIR and PR using data collected over 38 months in a rural area of Malawi. We then quantify the relationship between EIR and PR by using empirical and mechanistic statistical models. Results: Hotspots identified through the EIR and PR partly overlapped during high transmission seasons but not during low transmission seasons. The estimated relationship showed a 1-month delayed effect of EIR on PR such that at lower levels of EIR, increases in EIR are associated with rapid rise in PR, whereas at higher levels of EIR, changes in EIR do not translate into notable changes in PR. Conclusions: Our study emphasises the need for integrated malaria control strategies that combine vector and human host managements monitored by both entomological and parasitaemia indices. Funding: This work was supported by Stichting Dioraphte grant number 13050800.

Assessing the age- and gender-dependence of the severity and case fatality rates of COVID-19 disease in Spain.

Background: The assessment of the severity and case fatality rates of coronavirus disease 2019 (COVID-19) and the determinants of its variation is essential for planning health resources and responding to the pandemic. The interpretation of case fatality rates (CFRs) remains a challenge due to different biases associated with surveillance and reporting. For example, rates may be affected by preferential ascertainment of severe cases and time delay from disease onset to death. Using data from Spain, we demonstrate how some of these biases may be corrected when estimating severity and case fatality rates by age group and gender, and identify issues that may affect the correct interpretation of the results. Methods: Crude CFRs are estimated by dividing the total number of deaths by the total number of confirmed cases. CFRs adjusted for preferential ascertainment of severe cases are obtained by assuming a uniform attack rate in all population groups, and using demography-adjusted under-ascertainment rates. CFRs adjusted for the delay between disease onset and death are estimated by using as denominator the number of cases that could have a clinical outcome by the time rates are calculated. A sensitivity analysis is carried out to compare CFRs obtained using different levels of ascertainment and different distributions for the time from disease onset to death. Results: COVID-19 outcomes are highly influenced by age and gender. Different assumptions yield different CFR values but in all scenarios CFRs are higher in old ages and males. Conclusions: The procedures used to obtain the CFR estimates require strong assumptions and although the interpretation of their magnitude should be treated with caution, the differences observed by age and gender are fundamental underpinnings to inform decision-making.

Spatiotemporal patterns and integrated approach to prioritize areas for surveillance and control of visceral leishmaniasis in a large metropolitan area in Brazil.

Visceral leishmaniasis (VL) is a neglected disease with worldwide distribution. Brazil is the country with the largest number of cases in the Americas, and the state of Minas Gerais presents a high VL-related burden and a high case fatality rate. We aimed to analyse the spatial and spatiotemporal patterns of VL occurrence and to identify priority risk areas for surveillance and control in the metropolitan region of Belo Horizonte-MG, the third largest metropolitan area in Brazil. An ecological study was conducted considering all cases of VL in humans confirmed from 2006 to 2017. The crude and smoothed incidence rates were used to analyse the distribution patterns of the disease (dispersed, random, or clustered) based on global and local indicators of spatial association and space-time risk assessment. Positive spatial autocorrelation and spatial dependence were found between incidence rates. It was possible to observe a high concentration of VL cases in the metropolitan core area, with the identification of two high-risk clusters in strictly urban areas, showing an urban association with the disease. Ten municipalities were categorised as high risk for VL occurrence. Our results provide evidence for making decisions in surveillance programs, suggesting the prioritisation of the municipalities with more risk of transmission.