Correlation of the basic reproduction number ( R 0) and eco-environmental variables in Colombian municipalities with chikungunya outbreaks during 2014-2016

Abstract

Chikungunya virus (CHIKV) emerged in Colombia in 2014 into a population presumed fully susceptible. This resulted in a quick and intense spread across Colombia, resulting in an epidemic that affected an estimated 450,000 people. The reported Colombian cases accounted for over 49% of all the cases reported to the PAHO. Eco-environmental factors are known to be associated with the spread of arboviruses such as CHIKV, and likely contribute to the differences in transmission profiles that were observed across several municipalities. To determine the association of eco-environmental factors and CHIKV, the basic reproduction number (R0) in 85 municipalities, which accounted for 65.6% of reported CHIKV cases in Colombia, was estimated. Estimates of R0 ranged from 1 to 9, where over 76% of municipalities had R0 values between 1 and 2. When we looked at the distribution of R0, the cumulative proportions were 20% with R0>2, 14% with R0>3, and 9% with R0>4. Next, we determined that there were different patterns of correlation between environmental and/or ecological variables and R0 when we considered different R0 lower-thresholds. Broadly, we found that temperature-related variables are significantly and positively correlated to R0 regardless of the lower threshold, while other variables like duration of outbreak and size of the urban area are inversely related to R0. Specifically, we conclude that high values of temperature-related variables where R0 > 1 will result in a fast growth of cases in a shorter time period (with faster cessation of outbreak transmission) but will result overall in a fewer total cases compared to outbreak areas (R0 > 1, but classified as lower). Thus, in the absence of vector control, a less explosive outbreak may be more advantageous for the virus in terms of transmission.

About the Authors:

Víctor Hugo Peña-García

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (VHPG); [email protected] (RCC)

Affiliation: Programa de Estudio de Control de Enfermedades Tropicales (PECET), Facultad de Medicina, Universidad de Antioquia, Medellín, Colombia

ORCID logo http://orcid.org/0000-0002-4856-3166

Rebecca C. Christofferson

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (VHPG); [email protected] (RCC)

Affiliation: Department of Pathobiological Sciences, School of Veterinary Medicine, Louisiana State University, Baton Rouge, Los Angeles, United States of America

ORCID logo http://orcid.org/0000-0003-2806-1131

Introduction

Chikungunya (CHIKV) is an arbovirus that causes a febrile illness, and belongs to the family Togaviridae and genus Alphavirus. CHIKV is transmitted by vectors belonging to genus Aedes, primarily Ae. aegypti and Ae. albopictus [1]. CHIKV had been limited to Asia and Africa until 2007 when the virus reached Europe, causing sporadic outbreaks, and finally the Americas in 2013, where it quickly caused over 1 million cases [2]. CHIKV emerged in Colombia in 2014, where the population was generally naïve. The subsequent Colombian epidemic of CHIKV accounts for approximately 49.15% of cases in Latin America as reported by the Pan American Health Organization (PAHO) [3].

One of the metrics used to describe the transmission dynamics of infectious diseases is the basic reproduction number (R0), which is defined as the average number of secondary cases generated by a primary case introduced into a fully susceptible population [4, 5]. It has been suggested for vector-borne pathogens that R0 is the expected number of people that will be infected by one person initially infected by a vector [6, 7]. It is theorized that an epidemic occurs when R0 is greater than or equal to 1, although an outbreak is possible even with R0 lower than 1 [7].

The generation time, which is the time between the infection of a primary case and one of its secondary cases, is crucial to estimate of R0 but this value is often unknown [8]. Alternatively, a parameter called the serial interval has been described and is defined as the elapsed time between the onset of the symptoms of the primary and the secondary cases. This interval is more consistently observable in settings where clinical manifestation leads to case detection, but the serial interval is not readily extracted from epidemiological data [8]. Alternatively, some studies of arboviruses have assumed that the serial interval is the sum of estimates of the extrinsic incubation period of a virus (within the mosquito) and the intrinsic incubation periods of a virus (within the vertebrate host) [9, 10]. This is appropriate, as it takes into account the properties of virus transmission both within the vector as well as the vertebrate and is the definition for serial interval we use herein. However, it is important to acknowledge that while the intrinsic incubation period for CHIKV has been estimated to be approximately 3 days, there are several estimates of the extrinsic incubation period, which depend on other extrinsic factors, most notably temperature, leading to a range of possible serial intervals [11–15].

While R0 describes the epidemic potential of a pathogen, it can also be used to investigate modulation of this potential by external factors. For example, temperature is a known effector of mosquito-borne diseases, particularly as it relates to the vector-virus interaction [16–19]. Since arboviral transmission is affected by environmental and ecological variables, we wanted to determine whether any pattern could be discerned from administrative level incidence data in Colombia. Accordingly, we estimated R0 of CHIKV from the Colombian National Institute of Health (INS) and investigated the correlation with many potentially influential eco-environmental variables to determine if and how they influence the CHIKV spread and incidence in an outbreak scenario. Further, given the differences in transmission dynamics across the municipalities, we investigated the discrepancy between the magnitude of R0 and the number of cases.

Methods

Estimation of R0 and variables of interest

R0 was estimated using the White-Pagano methodology, which was shown to be consistent with other methodologies [9, 20]. However, there are advantages to the White-Pagano method when estimates must be made on administrative incidence data. An important advantage of White-Pagano method is that serial interval can be directly approximated from epidemiological data. Given that the White-Pagano methodology allows for the estimation of R0 without a priori observation of the exact serial interval, it is suitable for the estimation of R0 for different municipalities with different arboviruses transmission conditions, such as temperature profiles, where the serial interval likely varies. However, the WP method does require an estimated range of possible serial intervals. Herein we assumed a possible serial interval of 1 or 2 weeks so that the unit of time matched the incidence data (as we do not know the distribution of cases in days within each epidemiological week), which allows for some range in estimation of CHIKV serial interval that could be related to variability in EIP [21–24]. The likelihood-based method used for the estimation of the R0 and the serial interval (time between the onset of symptoms of first case and the onset of symptoms of the secondary case [8]) is as follows:where , Nt is the incidence data; Pj is the serial interval with j = 1, …, k; and t is the time of epidemiological data. Using this methodology, it is possible to estimate both R0 and Pj by using maximum likelihood techniques [20]. Thus, the White-Pagano methodology allows for estimation of the R0 even when the serial interval is not known a priori, which makes this methodology suitable for the estimation of R0 for arboviruses in different municipalities, i.e. with different temperatures which can modify the length of the extrinsic incubation period within the mosquito vector. However, this methodology is sensitive to the total number of cases in the system. Thus, we selected the Colombian municipalities with the highest number of cases in order to confidently apply this derivation of R0 to comparisons across environmental variables of the municipalities [20].

With the range of serial intervals used to initiate the model of 1–2 weeks, we assumed a binomial distribution for the serial interval for all municipalities [20]. The analysis determined the probability of having a 1- or 2-week serial interval for each municipality and assigned the serial interval with the highest probability (S1 Fig). Additionally, given the 2-week maximum serial interval, case counts were included in the determination of the epidemic curve per municipality if new cases occurred within two weeks of the previous case. If a case occurs outside of this two-week window, it was considered sporadic or possibly the results of a new introduction and not included. Having shaped the epidemic curve for each local epidemic, the estimation of R0 included cases starting from the week of the first case until the local peak of the epidemic. The outbreak decline phase of the epidemic is not taken into account for the R0 estimation using this methodology [20].

Epidemiological and demographic data

Epidemiological data was obtained from Colombian INS which publishes through its web page (www.ins.gov.co) the reported cases collected by the national system of epidemiological surveillance (SIVIGILA). Data for CHIKV cases was taken for the years 2014 to 2016. Data collected through SIVIGILA is reported on a weekly basis, i.e. weekly number of chikungunya cases per city reporting. Population projections developed by the Colombian National administrative department of Statistics (DANE according to its name in Spanish) were used to estimate population size per city and subsequently incidence. Incidence was estimated per 100,000 individuals per each city. Urban area size (termed “urbanicity”) and map resources were obtained from the National Geostatistical Framework (MGN, according to its name in Spanish) of DANE.

Environmental data

Temperature data of each city were extracted from Worldclim climatic layers [25]. We used data of minimum temperature, maximum temperature and average temperature by using the coordinates taken from the center of the urban area of each municipality that were obtained from the software Google Earth Pro version 7.3.2.5491 (Google LLC, Mountain view, California, USA). The difference of temperature (maximum temperature minus minimum temperature) was estimated from that data.

For each municipality, we calculated the number of cumulative cases until epidemic peak (“Peak cCases”), the time to onset of detected outbreak in the municipality (the first week with reported cases in the municipality relative to the first week with reported cases in Colombia, “outbreak onset”), and duration of the outbreak (in weeks).

Statistical analysis

All analyses were performed in R version 3.5.1. After R0 was estimated as above, we correlated the values of environmental and demographic variables with the R0 estimate using Spearman’s correlation across all municipalities. Correlations with p<0.05 were considered significant. Following this, we wanted to determine if there was a differential pattern of corollaries dependent on the intensity of transmission. As R0 is a threshold parameter (lower bound of 1 to indicate outbreak likelihood), we did this by altering the (lower) threshold to indicate differential intensity of transmission. Thus, we conducted the same Spearman’s correlation analysis for the following subsets of municipalities: R0 > 2, R0 > 2.5, and R0 > 3. These were chosen based on the distribution of R0 (see Results below).

We then wanted to determine if the variable values (such as temperature) were differentially associated depending on the R0 value. For this, we divided municipalities into three pairwise comparisons: 1 < R0 ≤ 2 versus 2 < R0; 1 < R0 < 2.5 versus 2.5 < R0; and 1 < R0 < 3 versus 3 < R0. Comparisons were done using Wilcoxon rank sum test.

Maps were created in ArcGIS 10.4.1 from publicly available dataset described above and from our calculation of R0 from those datasets.

Results

Estimation of R0 from CHIKV incidence data

From 846 municipalities that reported chikungunya cases between 2014 and 2016 we selected 85 that had a recognizable epidemic peak and reported more than 150 cases. The chosen municipalities comprised approximately 65.6% of Colombian reported cases (Fig 1). From those municipalities, the range of cases reported over the two years of peak CHIKV incidence in Colombia was 151 (city of Soledad) to 4179 cases (city of Cali). The city of Turbaco reported the highest number of cases in a single week, which was 589 cases.

[Figure omitted. See PDF.]

Fig 1. General behavior of the chikungunya epidemic that occurred in Colombia from 2014 to 2016.

(A) Cumulative incidence discriminated by municipalities in Colombia according to cases reported through SIVIGILA by INS, (B) R0 estimated for 85 municipalities with epidemic CHIKV in Colombia and their geographical locations, (C) Behavior in time of the total reported cases in Colombia (blue line) and the reported cases of the 85 municipalities chosen to estimate R0 (red line). All maps were generated using ArcGIS 10.4.1 for this publication by the authors and has not been published elsewhere.

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

The R0 in the 85 municipalities ranked from 1.11 in Floridablanca (Department of Santander) to 9.53 in the municipality of El Molino (La Guajira department). Most of the municipalities were estimated to have a one-week serial interval rather than a two-week serial interval (S1 Fig). The distribution of R0 values estimated are given in S2 Fig. Estimates of R0 ranged from 1 to 9, where over 76% of municipalities had R0 values between 1 and 2. When we looked at the distribution of R0, the cumulative proportions were 20% with R0>2, 14% with R0>3, and 9% with R0>4.

When we visualized the estimated R0 versus CHIKV incidence in each municipality, we found the relationship to be relatively dichotomous with an apparent break at approximately R0 > 2 (Fig 2). While the relationship between incidence and R0 is obvious, the stark separation shown by this relationship was surprising and further lead us to investigate alternate lower bounds of 2.5 and 3 (see below).

[Figure omitted. See PDF.]

Fig 2. Relationship of R0 and cumulative CHIKV incidence of 85 municipalities during their respective outbreaks between 2014 and 2016.

R0 higher than 2 are highlighted in red while blue indicates those municipalities with 1 < R0 < 2.

https://doi.org/10.1371/journal.pntd.0007878.g002

Correlation of R0 and environmental/ecological variables when lower bound or R0 is varied

Table 1 outlines the results of the correlation analysis of 13 variables. When all municipalities were considered, the number of cases, mean, max, and min temperature, as well as the duration of the epidemic and size of urban area were all significantly correlated to R0. Similar to what is seen in Fig 2, there was a surprising negative correlation between R0 and the duration and number of cases. Temperature was positively associated with R0 values.

[Figure omitted. See PDF.]

Table 1. Spearman correlation coefficients (ρ) between R0 and different variables when the lower threshold of R0 is 2, 2.5, or 3.

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

We determined that the magnitude and direction of the correlation coefficients of some variables changed according to the lower bound of R0 (Table 2). When the lower bound of R0 was greater than 2, the maximum temperature, total population, urban population, and incidence were significant. Number of cases and duration of outbreak was not significantly correlated with the value of R0, indicating that the relationship between the shape and size of the epidemic and R0 is driven mostly by those outbreaks where 1≤R0≤2. As the lower bound was further raised to 2.5 and 3, correlation with these variables was not statistically significant. The only exception was a statistically significant and negative correlation between R0 and the minimum temperature when the lower bound was 3, indicating that the biological trend of higher temperatures leading to more transmission is still apparent, but not captured with maximum temperature solely.

[Figure omitted. See PDF.]

Table 2. P-value of Wilcoxon sum rank test with different variables between municipalities separated according to groups where R0 > of 2, 2.5 or 3 compared to groups where R0 ≤ 2, 2.5, or 3.

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

Association of R0 when lower bound is modified with eco-environmental variables

The Wilcoxon rank sum test showed that there were differences in several variables when compared to groups of “low” versus “high” based on the lower bound of the R0 estimates (Table 2). We found consistent differences in the distribution of number of cases when all municipalities were considered (lower threshold R0 = 1) versus the modified lower thresholds of 2 or 2.5. Mean, minimum, and difference in temperature was not consistently different across all comparisons. Mean and minimum temperature were only significant between the comparison of R0>1 versus R0>3; difference in temperature was only significant in the comparisons of R0>1 versus R0>2.5. Maximum temperature was different among all comparisons, indicating that there is a consistent relationship between this variable and increasing R0 estimates.

Population size was not statistically associated with differences in R0 estimates, and neither was the number of people in the urban center, incidence during the epidemic, population density, or the outbreak onset. Interestingly, both the duration of the epidemic and the size of the urban area were statistically significant in all comparisons (Figs 3 and 4).

[Figure omitted. See PDF.]

Fig 3. Distribution of values of variables related to temperature of the municipalities with significant differences when they are separated by a R0 value threshold of 2, 2.5 or 3.

https://doi.org/10.1371/journal.pntd.0007878.g003

[Figure omitted. See PDF.]

Fig 4. Distribution of values of different variables of municipalities with significant differences when they are separated by a R0 value threshold of 2, 2.5 or 3.

https://doi.org/10.1371/journal.pntd.0007878.g004

Maximum temperature was positively correlated and consistently differential among comparisons. Fig 3 shows this relationship where those municipalities in the lower group had a tendency towards lower temperatures than those municipalities in the higher R0. Interestingly, the total number of cases and size of the urban area were negatively correlated and consistently statistically different among lower threshold R0 groupings. Surprisingly, lower R0 groupings had a tendency to have higher incidence of CHIKV. The size of the urban area was also negatively and consistently different, with higher R0 groupings having smaller urban areas. Finally, the duration of the epidemic tended to be shorter in those municipalities with higher R0 estimates (Fig 4).

Discussion

The White-Pagano methodology allows for estimation of both the serial interval and R0 using administrative level, epidemiological data [20]. This method is ideal for studies that rely on administrative level data, as it does not require a priori estimation of parameters related to mosquito populations, which may be needed on a per-municipality basis [26]. Using this method, we estimated that in municipalities that account for over 64% of the CHIKV cases in Colombia, the R0 ranged from 1.11 to over 9. This wide range indicates heterogeneity in transmission that, since Colombia was assumed to be a completely susceptible population, could have been related to other factors than the lack of herd immunity.

We investigated several variables that have been associated with arbovirus transmission (and thus R0) including temperature, urbanicity, and epidemiological characteristics [13, 27–30]. We confirmed what other studies have demonstrated, that higher temperatures (as measured by maximum temperature) is associated with several transmission factors that could lead to higher values of R0 [12, 18, 27, 31–37]. Though we consider a spectrum of environmental, ecological, and social variables, it is important to note that many other variables likely affect transmission and reporting of cases at different spatial resolutions. Variables such as housing type/quality, income levels, water availability, or access to health services likely all play a role in risk of exposure to arboviruses [38, 39]. However, mosquito density and other socio-ecological variables vary across and within municipalities and are not readily available at the country-level and thus the inclusion of these specific types of variables is of limited value to our model, but urban area size and density of population, included in our study, have been used as proxies for some of these other factors, such as mosquito population [38].

The finding that higher R0 is associated with smaller urban areas is interesting, especially given that there was a lack of consistently significant relationship between R0 and the total population of the urban area, the total municipality population, and the density of population indicated. Previous work had revealed that transmission by Aedes aegypti is positively associated with intermediate levels of human movement [40]. Particularly with CHIKV, there is evidence in Dhaka, Bangladesh that viral spread is mainly driven by transmission happening at distances that do not exceed neighborhood distances [41]. Our results indicate that higher R0 values are associated with smaller urban areas, which could be indicative of limited movement of humans in that area. The result is highly focal transmission, which has been reported for CHIKV [42].

Of interest was the negative correlation of R0 grouping with urbanicity, CHIKV case counts, and outbreak time-to-peak. That is, those municipalities with higher values of R0 showed an explosive outbreak with a fast increase of the number of CHIKV cases but a shorter time-to-peak. These higher R0 municipalities also had overall lower case counts at the time of peak compared to municipalities that displayed lower values of R0, where the epidemic continued for many weeks with smaller peak values, but ultimately higher case counts. This phenomenon can be seen in Fig 5, where we plotted the course of epidemic of the three municipalities with highest values of R0 and the three municipalities with the lowest values of R0 (but still ≥1).

[Figure omitted. See PDF.]

Fig 5. Epidemic peak of three municipalities with the highest (upper) and the three municipalities with lowest (below) R0 estimated.

Above each plot, the cumulative case count from the start until the highest outbreak peak is given.

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

In summary, this work suggests that CHIKV outbreaks are affected by high maximum temperatures and smaller urban areas, which leads to higher estimates of R0. These higher R0 values are then associated with explosive outbreaks with shorter times-to-peak and higher case counts at time of peak (Fig 5). This explosive nature of the initial outbreak leads to a lower overall case count compared to those outbreaks that had longer times-to-peak, lower peak-case counts, but cumulatively larger total number of cases. The tortoise-or-the-hare (TotH) model was proposed to also explain differences in infectious dynamics of arboviruses whereby arboviruses that adopt the tortoise strategy (long-lasting, moderate viremia) will persist longer than arboviruses that adopt the hare strategy (short-lived, high viremia) [43]. We show here that this TotH model, previously used to describe micro-dynamics, may also be suited to describe macro-level dynamics of arboviruses. In theory, the tortoise model would seem to be most advantageous for the virus, as it ultimately achieves the longest and most productive transmission. However, in the context of active surveillance and response where vector control may be activated and result in abrupt cessation or decrease in transmission, the hare model may ultimately represent the most fit population-level phenotype [40, 41]. That is, in areas where control is implemented–such as in Brazil during the Zika epidemic, or Miami during local transmission there–a faster “out-of-the-gate” phenotype (the hare) may result in more transmission, which could be interpreted as the most fit phenotype in this scenario [44, 45]. The TotH model of how to think about transmission represents an interesting relationship between perceived transmission intensity and the ultimate magnitude of transmission (i.e. number of cases). Thus, these findings could impact response prioritization, especially in areas where resources are limited. That is, more intense transmission may appear to be the most pressing concern, when in reality the less intense situations may need more intervention… And given the association of R0 with eco-environmental variables, the TotH relationship may further have use in predicting allocation of resources at a focal level.” Our findings suggest that the non-monotonic relationship between arbovirus transmission and temperature [46] can be at least in part detected with administrative level case counts and the relationships to the lower bounds of R0. And given the association of R0 with eco-environmental variables, the TotH relationship may further have use in predicting allocation of resources at a focal level.

Supporting information

S1 Fig. Frequency of probability values of the 85 municipalities of having one- and two-weeks serial interval.

https://doi.org/10.1371/journal.pntd.0007878.s001

(JPG)

S2 Fig. Frequency of CHIKV R0 values across 85 Colombian municipalities during epidemic peak in 2014–2016.

https://doi.org/10.1371/journal.pntd.0007878.s002

(JPG)

Acknowledgments

We would like to specially thank to Dr. Laura White for kindly providing R scripts to estimate R0 and to Dr. Michael Robert for his valuable feedback.

Citation: Peña-García VH, Christofferson RC (2019) Correlation of the basic reproduction number (R0) and eco-environmental variables in Colombian municipalities with chikungunya outbreaks during 2014-2016. PLoS Negl Trop Dis 13(11): e0007878. https://doi.org/10.1371/journal.pntd.0007878

References

1. Horwood PF, Buchy P. Chikungunya. Rev Sci Tech. 2015;34(2):479–89. pmid:26601450.

2. Weaver SC, Forrester NL. Chikungunya: Evolutionary history and recent epidemic spread. Antiviral Res. 2015;120:32–9. Epub 2015/05/12. pmid:25979669.

3. Oviedo-Pastrana M, Méndez N, Mattar S, Arrieta G, Gomezcaceres L. Lessons learned of emerging Chikungunya virus in two populations of social vulnerability of the Colombian tropics: epidemiological analysis. Arch Public Health. 2018;76:36. Epub 2018/07/23. pmid:30062010.

4. Diekmann O, Heesterbeek JA, Metz JA. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J Math Biol. 1990;28(4):365–82. pmid:2117040.

5. Heesterbeek JAP, Dietz K. The concept of Ro in epidemic theory. Stat Neerl. 1996;50(1):89–110.

6. Lopez LF, Coutinho FA, Burattini MN, Massad E. Threshold conditions for infection persistence in complex host-vectors interactions. C R Biol. 2002;325(11):1073–84. pmid:12506721.

7. Massad E, Coutinho FA, Burattini MN, Amaku M. Estimation of R0 from the initial phase of an outbreak of a vector-borne infection. Trop Med Int Health. 2010;15(1):120–6. Epub 2009/11/03. pmid:19891761.

8. Svensson A. A note on generation times in epidemic models. Math Biosci. 2007;208(1):300–11. Epub 2006/11/09. pmid:17174352.

9. Sanches RP, Massad E. A comparative analysis of three different methods for the estimation of the basic reproduction number of dengue. Infect Dis Model. 2016;1(1):88–100. Epub 2016/08/31. pmid:29928723.

10. Aldstadt J, Yoon IK, Tannitisupawong D, Jarman RG, Thomas SJ, Gibbons RV, et al. Space-time analysis of hospitalised dengue patients in rural Thailand reveals important temporal intervals in the pattern of dengue virus transmission. Trop Med Int Health. 2012;17(9):1076–85. Epub 2012/07/19. pmid:22808917.

11. Rudolph KE, Lessler J, Moloney RM, Kmush B, Cummings DA. Incubation periods of mosquito-borne viral infections: a systematic review. Am J Trop Med Hyg. 2014;90(5):882–91. Epub 2014/03/17. pmid:24639305.

12. Chan M, Johansson MA. The incubation periods of Dengue viruses. PLoS One. 2012;7(11):e50972. pmid:23226436

13. Alto BW, Wiggins K, Eastmond B, Ortiz S, Zirbel K, Lounibos LP. Diurnal Temperature Range and Chikungunya Virus Infection in Invasive Mosquito Vectors. J Med Entomol. 2018;55(1):217–24. Epub 2017/10/19. pmid:29040730.

14. Gubler DJ, Reiter P, Ebi KL, Yap W, Nasci R, Patz JA. Climate variability and change in the United States: potential impacts on vector- and rodent-borne diseases. Environ Health Perspect. 2001;109 Suppl 2:223–33. pmid:11359689.

15. Kakarla SG, Mopuri R, Mutheneni SR, Bhimala KR, Kumaraswamy S, Kadiri MR, et al. Temperature dependent transmission potential model for chikungunya in India. Sci Total Environ. 2018;647:66–74. Epub 2018/08/01. pmid:30077856.

16. Tesla B, Demakovsky LR, Mordecai EA, Ryan SJ, Bonds MH, Ngonghala CN, et al. Temperature drives Zika virus transmission: evidence from empirical and mathematical models. Proc Biol Sci. 2018;285(1884). Epub 2018/08/17. pmid:30111605.

17. Heitmann A, Jansen S, Luhken R, Leggewie M, Badusche M, Pluskota B, et al. Experimental transmission of Zika virus by mosquitoes from central Europe. Euro Surveill. 2017;22(2). pmid:28106528.

18. Carrington LB, Armijos MV, Lambrechts L, Scott TW. Fluctuations at a Low Mean Temperature Accelerate Dengue Virus Transmission by Aedes aegypti. PLoS Negl Trop Dis. 2013;7(4):e2190. pmid:23638208.

19. Liu Z, Zhang Z, Lai Z, Zhou T, Jia Z, Gu J, et al. Temperature Increase Enhances Aedes albopictus Competence to Transmit Dengue Virus. Front Microbiol. 2017;8:2337. pmid:29250045.

20. White LF, Pagano M. A likelihood-based method for real-time estimation of the serial interval and reproductive number of an epidemic. Stat Med. 2008;27(16):2999–3016. pmid:18058829.

21. Mbaika S, Lutomiah J, Chepkorir E, Mulwa F, Khayeka-Wandabwa C, Tigoi C, et al. Vector competence of Aedes aegypti in transmitting Chikungunya virus: effects and implications of extrinsic incubation temperature on dissemination and infection rates. Virol J. 2016;13:114. Epub 2016/07/01. pmid:27357190.

22. Alto BW, Wiggins K, Eastmond B, Velez D, Lounibos LP, Lord CC. Transmission risk of two chikungunya lineages by invasive mosquito vectors from Florida and the Dominican Republic. PLoS Negl Trop Dis. 2017;11(7):e0005724. Epub 2017/07/28. pmid:28749964.

23. Wimalasiri-Yapa MCRB, Stassen L, Hu W, Yakob L, McGraw AE, Pyke TA, et al. Chikungunya Virus Transmission at Low Temperature by Aedes albopictus Mosquitoes. Pathogens. 2019;8(3). pmid:31547257

24. Kakarla SG, Mopuri R, Mutheneni SR, Bhimala KR, Kumaraswamy S, Kadiri MR, et al. Temperature dependent transmission potential model for chikungunya in India. Sci Total Environ. 2019;647:66–74. Epub 2018/08/06. pmid:30077856.

25. Fick SE, Hijmans RJ. WorldClim 2: new 1‐km spatial resolution climate surfaces for global land areas. Int J Climatol. 2017;37(12):4302–15.

26. Peña-García VH, Sánchez-Vargas I, Christofferson R, Black WC Iv, Arboleda S, Triana-Chavez O. Estimation of DENV-2 Transmission as a Function of Site-Specific Entomological Parameters from Three Cities in Colombia. Ann Glob Health. 2019;85(1). Epub 2019/03/13. pmid:30873777.

27. Mordecai EA, Cohen JM, Evans MV, Gudapati P, Johnson LR, Lippi CA, et al. Detecting the impact of temperature on transmission of Zika, dengue, and chikungunya using mechanistic models. PLoS Negl Trop Dis. 2017;11(4):e0005568. Epub 2017/04/27. pmid:28448507.

28. Overgaard HJ, Olano VA, Jaramillo JF, Matiz MI, Sarmiento D, Stenström TA, et al. A cross-sectional survey of Aedes aegypti immature abundance in urban and rural household containers in central Colombia. Parasit Vectors. 2017;10(1):356. Epub 2017/07/27. pmid:28750651.

29. Ndenga BA, Mutuku FM, Ngugi HN, Mbakaya JO, Aswani P, Musunzaji PS, et al. Characteristics of Aedes aegypti adult mosquitoes in rural and urban areas of western and coastal Kenya. PLoS One. 2017;12(12):e0189971. Epub 2017/12/19. pmid:29261766.

30. Fuller TL, Calvet G, Genaro Estevam C, Rafael Angelo J, Abiodun GJ, Halai UA, et al. Behavioral, climatic, and environmental risk factors for Zika and Chikungunya virus infections in Rio de Janeiro, Brazil, 2015–16. PLoS One. 2017;12(11):e0188002. Epub 2017/11/16. pmid:29145452.

31. Yang HM, Macoris ML, Galvani KC, Andrighetti MT, Wanderley DM. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiol Infect. 2009;137(8):1188–202. pmid:19192322.

32. Grech MG, Sartor PD, Almirón WR, Ludueña-Almeida FF. Effect of temperature on life history traits during immature development of Aedes aegypti and Culex quinquefasciatus (Diptera: Culicidae) from Córdoba city, Argentina. Acta Trop. 2015;146:1–6. Epub 2015/02/27. pmid:25733491.

33. Watts DM, Burke DS, Harrison BA, Whitmire RE, Nisalak A. Effect of temperature on the vector efficiency of Aedes aegypti for dengue 2 virus. Am J Trop Med Hyg. 1987;36(1):143–52. pmid:3812879.

34. Lambrechts L, Paaijmans KP, Fansiri T, Carrington LB, Kramer LD, Thomas MB, et al. Impact of daily temperature fluctuations on dengue virus transmission by Aedes aegypti. Proc Natl Acad Sci U S A. 2011;108(18):7460–5. pmid:21502510.

35. Carrington LB, Seifert SN, Armijos MV, Lambrechts L, Scott TW. Reduction of Aedes aegypti vector competence for dengue virus under large temperature fluctuations. Am J Trop Med Hyg. 2013;88(4):689–97. pmid:23438766.

36. Huang YJ, Hsu WW, Higgs S, Vanlandingham DL. Temperature Tolerance and Inactivation of Chikungunya Virus. Vector Borne Zoonotic Dis. 2015;15(11):674–7. pmid:26565772.

37. Park SL, Huang YJ, Hsu WW, Hettenbach SM, Higgs S, Vanlandingham DL. Virus-specific thermostability and heat inactivation profiles of alphaviruses. J Virol Methods. 2016;234:152–5. Epub 2016/04/11. pmid:27079828.

38. Piedrahita LD, Agudelo Salas IY, Marin K, Trujillo AI, Osorio JE, Arboleda-Sanchez SO, et al. Risk Factors Associated with Dengue Transmission and Spatial Distribution of High Seroprevalence in Schoolchildren from the Urban Area of Medellin, Colombia. Can J Infect Dis Med Microbiol. 2018;2018:2308095. Epub 2018/09/25. pmid:30245759.

39. Delmelle E, Hagenlocher M, Kienberger S, Casas I. A spatial model of socioeconomic and environmental determinants of dengue fever in Cali, Colombia. Acta Trop. 2016;164:169–76. Epub 2016/10/30. pmid:27619189.

40. Stone CM, Schwab SR, Fonseca DM, Fefferman NH. Human movement, cooperation and the effectiveness of coordinated vector control strategies. J R Soc Interface. 2017;14(133). pmid:28855386.

41. Salje H, Lessler J, Paul KK, Azman AS, Rahman MW, Rahman M, et al. How social structures, space, and behaviors shape the spread of infectious diseases using chikungunya as a case study. Proc Natl Acad Sci U S A. 2016;113(47):13420–5. Epub 2016/11/07. pmid:27821727.

42. Salje H, Cauchemez S, Alera MT, Rodriguez-Barraquer I, Thaisomboonsuk B, Srikiatkhachorn A, et al. Reconstruction of 60 Years of Chikungunya Epidemiology in the Philippines Demonstrates Episodic and Focal Transmission. J Infect Dis. 2016;213(4):604–10. Epub 2015/09/25. pmid:26410592.

43. Althouse BM, Hanley KA. The tortoise or the hare? Impacts of within-host dynamics on transmission success of arthropod-borne viruses. Philos Trans R Soc Lond B Biol Sci. 2015;370(1675). pmid:26150665.

44. Heukelbach J, Alencar CH, Kelvin AA, de Oliveira WK, Pamplona de Goes Cavalcanti L. Zika virus outbreak in Brazil. J Infect Dev Ctries. 2016;10(2):116–20. Epub 2016/03/02. pmid:26927450.

45. Philip C, Novick CG, Novick LF. Local Transmission of Zika Virus in Miami-Dade County: The Florida Department of Health Rises to the Challenge. J Public Health Manag Pract. 2019;25(3):277–87. Epub 2019/04/02. pmid:30933006.

46. Mordecai EA, Caldwell JM, Grossman MK, Lippi CA, Johnson LR, Neira M, et al. Thermal biology of mosquito-borne disease. Ecol Lett. 2019. Epub 2019/07/10. pmid:31286630.