Global transmission suitability maps for dengue virus transmitted by Aedes aegypti from 1981 to 2019 ==================================================================================================== * Taishi Nakase * Marta Giovanetti * Uri Obolski * José Lourenço ## ABSTRACT Mosquito-borne viruses increasingly threaten human populations due to accelerating changes in climate, human and mosquito migration, and land use practices. Over the last three decades, the global distribution of dengue has rapidly expanded, causing detrimental health and economic problems in many areas of the world. To develop effective disease control measures and plan for future epidemics, there is an urgent need to map the current and future transmission potential of dengue across both endemic and emerging areas. Expanding and applying Index P, a previously developed mosquito-borne viral suitability measure, we map the global climate-driven transmission potential of dengue virus transmitted by *Aedes aegypti* mosquitoes from 1981 to 2019. This database of dengue transmission suitability maps and an R package for Index P estimations are offered to the public health community as resources towards the identification of past, current and future transmission hotspots. These resources and the studies they facilitate can contribute to the planning of disease control and prevention strategies, especially in areas where surveillance is unreliable or non-existent. ## Background & Summary Global changes in environmental conditions that favor a closer proximity between vector and human populations are facilitating the spread of mosquito-borne viruses (MBVs)1. These processes, which include climate change, urbanization, deforestation and migration are interacting and evolving, creating a complex landscape of current and future disease risk. For example, the recent emergence and epidemic spread of Zika virus (ZIKV) caused severe health and economic problems in Latin America2. The public health impact of chikungunya (CHIKV) has also significantly increased in recent years with its rapid emergence in the Americas in late 20133–5. The threat of these events is further magnified by the potential of their pathogens to cause severe health complications such as dengue haemorrhagic fever, ZIKV-induced congenital and neurological disorders, and CHIKV-associated Guillain-Barré Syndrome6, 7. Dengue virus (DENV), transmitted by *Aedes spp*. mosquitoes, accounts for an increasingly large proportion of all vector-borne disease, with an estimated 400 million infections per year8, 9. The geographical range of DENV has expanded in recent years, becoming endemic in much of Central America, South America and Southeast Asia and establishing epidemic cycles in parts of Africa and North America10. In the last decade, Europe has also seen a rise in autochthonous DENV transmission as mosquito populations advance north from the Mediterranean and travel to endemic areas increases11. The only licensed dengue vaccine remains in limited use because of the higher risk of severe disease from breakthrough infections in seronegative recipients12. Current control efforts are focused on mosquito control and viral surveillance in human and mosquito populations13. Given the increasing threat posed by DENV and the reliance on vector control to mitigate transmission, understanding and mapping the transmission potential of DENV by *Aedes aegypti* mosquitoes, its primary vector, is essential for control and health systems planning. The spatiotemporal dynamics of MBVs are primarily governed by the interplay among three factors: the physiological interactions between virus and vector, that between virus and host, and the population dynamics of the vector and host. The carrying capacity and seasonal oscillations of mosquito populations are influenced by many environmental and ecological factors including temperature14–16, humidity16, host population density16, 17, vegetation types18 and the presence of water sources16, 17. Climatic factors are particularly important because they also alter each mosquito’s potential to transmit the virus to new hosts by causing changes in physical and behavioral traits such as life span19, incubation period20 and biting rate21. With data from experimental studies, it is possible to mathematically characterize these relationships between climatic variables and mosquito-viral traits and incorporate them into mechanistic transmission models22. The basic reproductive number (*R*), which measures the average number of secondary infections generated by a single infectious host in a fully susceptible host population, is an important measure of transmission potential. For MBVs, expressions of *R* involve many interacting host, pathogen and vector variables, some of which are difficult to parameterize due to limited data (e.g. geographical mosquito distribution and density). To overcome this challenge, different measures of transmission potential, called suitability measures, have been developed18, 21, 23, 24. Vectorial capacity, one of the most widely known suitability measures, uses purely entomological variables to estimate the potential number of infectious bites that would arise from a single infectious person over one day23. A number of studies have built upon the concepts in vectorial capacity to develop new measures of transmission potential that consider the interactions between climate and mosquito-viral traits22, 24, 25. Here, we update and apply a previously developed mosquito-borne viral suitability measure24, referred to as Index P, to estimate geographical maps and time series for climate-driven DENV transmission potential of the *Ae. aegypti* vector. Index P is a proxy for the transmission potential of a single adult female mosquito under conditions where susceptible hosts, the virus and its vectors are assumed to be present26. Unlike previously developed suitability measures18, 25, Index P has a direct biological interpretation and takes into account factors that are not purely entomological (e.g. human infectious period 1*/σ**h* and the transmission probability from infected human to mosquito *ϕ**h*→*v*). Local temperature and humidity time series are its main inputs, making its framework sufficiently general to be applied to any location for which such climate data exists and for any MBV for which there is empirical data on the relationship between climate and vector-viral traits. It has been successfully used to characterize the transmission potential and epidemiology of West Nile virus in Israel27, Portugal28, Brazil29, and Italy30, CHIKV and ZIKV in the Dominican Republic5 and Mexico31, and DENV in Myanmar26, Brazil24, 32and Mexico31. So far, the application of Index P has been performed at a small scale driven by specific and limited research goals. Given the sparsity of detailed epidemiological and ecological data for DENV (e.g. infection incidence and mosquito density) and the disproportionate impact of DENV in developing countries, there is a need to extend these analyses globally. Hence, we offer a database of spatiotemporal maps of Index P for 186 different countries and territories and an easy-to-use R package for Index P estimation as ready-to-use resources for the visualization and analysis of climate-driven DENV transmission potential globally over the last four decades. In combination with previous work to map the distribution of DENV vectors18, 33, 34, these maps improve our understanding of past and current DENV transmission potential and have the capacity to inform disease control and prevention strategies. For example, in Central Africa where *Ae. aegypti* is predicted to be widely disseminated18 but reported DENV cases remain low35, these maps can help determine whether the apparent mismatch is due to low transmission potential of the vector or under-reporting and poor surveillance. This information can then be used to identify high-risk, low-surveillance areas where seroprevalence surveys should be targeted. It thus addresses a central limitation of occurrence-based predictions of dengue infection risk8, which can underestimate risk in areas where reported incidence is low or absent due to limited surveillance rather than low transmission. The maps may also be used to identify regions where *Ae. aegypti* might not be present but due to high DENV transmission potential there is a high likelihood of spillover of sylvatic DENV into human populations36, 37. ## Methods ### Fundamental theory of Index P Index P is a climate-driven suitability measure for mosquito-borne viruses derived from a mechanistic model previously developed to study the transmission dynamics of Zika38 and dengue26. The equation for *R* in this model (equation (1)) can be decomposed into two components: the number of female mosquitoes per host (*M*) and the transmission potential of each female mosquito (*P*). Index P measures transmission potential solely based on *P* (equation (2)), which is interpreted as the reproductive potential of a single adult female mosquito in a fully susceptible host population. A description of the parameters in the equation for Index P can be found in Table 1. Given the availability of empirical studies that quantify the relationships between meteorological variables and viral-vector traits (e.g. extrinsic incubation period), it is possible to define equations for several parameters in the expression for *P* in terms of temperature (*t*) and humidity (*u*). Temperature and humidity time series data can be input into these equations to derive time series of the entomological parameters, which when combined with prior information on the selected vector/host/virus system can provide estimates of Index P over time. View this table: [Table 1.](http://medrxiv.org/content/early/2022/11/05/2022.11.04.22281958/T1) Table 1. Descriptions and probability distributions of the biological parameters used in the estimation of Index P for DENV transmitted by Aedes aegypti mosquitoes. SD is standard deviation. The transmission probability per mosquito bite (*Ae. aegypti*-to-human and human-to-*Ae. aegypti*) are defined by deterministic temperature-dependent equations. ![Formula][1] where ![Formula][2] ### Proposal of informative distributions of the biological parameters The estimation of Index P requires some prior knowledge of the biology of the selected host, vector and viral system. In this study, we defined probability distributions for human hosts, *Ae. aegypti* vectors and dengue virus. Several mosquito species in the genus *Aedes* are known to transmit DENV to humans. We focused on *Ae. aegypti* because it accounts for the majority of vector-human transmission39, and it is the species for which the most empirical data exists on the relationship between climate and vector-viral traits. After a review of relevant literature, we defined probability distributions for six host-virus and vector-virus parameters: extrinsic *Ae. aegypti*-DENV incubation period, adult *Ae. aegypti* lifespan, adult *Ae. aegypti* biting rate, human lifespan, intrinsic human-DENV incubation period, and human-DENV infectious period (Table 1). These probability distributions are either directly sampled (*µ**h*, *γ**h*, *σ**h* and ![Graphic][3]) or used as likelihood functions for the estimation of climate-dependent probability distributions (![Graphic][4] and ![Graphic][5]). ### Climate-dependent functions for entomological parameters Weather-dependent functions are defined for each of the entomological parameters in the expression for Index P (equation (2)). The adult human mortality rate (*µ**h*), intrinsic incubation period (1/*γ**h*) and the human infectious period (1/*σ**h*) are taken to be climate-independent. As detailed in Obolski et al.24, adult vector mortality (*µ**v*) and the probability of transmission from vector to human per bite (*ϕ**v*→*h*) are modeled by temperature-dependent functions estimated in experiments on laboratory strains of *Ae. aegypti* by Yang et al.40 (equation (3)) and Lambrechts et al.41 (equation (4)), respectively. In the original implementation of Index P, the probability of transmission from human to mosquito per bite (*ϕ**h*→*v*) was assumed to be constant. To account for the temperature dependence of *ϕ**h*→*v* and to capture threshold effects at extreme temperatures, we expanded the framework by modeling *ϕ**h*→*v* with a ramp model estimated by Lambrechts et al.41 (equation (5)). ![Formula][6] ![Formula][7] ![Formula][8] Following the original approach by Obolski et al.24, the influence of humidity on adult mosquito mortality rate (*µ**v*) and mosquito biting rate (*a**v*) is introduced according to the following expressions. The time series for relative humidity is normalized to [0, 1] with the ecological variables centered around the local average (*ū*). ![Formula][9] ![Formula][10] The total effect of climate on the entomological parameters is then modeled by combining the estimated temperature-dependent and humidity-dependent functions. ![Formula][11] ![Formula][12] ![Formula][13] ![Formula][14] The scaling coefficients *η* and *ρ* in the expressions for mosquito biting rate (equation (11)) and adult mosquito mortality rate (equation (8)) are estimated from the user-defined probability distributions of *a**v* and *µ**v*, respectively. These distributions permit deviations from the ideal laboratory conditions used in the original empirical studies. The multiplicative coefficient *η* determines the magnitude by which the effect of temperature on the adult mosquito mortality deviates from that observed under laboratory conditions. The effect of temperature under field conditions is identical to that under laboratory conditions when *η* = 0, stronger when *η >* 1 and weaker when *η <* 1. The exponential coefficient *ρ* modulates the relative influence that humidity has on adult mosquito mortality and mosquito biting rate: humidity has no effect when *ρ* = 0 and a stronger effect when *ρ >* 0. For further information on the biological interpretation of these coefficients and their estimation, please refer to Obolski et al.24. Using recent, high-resolution empirical data, we also updated the probability distribution of the extrinsic incubation period (EIP) of *Ae. aegypti* mosquitoes. We introduced a temperature-dependent probability distribution of *γ**v* derived from observational studies of EIP in human subjects experimentally infected by wild *Ae. aegypti* mosquitoes20. Chan and Johansson20 estimate the temperature-dependent distribution of the EIP of *Ae. aegypti* mosquitoes (1*/γ**v*) by directly fitting a log-normal time-to-event model to EIP observational data: ![Formula][15] Where *α*(*t*) and *β* (*t*) are temperature-dependent functions of the mean and standard derivation of the natural logarithm of the EIP (Table 1). ### Estimation of Index P Index P estimation is based on a Markov chain Monte Carlo (MCMC) sampling framework which derives the distribution of the Index P time series using the climate-driven functions for the entomological parameters (*µ**v*, *ϕ**v*→*h*, *ϕ**h*→*v*, *a**v*), the user-defined biological probability distributions (Table 1), and the climate time series data. Broadly, the estimation of Index P can be broken down into three steps. First, the posterior distributions of *µ**v* and *a**v* are estimated by fitting the scaling coefficients *η* and *ρ* to the corresponding probability distributions proposed by the user and the time series of the climatic variables. For each iteration of the MCMC procedure, the scaling coefficients are sampled and the climate-driven time series of *µ**v* and *a**v* are estimated. The posterior probability of the resultant time series is assessed through the product of the user-defined likelihood functions of *µ**v* and *a**v* with the prior probabilities of those parameters. Then, for each time point in the climate time series, the user-defined distributions of *µ**h*, *γ**h*, *σ**h* and *γ**v* and the derived posterior distributions of *µ**v* and *a**v* are independently sampled, and the temperature-dependent equations for *ϕ**v*→*h* and *ϕ**h*→*v* are solved. Finally, the sampled parameter values are plugged into the Index P expression, which yields a distribution of Index P time series bound according to the user-defined biological characteristics of the virus, vector, and host under investigation. An example of the workflow for the estimation of Index P time series at a given spatial pixel is provided in Figure 1. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/11/05/2022.11.04.22281958/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/11/05/2022.11.04.22281958/F1) Figure 1. Summary of the steps required to estimate Index P at a single spatial pixel. ### Estimation and sampling of Index P globally We used monthly temperature and humidity data, as well as informed probability distributions of the biological parameters as inputs to estimate the monthly Index P time series for each spatial pixel (∼ 28 km2) of 186 individual countries and territories from 1981 to 2019. For each country, spatial polygons were used to curate the monthly temperature and relative humidity data sourced from Copernicus.eu42 for all spatial pixels within its geographical boundary. The following steps were then performed for each spatial pixel. Using the climate data and the user-defined probability distributions for mosquito biting rate (*a**v*) and adult mosquito mortality rate (*µ**v*), we estimated the climate-dependent posterior distributions of these two parameters by fitting the scaling coefficients *η* and *ρ*. For each MCMC chain, we drew 20, 000 samples with the first 20% of samples treated as burn-in. Given a lack of information on the possible estimates of *ρ* and *η*, we assumed uninformative, flat priors for the two scaling coefficients in the ranges 0 − 10 and 0 − 20 respectively. We then drew 1, 000 samples from the climate-independent distributions of *µ**h*, *γ**h* and *σ**h*, 1, 000 samples from the climate-dependent distributions of the *γ**v*, *µ**v* and *a**v* time series, and solved the deterministic temperature-dependent equations of *ϕ**v*→*h* and *ϕ**h*→*v*. These sampled values were plugged into the equation for Index P, resulting in 1, 000 monthly Index P time series. Finally, we calculated various summary statistics of the distribution of monthly Index P time series for each spatial pixel. Details on the calculated summary statistics are provided in the Data Records section. These estimates were collated at the country/territory level into TIFF files43. ## Data Records ### Meteorological data The meteorological data from 1981 to 2019 that were used to estimate the spatiotemporal maps of Index P are freely available from the “Essential Climate Variables for assessment of climate variability from 1979 to present” dataset published and maintained by Copernicus.eu42. We used the average monthly surface air temperature (*K*) and surface air relative humidity (%). The spatial pixel resolution used was 0.25° ×0.25° (∼28 km2), the maximum resolution available for the Copernicus.eu dataset. ### DENV case data Weekly notified dengue case data per municipality in Brazil for the years 2000 to 2014 were provided by the Brazilian Ministry of Health ([https://www.gov.br/pt-br](https://www.gov.br/pt-br)) and by the SINAN (Sistema de Informação de Agravos de Notificação; [http://portalsinan.saude.gov.br/](http://portalsinan.saude.gov.br/))44. It includes clinically suspected infections without laboratory confirmation for 5,570 municipalities. Four municipalities had incomplete time series data and were excluded from the analysis: Mojuí dos Campos (Pará), Pescaria Brava (Santa Catarina), Balneário Rincão (Santa Catarina) and Paraíso das Águas (Mato Grosso do Sul). ### Population data The yearly population size per municipality in Brazil for the years 2001 to 2014 was sourced from IBGE - Instituto Brasileiro de Geografia e Estatística ([https://www.ibge.gov.br/en/home-eng.html](https://www.ibge.gov.br/en/home-eng.html))45. Missing population data for the year 2000 was assumed to be the same as that for the year 2001. ### Global Index P data For each country or territory, we calculated epidemiologically relevant summary statistics for the spatiotemporal distribution of Index P. These include: spatiotemporal map(s) for each month (468 layers), each year (39 layers), a typical year (12 layers) and the entire period (1 layer); spatiotemporal maps of the number of months Index P is above 1.0 for each year (39 layers), and for a typical year (12 layers); and spatiotemporal maps of the peak and trough timing during a typical year (12 layers each). The spatiotemporal maps are provided as TIFF files43. The Index P spatiotemporal maps are projected with World Geodetic System (WGS)84, latitude/longitude coordinate system (EPSG: 4326; [https://epsg.io/4326](https://epsg.io/4326)). The complete dataset and a description of its contents are available from a figshare repository ([http://www.figshare.com/](http://www.figshare.com/))43. The Index P maps can be processed with packages such as *raster*46 in the R programming language. A sample of the Index P maps and time series is in Figure 2. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/11/05/2022.11.04.22281958/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/11/05/2022.11.04.22281958/F2) Figure 2. Summary of global spatiotemporal maps of estimated Index P for DENV transmitted by *Aedes aegypti* mosquitoes. Global map of mean Index P during a typical year at a spatial resolution of 0.25° ×0.25° (∼28km2). Includes monthly time series of average monthly Index P from 2000 to 2019 for six cities: Bangkok, Thailand; Brisbane, Australia; Dakar, Senegal; Miami, United States; Rio de Janeiro, Brazil; and Rome, Italy. Cities were selected to capture the diversity of Index P dynamics observed globally across regions where DENV is endemic, emerging or absent. ### Technical Validation Using historical DENV case data recorded in Brazil from 2000 to 2014, we performed a validation of Index P as a climate-based measure of transmission potential. We focused on Brazil for several reasons: (1) DENV is endemic with high disease burden8, 47, (2) the country is climatically diverse with both temperate and tropical areas48 and there is large variation in the intensity and seasonality of observed DENV dynamics49. Given that DENV incidence data are available at the municipality level for Brazil, we calculated and validated the relevant summary statistics of Index P at the municipality level rather than at the spatial pixel level. Here, we included several analyses that support the application of Index P to explain past spatiotemporal dynamics of DENV, reaffirming the results of previous studies5, 24, 26–32. ### Local transmission intensity To quantify how well the values of Index P characterize spatial variation in local DENV incidence, we compared the mean Index P during a typical year and the mean yearly incidence in each municipality across Brazil. A typical yearly time series is estimated by calculating the mean value of a chosen variable (e.g. Index P, log(incidence+1), incidence, etc.) independently for each of the twelve months of the year across all available years. The map of transmission potential was consistent with the spatial distribution of incidence (Figure 3a-b). The hyperendemic regions in the Northeast and the Midwest had higher estimated transmission potential, while the low incidence areas in the South stretching from Rio Grande do Sul to Minas Gerais had lower estimated transmission potential. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/11/05/2022.11.04.22281958/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2022/11/05/2022.11.04.22281958/F3) Figure 3. Comparison of reported DENV incidence and mean Index P for a typical year (**a**) Mean yearly incidence in each municipality. (**b**) Mean annual Index P for a typical year. (**c**) The relationship between mean annual Index P and mean yearly incidence. The area below the lower fence (*Q*1 − 1.5 * *IQR* for (0, 0.5) and (0.5, ∞)) which encompasses outliers is shaded grey with municipalities that form part of Brazil’s Legal Amazon highlighted black (these municipalities are also highlighted in panel **a**). Estimated transmission potential and incidence exhibited a nonlinear relationship. Incidence increased with Index P until reaching a plateau in which increases in transmission potential did not yield higher incidence (Figure 3c). This nonlinear relationship is consistent with the interpretation of Index P as the absolute transmission potential of an adult mosquito, where increases in Index P are tied to increases in *R* and the probability and scale of epidemic growth50. The observed plateau at high Index P values can be explained by the dynamics of viral transmission, where the size of epidemics is constrained by factors other than transmission potential, particularly the development of herd immunity over time. As such, municipalities estimated to have very high transmission potential presented similar levels of high incidence. The reported incidence was lower than expected from the estimated transmission potential in several municipalities in the Northwest states of Amazonas and Pará (Figure 3a-b). Given that these areas form part of the Amazon rainforest, these inconsistencies might be due to a combination of poor surveillance, lack of connectivity to allow frequent viral introduction and/or low population density. Indeed, the municipalities within Brazil’s Legal Amazon that had lower than expected incidence had some of the largest areas of all municipalities, and relatively low population densities (Figure 3A, 3c). ### Typical year seasonality We compared the Index P dynamics to the observed DENV seasonality by calculating Spearman’s correlation coefficient between the time series of Index P and log-transformed incidence in each municipality during a typical year (Figure 4). Municipalities with an average of fewer than 12 cases per year (2725/5570) were excluded since their incidence dynamics did not have sufficient information for a measurable seasonal signal. To quantify the time delay between the Index P and incidence dynamics, we calculated the correlation coefficient for all possible monthly lags (i.e. -5 to 6 months inclusive) and selected the lag with the highest coefficient. The lag with the highest correlation coefficient and the coefficient itself are referred to as the lag and lag-adjusted correlation of a municipality respectively. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/11/05/2022.11.04.22281958/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2022/11/05/2022.11.04.22281958/F4) Figure 4. Correlation between the time series of Index P and incidence during a typical year at the municipality level. (**a**) Lag-unadjusted Spearman’s correlation coefficient between Index P and log-transformed incidence during a typical year. The histograms show the distribution of coefficients (mean represented by blue dashed line). Municipalities with fewer than 12 cases per year on average are excluded and coloured grey in the maps (N=2725) given that they have insufficient case data to have detectable seasonal signals. (**b**) Map of the predicted lag of each municipality. The circular barplot shows the distribution of lags from -5 months to +6 months. The proportion of municipalities observed for each lag is noted (proportions less than 1% are not shown). (**c**) Lag-adjusted Spearman’s correlation coefficient between Index P and log-transformed incidence during a typical year. The distribution of lag-unadjusted correlation coefficients was left skewed with 42% of municipalities showing relatively high coefficients greater than 0.6 (Figure 4a). Given the development time of *Aedes* mosquitoes (∼2 weeks), the generation time of DENV in human populations (∼2 −3 weeks)51, 52and the time for mosquitoes to respond to climatic signals, it can take many weeks for increases in transmission potential to translate into increases in incidence5, 53, 54. Indeed, the Index P dynamics in 88% of municipalities preceded the incidence dynamics by 1-3 months, with the proportion of municipalities with time delays of 1, 2 and 3 months being 25%, 47%, and 15% respectively (Figure 4b). When we corrected for this time delay, the dynamics for Index P and incidence were highly correlated with 82% of the municipalities with an average of 12 or more cases per year showing correlation coefficients greater than 0.90 (Figure 4c). Some municipalities (4%) had extreme lags (i.e. negative or greater than 3 months) inconsistent with the biological interpretation of Index P (in the North, spanning the states of Roraima, Pará, Amapá and Maranhão, Figure 4b). Given that these regions have monsoon-driven climates characterized by intense timely annual rainfall, a factor not considered in Index P but known to influence mosquito abundance, precipitation may play an outside role in determining seasonality in these regions. ### Usage Notes Researchers can use Index P maps to explore monthly transmission potential of DENV spread by *Ae. aegypti* mosquitoes at the spatial pixel level (∼28 km2) from 1981 to 2019. There are, however, several caveats to the application of these maps. Firstly, from our technical validation, it can be concluded that areas with low estimated transmission potential can be expected to have lower incidence. However, areas with high estimated potential are susceptible to more uncertainty, as non-climatic factors, particularly accumulated herd immunity and host/vector densities, can influence whether favorable climatic conditions translate into higher disease incidence. High Index P thus should not be interpreted as a guarantee of transmission but rather as one indicator of favorable climatic conditions for DENV circulation. Secondly, while Index P dynamics reliably capture the seasonality of dengue transmission, there is likely to be a time delay between the proposed climatic effects and the observed changes in incidence, the length of which can vary between locations. Our analyses, and those of others53, 54, suggest that this time delay is between about 1 to 3 months, but a comparison of past incidence and Index P dynamics at a local level is necessary for more precise, locally-specific estimates. Finally, Index P is unlikely to mimic incidence patterns at very high spatiotemporal resolutions due to obfuscation of the seasonal signals of incidence by factors such as imperfect surveillance and local stochasticity. Thus Index P is most informative at aggregated spatiotemporal resolutions where incidence is expected to demonstrate clear seasonal signals. To encourage reuse and deepen understanding of the Index P methodology, we developed a new version of the **M**osquito-borne **V**iral **S**uitability **E**stimator (MVSE) software package for the R programming environment24. This R package provides a set of related functions that can be used to estimate Index P time series given temperature and humidity time series and user-defined probability distributions of the biological parameters for the selected host/vector/virus system. A short tutorial on the features of the MVSE package is also made available. Furthermore, we have provided an R Markdown document that highlights useful R software packages and functions for the visualization and analysis of the dataset of Index P maps. Though all code used to generate this dataset was written in R and C++, the maps themselves are provided as TIFF files to facilitate cross-compatibility. ## Data Availability The meteorological data from 1981 to 2019 that were used to estimate the spatiotemporal maps of Index P are freely available from the Essential Climate Variables for assessment of climate variability from 1979 to present dataset published and maintained by Copernicus.eu ([https://cds.climate.copernicus.eu/cdsapp#!/dataset/ecv-for-climate-change](https://cds.climate.copernicus.eu/cdsapp#!/dataset/ecv-for-climate-change)). Weekly notified dengue case data per municipality in Brazil for the years 2000 to 2014 were provided by the Brazilian Ministry of Health ([https://www.gov.br/pt-br](https://www.gov.br/pt-br)) and by the SINAN (Sistema de InformaÇão de Agravos de NotificaÇão - [http://portalsinan.saude.gov.br/](http://portalsinan.saude.gov.br/)). The yearly population size per municipality in Brazil for the years 2001 to 2014 was sourced from IBGE - Instituto Brasileiro de Geografia e Estatística ([https://www.ibge.gov.br/en/home-eng.html](https://www.ibge.gov.br/en/home-eng.html)). The complete dataset of spatiotemporal index P maps and a description of its contents are available from a figshare repository at [https://doi.org/10.6084/m9.figshare.21502614](https://doi.org/10.6084/m9.figshare.21502614). [https://cds.climate.copernicus.eu/cdsapp#!/dataset/ecv-for-climate-change](https://cds.climate.copernicus.eu/cdsapp#!/dataset/ecv-for-climate-change) [https://www.ibge.gov.br/en/statistics/downloads-statistics.html](https://www.ibge.gov.br/en/statistics/downloads-statistics.html) [http://portalsinan.saude.gov.br/](http://portalsinan.saude.gov.br/) [https://doi.org/10.6084/m9.figshare.21502614](https://doi.org/10.6084/m9.figshare.21502614) ## Code availability The MVSE R package can be installed from the GitHub repository at [https://github.com/TaishiNakase/MVSE](https://github.com/TaishiNakase/MVSE). The MVSE package tutorial and the R Markdown document that provides some example code for the visualization of the Index P maps can be found at the GitHub repository [https://github.com/TaishiNakase/Index-P-estimation-and-applications](https://github.com/TaishiNakase/Index-P-estimation-and-applications). ## Author contributions statement T.N. generated and analyzed the data, wrote the code and led the writing of the manuscript. J.L. conceived the idea, supervised the project and contributed extensively to all ideas and writing. M.G. and U.O. read and critically revised the manuscript. ## Competing interests The authors declare no competing interests. * Received November 4, 2022. * Revision received November 4, 2022. * Accepted November 5, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. 1.Franklinos, L. H. V., Jones, K. E., Redding, D. W. & Abubakar, I. The effect of global change on mosquito-borne disease. The Lancet infectious diseases 19, e302–e312 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(19)30161-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F05%2F2022.11.04.22281958.atom) 2. 2.Faria, N. R. et al. Zika virus in the americas: Early epidemiological and genetic findings. Science 352, 345–349 (2016). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTIvNjI4My8zNDUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8xMS8wNS8yMDIyLjExLjA0LjIyMjgxOTU4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 3. 3.Cauchemez, S. et al. Local and regional spread of chikungunya fever in the americas. Euro surveillance 19, 1 (2014). 4. 4.Yactayo, S., Staples, J. E., Millot, V., Cibrelus, L. & Ramon-Pardo, P. Epidemiology of chikungunya in the americas. The J. infectious diseases 214, S441–S445 (2016). 5. 5.Petrone, M. E. et al. Asynchronicity of endemic and emerging mosquito-borne disease outbreaks in the dominican republic. Nat. communications 12, 151–151 (2021). 6. 6.Moore, C. A. et al. Characterizing the pattern of anomalies in congenital zika syndrome for pediatric clinicians. JAMA pediatrics 171, 288–295 (2016). 7. 7.Stegmann-Planchard, S. et al. Chikungunya, a Risk Factor for Guillain-Barré Syndrome. Clin. Infect. Dis. 70, 1233–1235 (2019). 8. 8.Bhatt, S. et al. The global distribution and burden of dengue. Nat. (London) 496, 504–507 (2013). 9. 9.Yang, X., Quam, M. B. M., Zhang, T. & Sang, S. Global burden for dengue and the evolving pattern in the past 30 years. J. Travel. Medicine 28 (2021). 10. 10.Messina, J. P. et al. Global spread of dengue virus types: mapping the 70 year history. Trends microbiology (Regular ed.) 22, 138–146 (2014). 11. 11.European Centre for Disease Prevention and Control (ECDC). Dengue. In: ECDC. Annual epidemiological report for 2019. Available from [https://www.ecdc.europa.eu/en/publications-data/dengue-annual-epidemiological-report-2019](https://www.ecdc.europa.eu/en/publications-data/dengue-annual-epidemiological-report-2019) (2021). 12. 12.Sridhar, S. et al. Effect of dengue serostatus on dengue vaccine safety and efficacy. The New Engl. journal medicine 379, 327–340 (2018). 13. 13.World Health Organisation (WHO). Global strategy for dengue prevention and control 2012-2020. Available from [https://apps.who.int/iris/handle/10665/75303](https://apps.who.int/iris/handle/10665/75303) (2020). 14. 14.Brady, O. J. et al. Modelling adult aedes aegypti and aedes albopictus survival at different temperatures in laboratory and field settings. Parasites & vectors 6, 351–351 (2013). 15. 15.Brady, O. J. et al. Global temperature constraints on aedes aegypti and ae. albopictus persistence and competence for dengue virus transmission. Parasites & vectors 7, 338–338 (2014). 16. 16.Custódio, J. M. d. O. et al. Abiotic factors and population dynamic of aedes aegypti and aedes albopictus in an endemic area of dengue in brazil. Revista do Instituto de Medicina Trop. de Sao Paulo 61, e18–9 (2019). 17. 17.Ramos, M. M. et al. Epidemic dengue and dengue hemorrhagic fever at the texas-mexico border: Results of a household-based seroepidemiologic survey, december 2005. The Am. journal tropical medicine hygiene 78, 364–369 (2008). 18. 18.Kraemer, M. U. G. et al. The global distribution of the arbovirus vectors aedes aegypti and ae. albopictus. eLife 4, e08347–e08347 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.08347&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26126267&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F05%2F2022.11.04.22281958.atom) 19. 19.Hugo, L. E. et al. Adult survivorship of the dengue mosquito aedes aegypti varies seasonally in central vietnam. PLoS neglected tropical diseases 8, e2669–e2669 (2014). 20. 20.Chan, M. & Johansson, M. A. The incubation periods of dengue viruses. PloS one 7, e50972–e50972 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0050972&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23226436&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F05%2F2022.11.04.22281958.atom) 21. 21.Hamlet, A. et al. The seasonal influence of climate and environment on yellow fever transmission across africa. PLoS neglected tropical diseases 12, e0006284–e0006284 (2018). 22. 22.Mordecai, E. A. et al. Detecting the impact of temperature on transmission of zika, dengue, and chikungunya using mechanistic models. PLoS neglected tropical diseases 11, e0005568–e0005568 (2017). 23. 23.Smith, D. L. et al. Ross, macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogens. PLoS pathogens 8, e1002588–e1002588 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1002588&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22496640&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F05%2F2022.11.04.22281958.atom) 24. 24.Obolski, U. et al. Mvse: An r-package that estimates a climate-driven mosquito-borne viral suitability index. Methods ecology evolution 10, 1357–1370 (2019). 25. 25.Messina, J. P. et al. Mapping global environmental suitability for zika virus. eLife 5 (2016). 26. 26.Perez-Guzman, P. N. et al. Measuring mosquito-borne viral suitability in myanmar and implications for local zika virus transmission. PLoS currents 10 (2018). 27. 27.Lourenço, J., Obolski, U., Thompson, R. N. & Thézé, J. Characterising west nile virus epidemiology in israel using a transmission suitability index. Eurosurveillance 25 (2020). 28. 28.Lourenço, J. et al. West nile virus transmission potential in portugal. Commun. biology 5, 6–6 (2022). 29. 29.Costa, E. A. et al. West nile virus in brazil. Pathogens 10 (2021). 30. 30.Lourenço, J., Pinotti, F., Nakase, T., Giovanetti, M. & Obolski, U. Letter to the editor: Atypical weather is associated with the 2022 early start of west nile virus transmission in italy. Eurosurveillance 27 (2022). 31. 31.Carreto, C., Gutiérrez-Romero, R. & Rodríguez, T. Climate-driven mosquito-borne viral suitability index: measuring risk transmission of dengue, chikungunya and zika in mexico. Int. journal health geographics 21 (2022). 32. 32.Adelino, T. E. R. et al. Field and classroom initiatives for portable sequence-based monitoring of dengue virus in brazil. Nat. communications 12, 2296–12 (2021). 33. 33.Ryan, S. J., Carlson, C. J., Mordecai, E. A. & Johnson, L. R. Global expansion and redistribution of aedes-borne virus transmission risk with climate change. PLOS Neglected Trop. Dis. 13, 1–20 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0007012&link_type=DOI) 34. 34.Variza, P. F. et al. Updated spatio-temporal distribution of aedes (stegomyia) albopictus in brazil. Acta tropica 232, 106511–106511 (2022). 35. 35.Gainor, E. M., Harris, E. & Labeaud, A. D. Uncovering the burden of dengue in africa: Considerations on magnitude, misdiagnosis, and ancestry. Viruses 14, 233 (2022). 36. 36.Vasilakis, N., Cardosa, J., Hanley, K. A., Holmes, E. C. & Weaver, S. C. Fever from the forest: prospects for the continued emergence of sylvatic dengue virus and its impact on public health. Nat. reviews. Microbiol. 9, 532–541 (2011). 37. 37.Young, K. I. et al. Abundance and distribution of sylvatic dengue virus vectors in three different land cover types in sarawak, malaysian borneo. Parasites & vectors 10, 406–406 (2017). 38. 38.Lourenço, J. et al. Epidemiological and ecological determinants of zika virus transmission in an urban setting. eLife 6, e29820 (2017). 39. 39.Simmons, C., Farrar, J., Nguyen, V. & Wills, B. Dengue. The New Engl. journal medicine 366 (2012). 40. 40.Yang, H. M., Macoris, M. L. G., Galvani, K. C., Andrighetti, M. T. M. & Wanderley, D. M. V. Assessing the effects of temperature on the population of aedes aegypti, the vector of dengue. Epidemiol. infection 137, 1188–1202 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268809002040&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19192322&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F05%2F2022.11.04.22281958.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268524000014&link_type=ISI) 41. 41.Lambrechts, L. et al. Impact of daily temperature fluctuations on dengue virus transmission by aedes aegypti. Proc. Natl. Acad. Sci. - PNAS 108, 7460–7465 (2011). 42. 42.Essential climate variables for assessment of climate variability from 1979 to present. Copernicus Climate Data Store [https://cds.climate.copernicus.eu/cdsapp#!/dataset/ecv-for-climate-change](https://cds.climate.copernicus.eu/cdsapp#!/dataset/ecv-for-climate-change) (2018). 43. 43.Global climate-driven transmission suitability maps for dengue virus transmitted by aedes aegypti mosquitoes from 1981 to 2019. figshare [https://doi.org/10.6084/m9.figshare.21502614](https://doi.org/10.6084/m9.figshare.21502614) (2022). 44. 44.Sistema de Informação de Agravos de Notificação (SINAN) [http://portalsinan.saude.gov.br/](http://portalsinan.saude.gov.br/) (2016). 45. 45.Instituto Brasileiro de Geografia e Estatística (IBGE) [https://www.ibge.gov.br/en/statistics/downloads-statistics.html](https://www.ibge.gov.br/en/statistics/downloads-statistics.html) (2018). 46. 46.Hijmans, R. J. raster: Geographic Data Analysis and Modeling (2022). R package version 3.6-3. 47. 47.Araujo, V. E. M. d., Bezerra, J. M. T., Amâncio, F. F., Passos, V. M. d. A. & Carneiro, M. Increase in the burden of dengue in brazil and federated units, 2000 and 2015: analysis of the global burden of disease study 2015. Revista brasileira de epidemiologia 20 Suppl 01, 205–216 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1590/1980-5497201700050017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28658384&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F11%2F05%2F2022.11.04.22281958.atom) 48. 48.Peel, M., Finlayson, B. & McMahon, T. Updated world map of the köppen-geiger climate classification. Hydrol. earth system sciences 11, 1633–1644 (2007). 49. 49.Codeco, C. T. et al. Fast expansion of dengue in brazil. The Lancet Reg. Heal. - Am. 12, 100274 (2022). 50. 50.Keeling, M. J. & Rohani, P. Modeling infectious diseases in humans and animals (Princeton University Press, Princeton, N.J., 2008). 51. 51.Aldstadt, J. et al. Space-time analysis of hospitalised dengue patients in rural thailand reveals important temporal intervals in the pattern of dengue virus transmission. Trop. medicine & international health 17, 1076–1085 (2012). 52. 52.Guzzetta, G., Marques-Toledo, C. A., Rosà, R., Teixeira, M. & Merler, S. Quantifying the spatial spread of dengue in a non-endemic brazilian metropolis via transmission chain reconstruction. Nat. communications 9, 2837–8 (2018). 53. 53.Lowe, R. et al. Nonlinear and delayed impacts of climate on dengue risk in barbados: A modelling study. PLoS medicine 15, e1002613 (2018). 54. 54.Yuan, H.-Y. et al. The effects of seasonal climate variability on dengue annual incidence in hong kong: A modelling study. Sci. reports 10, 4297–4297 (2020). 55. 55.Macial-de Freitas, R., Codeço, C. T. & Lourenço-de Oliveira, R. Daily survival rates and dispersal of aedes aegypti females in rio de janeiro, brazil. The Am. journal tropical medicine hygiene 76, 659–665 (2007). 56. 56.Yasuno, M. & Tonn, R. J. A study of biting habits of aedes aegypti in bangkok, thailand. Bull. World Heal. Organ. 43, 319–325 (1970). 57. 57.Trpis, M. & Hausermann, W. Dispersal and other population parameters of aedes aegypti in an african village and their possible significance in epidemiology of vector-borne diseases. The Am. journal tropical medicine hygiene 35, 1263–1279 (1986). 58. 58.Duong, V. et al. Asymptomatic humans transmit dengue virus to mosquitoes. Proc. Natl. Acad. Sci. - PNAS 112, 14688–14693 (2015). 59. 59.Nguyet, M. N. et al. Host and viral features of human dengue cases shape the population of infected and infectious aedes aegypti mosquitoes. Proc. Natl. Acad. Sci. - PNAS 110, 9072–9077 (2013). 60. 60.Nishiura, H. & Halstead, S. B. Natural history of dengue virus (denv)—1 and denv—4 infections: Reanalysis of classic studies. The J. infectious diseases 195, 1007–1013 (2007). [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/inline-graphic-1.gif [4]: /embed/inline-graphic-2.gif [5]: /embed/inline-graphic-3.gif [6]: /embed/graphic-4.gif [7]: /embed/graphic-5.gif [8]: /embed/graphic-6.gif [9]: /embed/graphic-7.gif [10]: /embed/graphic-8.gif [11]: /embed/graphic-9.gif [12]: /embed/graphic-10.gif [13]: /embed/graphic-11.gif [14]: /embed/graphic-12.gif [15]: /embed/graphic-13.gif