- Open Access
Bionomic response of Aedes aegypti to two future climate change scenarios in far north Queensland, Australia: implications for dengue outbreaks
© Williams et al.; licensee BioMed Central Ltd. 2014
- Received: 10 June 2014
- Accepted: 16 September 2014
- Published: 19 September 2014
Dengue viruses are transmitted by anthropophilic mosquitoes and infect approximately 50 million humans annually. To investigate impacts of future climate change on dengue virus transmission, we investigated bionomics of the mosquito vector, Aedes aegypti.
Using a dynamic life table simulation model (the Container inhabiting mosquito simulation CIMSiM) and statistically downscaled daily values for future climate, we assessed climate change induced changes to mosquito bionomics. Simulations of Ae. aegypti populations for current (1991-2011) and future climate (2046-2065) were conducted for the city of Cairns, Queensland, the population centre with most dengue virus transmission in Australia. Female mosquito abundance, wet weight, and the extrinsic incubation period for dengue virus in these mosquitoes were estimated for current and future climate (MPI ECHAM 5 model, B1 and A2 emission scenarios).
Overall mosquito abundance is predicted to change, but results were equivocal for different climate change scenarios. Aedes aegypti abundance is predicted to increase under the B1, but decrease under the A2 scenario. Mosquitoes are predicted to have a smaller body mass in a future climate. Shorter extrinsic incubation periods are projected.
It is therefore unclear whether dengue risk would increase or decrease in tropical Australia with climate change. Our findings challenge the prevailing view that a future, warmer climate will lead to larger mosquito populations and a definite increase in dengue transmission. Whilst general predictions can be made about future mosquito borne disease incidence, cautious interpretation is necessary due to interaction between local environment, human behaviour and built environment, dengue virus, and vectors.
- Future Climate
- Dengue Virus
- Future Climate Scenario
- Vectorial Capacity
- Mosquito Abundance
Two point five billion people globally are at risk of dengue virus (DENV) infection, transmitted by Aedes mosquitoes mostly in urban and peri-urban areas. Approximately 50 million DENV infections occur annually, with a disease burden of around 1300 disability-adjusted life years .
The principal DENV vectors are Aedes aegypti (L.) and Aedes albopictus (Skuse); mosquitoes that utilize water-filled structures in human settlements for breeding. The contribution of rainfall, humidity and temperature to the provision of mosquito habitat and the growth and development of these ectotherms means there is a strong link between climate and dengue disease, and changes to climate are projected to alter future transmission. A rise in global mean surface temperature of 0.4-2.6°C has been forecast for the mid-21st century, with spatially heterogeneous changes to precipitation and sea level . Thus, the intensity and distribution of mosquito-borne diseases, including dengue, is expected to alter in a future climate.
A number of studies predicting changes to future disease range and intensity exist, and there is a body of evidence that suggests increased geographic range for dengue virus transmission in a future climate. Some of these are based on the relationships between climate variables and disease incidence. There are models that describe a positive relationship between notified dengue cases and increasing temperature  and humidity . Based on these models projected climate change has been predicted to increase dengue transmission by expansion of geographic range [3–5] and increased transmission intensity and seasonal duration in areas already at risk [3, 6]. Generally a global increase in the ability of mosquitoes to transmit dengue viruses has been predicted. But more recently the impact of social factors such as vector control, housing quality and income has been predicted to mitigate dengue risk. Use of a statistical model incorporating gross domestic product demonstrated that climate change effects alone would increase the probability of dengue transmission, but increasing gross domestic product will decrease future risk .
The geographic range of dengue has generally decreased in recent decades , particularly in developed nations such as Australia [8, 9], despite a warming climate since the mid 20th century . Recent warming climate has not been associated strongly with increased dengue transmission .
Superficially observations and projections are discordant. This may be attributable in part to the simplicity of dengue models. Statistical models that relate climate variables to disease presence and incidence over large geographic areas may produce projections that are in error, as they are based on historic associations between weather and disease. Conversely, models based on observations from smaller spatial scales have the facility to incorporate locally-relevant aspects of dengue ecology, as do models that incorporate mechanistic processes determining the ecology of the vector . Mechanistic modeling with the dengue simulation model (DENSiM) has been used to model inter-annual variability in dengue activity in Queensland , demonstrating that local meteorological phenomena impact dengue transmission.
Relationships between temperature and mosquito fitness have been established for a number of mosquito species , and thermal optima (ranges of temperatures at which fitness is maximized) have been established. In general, climate warming is expected to favour mosquitoes that have broad thermal optima and can adapt to higher temperatures . Studies of reproduction of Ae. aegypti indicate that it will adapt well to moderate increases in temperature .
The relationship between climate and DENV vector mosquito productivity and bionomics has not been investigated. There have been predictions that climate change will increase the abundance of disease vectors because mosquitoes will develop faster and complete more generations annually, and that warmer temperatures will decrease virus incubation time , with both phenomena contributing to increased dengue incidence. However, these predictions have not been supported by quantitative evidence from mosquito population dynamics models.
As part of a larger study to understand the impact of future climate change on dengue transmission in Australia, we investigated likely changes to Ae. aegypti bionomics under current and future climate. Climate has changed through time and will continue to do so. Further, changes to society, economy, housing, and travel patterns will also go on. Therefore, the causes for change in dengue epidemiology are many, and disentangling the predominant risk determinants is challenging.
In an effort to do this, we investigated how mosquito abundance, body size and virus incubation rate are projected to alter utilizing a field-validated mechanistic mosquito population model and climate projections. This was done for an area of annual dengue transmission and high receptivity to dengue importation, Cairns, in Far North Queensland, Australia. This location is the hub of local dengue transmission in Australia, with the largest epidemic occurring there in 2008-2009 (over 900 cases and one death) . Given the history of dengue vector and epidemiology studies performed there, Cairns presents an ideal location to investigate potential impacts of climate change on dengue vector bionomics.
The Container-inhabiting Mosquito Simulation (CIMSiM)
The modelling software CIMSiM 3.27 (University of California) was used to simulate Ae. aegypti population dynamics. CIMSiM  accurately models Ae. aegypti population dynamics in Queensland , and has been used to study Ae. aegypti productivity [9, 17] and persistence . CIMSiM was also used to provide entomological data for modeling inter-annual variation in dengue transmission in Australia  using the DENSiM (dengue simulation) model attached to it . CIMSiM generates daily estimates of egg, larval, pupal and adult numbers per hectare by integrating daily meteorological observations with information about available breeding habitats. In this way, CIMSiM can be applied to any locality in which the container type and frequency are known and for which meteorological data are available.
Breeding container parameters as used in Aedes aegypti population modelling using CIMSiM
Pot plant saucers
18.0 d, 17.4 h
6.1 × 2.0 × 26.0
15.5 × 22.6 × 4.2
35.3 d, 11.0 w
64.6 × 72.7 × 172.0
200 d × 200 h
Water shed ratio
Draw down (L)
Initial food (mg)
density per ha.
Meteorological data for current and future climates
Carbon cycle models generally give estimates of atmospheric carbon dioxide concentrations ranging from 500 to 1200 ppm for the mid 21st century compared to concentrations of 352 ppm in 1990 . The two emission scenarios A2 and B1 differ in their estimated atmospheric carbon dioxide concentration and consequent temperature rises. The A2 emission scenario is estimated to lead to an atmospheric carbon dioxide concentration between 372 and 572 ppm. This scenario assumes a relatively slower transition to renewable energy sources with consequent high carbon emissions, and likely global temperature increases from 2.0-5.4°C by the end of the 21st C . In comparison the B1 emission scenario assumes carbon dioxide concentrations will stabilise at 550 ppm by the mid-21st century. The B1 climate scenario includes rapid economic growth with a population rise to 9 billion in 2050, with slow progression towards a service and information economy with some reductions in material intensity, with likely global temperature increases from 1.1-2.9°C by the end of the 21st C . These scenarios were chosen to reflect a range of possible futures with varying carbon emissions and temperature rises to determine the likely impact of future mitigation strategies on mosquito bionomics.
Daily average meteorological values (standard deviation in parentheses) for Cairns, Queensland for both current climate and future climate (MPI ECHAM 5 model) under two emissions scenarios
Max T (°C)
Ave T (°C)
Min T (°C)
Future B1 (2046-65)
Future A2 (2046-65)
For a focal site of dengue transmission in tropical Australia (Cairns, Queensland), female mosquito abundance, wet weight, and the extrinsic incubation period for dengue virus in Ae. aegypti were estimated using simulations for the recent climate (1990-2011, as an estimate for ‘current climate’) and for both future climate change scenarios using the MPI ECHAM 5 model for the period 2046-2064). The monthly averages for each biological parameter (adult female mosquito abundance, wet weight) were calculated based on CIMSiM outputs. Extrinsic Incubation Period (EIP) of dengue virus in the mosquitoes was obtained for simulation periods using the DENSiM module attached to CIMSiM 3.27. Thirty replicate simulation runs were performed for each scenario.
Monthly data for abundance, wet weight and EIP were compared between the three scenarios (current climate and two future climate scenarios) using Analysis of Variance (STATA Ver 11. Statacorp, College Station TX). Bartlett’s test was used to check for homogeneity of variance  and Scheffe’s test was used a posteriori to determine which scenarios differed . Comparisons were made for the months of January (wet season, during which most DENV transmission occurs in Cairns) and July (dry season, typical low point in transmission).
Average temperature increase of about 0.6°C and approximately 1 mm additional rain per day are projected for Cairns 2046-2064. These projections are not significantly different between the A2 and B1 scenarios, although less rainfall increase was predicted for the A2 scenario (Table 2).
Changes to future climate have been predicted to influence mosquito-borne disease incidence, primarily through alterations to mosquito ecology and pathogen incubation rates and survival. To characterize such changes for tropical Australia, we used a mechanistic modeling approach to explore changes to Ae. aegypti bionomics and dengue virus incubation rate.
A warmer and more humid climate will not necessarily increase mosquito abundance. Predictions for future mosquito abundance are inconsistent, with an increase or decrease under two scenarios (B1) and A2, respectively). Given the uncertainty in magnitude of temperature change under future climate, it makes sense to interpret these as equivocal results, neither predicting an increase or a decrease. A decrease can be predicted with climate warming indicating that simplistic extrapolation based on models incorporating few biological and environmental factors may lead to spurious conclusions regarding vector abundance, and therefore disease incidence.Faster development rates enabling a greater number of generations to be completed per year may result from increased temperature. However, higher temperatures may lead to decreased abundance by increasing the probability of exceeding the optimal temperature range for growth and development, and thereby decreasing survival probabilities. Oviposition rates are predicted to increase under both future climate scenarios (Figure 3), and water levels for rain-filled containers are not predicted to decrease (data not shown). Therefore, under the range of climate encompassed by our modelling, mosquito abundance in a future climate is related to adult survivorship and the probability of successful completion of larval development.
The projected decrease in mosquito abundance under the A2 scenario may seem counter-intuitive, particularly given the projected increase in abundance for the B1 scenario (Figure 1). Higher temperatures are expected to increase the rate of larval development and boost mosquito populations, but if this is not matched by sufficient rainfall (less rainfall increase is predicted for the A2 scenario), the probability of eggs being hatched may decrease. For Ae. aegypti, eggs are laid at the waterline of containers, and high temperatures will cause greater evaporation of this water. Eggs are only hatched once the water level reaches that level again. An increase in temperature not matched by a sufficient rainfall increase may be just one mechanism leading to reduced mosquito abundance under the A2 scenario.
We do not claim here that our modelled representations of mosquito ecology will identically match every possible permutation in the field. Rather, our varying results predicting increases and decreases in mosquito abundance demonstrate that under future climate, both results are possible. This calls into question the overarching assumption that a warmer climate will lead to greater mosquito abundance.
Our modelling predicts Ae. aegypti of smaller body size, regardless of the particular emission scenario. This was not surprising, as mosquitoes will develop more quickly in higher temperatures, and thus have less opportunity to acquire nutrients. The decreased body size in future climate was not quite as evident for the A2 scenario in the late dry season, however (Figure 2), the result of cooler October temperatures predicted for this scenario compared with the B1, permitting longer larval periods. For Ae. aegypti, past research has indicated both decreased  and increased  host feeding with small body size. Whilst fecundity declines with decreasing body size in Ae. aegypti, host contact rates are thought to increase . However, the role of host defensive behavior in these studies was not taken into account. Anopheline mosquitoes with smaller body size take blood more frequently , but the contribution to malaria transmission is thought to decrease with smaller body size . Our current climate simulated body sizes were consistent with those measured in the field in the same region; 3.08 mm simulated compared with 3.01 mm recorded from the field , giving us confidence in our predictions. However, notwithstanding our predictions of smaller Ae. aegypti in a future climate, it is worth noting that mosquito body sizes will vary in the field (standard error of 0.01-0.02 mm ). Thus, whilst we predict overall smaller Ae. aegypti in a future climate, variation will still persist.
Future dengue transmission intensity and risk will be significantly, though certainly not solely, determined by vector ecology. Greater mosquito abundance may increase mosquito-host contact and increase dengue transmission. With decreased extrinsic incubation period, dengue transmission risk would be expected to increase. However, decreased mosquito body size is likely to decrease blood feeding frequency and may reduce dispersal and survival, with consequent lowered dengue incidence. Further, decreased body size could enhance dengue virus infectivity in mosquitoes and dissemination probability . Conversely, whilst decreased mosquito body size in Ae. aegypti does not appear to impact oviposition success  it is likely to decrease fecundity (at least in some months of the year) and have some impact on blood feeding frequency. Smaller body size is also likely to negatively impact dispersal and survival .
The vectorial capacity (VC) equation first proposed by Macdonald  describes the overall ability of a vector species to transmit a pathogen at a particular time and place, with extrinsic virus incubation duration in the mosquito related strongly to increases in VC. Based on these relations a decrease in virus incubation time as predicted here (Figure 4) would increase VC. However, contradictory findings on future mosquito abundance, coupled with smaller body size could both reduce VC.
It is unclear whether dengue incidence will increase in a future climate for our case study location in tropical north Queensland, Australia, if projections are based solely on vector factors. Indeed, projections of future vector numbers and abundance does not lead, inevitably, to projected increase in dengue incidence, and that therefore models require greater sophistication (and perhaps wholly new approaches to determining future risk for vector-borne disease risks).
Future work should focus on modeling future dengue transmission risk using a mechanistic approach which takes into account these predicted changes in mosquito body size and virus incubation rate. The Dengue Simulation model (DENSiM), which utilises outputs from the CIMSiM model used here, could be used for this purpose. DENSiM has been previously used to model inter-annual dengue activity in Australia .
By refining predictions of how dengue transmission may alter in a future climate we are better placed to determine the likely impacts of climate change on human health, so that resources to manage disease can be most effectively allocated. Such information may also be used to garner support for carbon emission reduction strategies. The information presented in this paper permits a more refined examination of how climate change may alter dengue in Australia, and may add rigour to the process of developing adaptive strategies and may inform government policy change.
Tom Scott and Andy Garcia (University of California) provided CIMSiM Version 3.27 and much technical advice and support for its use. Alex Evans, Bertrand Timbal and Yan Wang at the Bureau of Meteorology (Melbourne) were very helpful in obtaining statistically downscaled data for future climate. Melanie Bannister-Tyrrell was pivotal in the continued application of CIMSiM in Australia. This project was funded by the Commonwealth Department for Climate Change, via the NH&MRC (project 1003371).
- Guzman MG, Halstead SB, Artsob H, Buchy P, Farrar J, Gubler DJ, Peeling RW: Dengue: a continuing global threat. Nat Rev Microbiol. 2010, 8: S7-S16. 10.1038/nrmicro2460.PubMed CentralView ArticlePubMedGoogle Scholar
- IPCC: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Edited by: Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM. 2013, Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press, 1535-Google Scholar
- Patz JA, Campbell-Lendrum D, Holloway T, Foley JA: Impact of regional climate change on human health. Nature. 2005, 438 (7066): 310-317. 10.1038/nature04188.View ArticlePubMedGoogle Scholar
- Hales S, De Wet N, Maindonald J, Woodward A: Potential effect of population and climate changes on global distribution of dengue fever: an empirical model. Lancet. 2002, 360: 830-834. 10.1016/S0140-6736(02)09964-6.View ArticlePubMedGoogle Scholar
- Åström C, Rocklöv J, Hales S, Béguin A, Louis V, Sauerborn R: Potential distribution of dengue fever under scenarios of climate change and economic development. Ecohealth. 2012, 9: 448-454. 10.1007/s10393-012-0808-0.View ArticlePubMedGoogle Scholar
- Banu S, Hu W, Hurst C, Tong S: Dengue transmission in the Asia‒Pacific region: impact of climate change and socio‒environmental factors. Trop Med Int Health. 2011, 16: 598-607. 10.1111/j.1365-3156.2011.02734.x.View ArticlePubMedGoogle Scholar
- Van Kleef E, Bambrick H, Hales S: Systematic review of the geographic distribution of dengue fever and the potential influence of global climate change. 2010, TropIKA.net, Genève, 2010. Available from http://journal.tropika.net/scielo.php?script=sci_arttext&pid=S2078-86062010005000001&lng=en&nrm=iso. Accessed 31 Dec 2012Google Scholar
- Russell RC, Currie BJ, Lindsay MD, Mackenzie JS, Ritchie SA, Whelan PI: Dengue and climate change in Australia: predictions for the future should incorporate knowledge from the past. Med J Aust. 2009, 190: 265-268.PubMedGoogle Scholar
- Williams CR, Bader CA, Kearney MR, Ritchie SA, Russell RC: The extinction of dengue through natural vulnerability of its vectors. PLoS Negl Trop Dis. 2010, 4: e922-10.1371/journal.pntd.0000922.PubMed CentralView ArticlePubMedGoogle Scholar
- Kearney M, Porter WP, Williams C, Ritchie S, Hoffmann AA: Integrating biophysical models and evolutionary theory to predict climatic impacts on species’ ranges: the dengue mosquito Aedes aegypti in Australia. Funct Ecol. 2009, 23: 528-538. 10.1111/j.1365-2435.2008.01538.x.View ArticleGoogle Scholar
- Bannister-Tyrrell M, Williams C, Ritchie SA, Rau G, Lindesay J, Mercer G, Harley D: Weather-driven variation in dengue activity in Australia examined using a process-based modeling approach. Am J Trop Med Hyg. 2013, 88: 65-72. 10.4269/ajtmh.2012.11-0451.PubMed CentralView ArticlePubMedGoogle Scholar
- Walton WE, Reisen WK: Influence of climate change on mosquito development and blood-feeding patterns. Viral Infections and Global Change. Edited by: Singh S. 2014, New Jersey: Wiley, 35-56.Google Scholar
- Bader CA, Williams CR: Mating, ovariole number and sperm production of the dengue vector mosquito Aedes aegypti (L.) in Australia: broad thermal optima provide the capacity for survival in a changing climate. Physiol Entomol. 2012, 37: 136-144. 10.1111/j.1365-3032.2011.00818.x.View ArticleGoogle Scholar
- Ritchie SA: Dengue: Australia’s other pandemic. Microbiol Aust. 2009, 30: 114-117.Google Scholar
- Focks DA, Haile DG, Daniels E, Mount GA: Dynamic life table model for Aedes aegypti (Diptera: Culicidae): analysis of the literature and model development. J Med Entomol. 1993, 30: 1003-1017.View ArticlePubMedGoogle Scholar
- Williams CR, Johnson PH, Long SA, Rapley LP, Ritchie SA: Rapid estimation of Aedes aegypti population size using simulation modeling, with a novel approach to calibration and field validation. J Med Entomol. 2008, 45: 1173-1179. 10.1603/0022-2585(2008)45[1173:REOAAP]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Williams CR, Johnson PH, Ball TS, Ritchie SA: Productivity and population density estimates of the dengue vector mosquito Aedes aegypti (Stegomyia aegypti) in Australia. Med Vet Entomol. 2013, 27: 313-322. 10.1111/j.1365-2915.2012.01051.x.View ArticlePubMedGoogle Scholar
- Rašić G, Endersby NM, Williams C, Hoffmann AA: Using Wolbachia-based releases for suppression of Aedes mosquitoes: insights from genetic data and population simulations. Ecol Appl. 2014, in pressGoogle Scholar
- Focks DA, Daniels E, Haile DG, Keesling JE: A simulation model of the epidemiology of urban dengue fever: literature analysis, model development, preliminary validation, and samples of simulation results. Am J Trop Med Hyg. 1995, 53: 489-506.PubMedGoogle Scholar
- Tun‒Lin W, Burkot TR, Kay BH: Effects of temperature and larval diet on development rates and survival of the dengue vector Aedes aegypti in north Queensland, Australia. Med Vet Entomol. 2000, 14: 31-37. 10.1046/j.1365-2915.2000.00207.x.View ArticleGoogle Scholar
- Kay BH, Ryan PA, Lyons SA, Foley PN, Pandeya N, Purdie D: Winter intervention against Aedes aegypti (Diptera: Culicidae) larvae in subterranean habitats slows surface recolonization in summer. J Med Entomol. 2002, 39: 356-361. 10.1603/0022-2585-39.2.356.View ArticlePubMedGoogle Scholar
- Bartlett MS: Properties of sufficiency and statistical tests. Proc R Stat Soc A. 1937, 160: 268-282. 10.1098/rspa.1937.0109.View ArticleGoogle Scholar
- Scheffé H: The Analysis of Variance. 1959, New York: WileyGoogle Scholar
- Siegel JP, Novak RJ, Lampman RL, Steinly BA: Statistical appraisal of the weight–wing length relationship of mosquitoes. J Med Entomol. 1992, 29: 711-714.View ArticlePubMedGoogle Scholar
- Christophers SR: Aedes aegypti (L.) the yellow fever mosquito. 1960, London, UK: Cambridge University PressGoogle Scholar
- Hugo LE, Cook PE, Johnson PH, Rapley LP, Kay BH, Ryan PA, Ritchie SA, O'Neill SL: Field validation of a transcriptional assay for the prediction of age of uncaged Aedes aegypti mosquitoes in northern Australia. PLoS Negl Trop Dis. 2010, 4: e608-10.1371/journal.pntd.0000608.PubMed CentralView ArticlePubMedGoogle Scholar
- Xue RD, Edman JD, Scott TW: Age and body size effects on blood meal size and multiple blood feeding by Aedes aegypti (Diptera: Culicidae). J Med Entomol. 1995, 32: 471-474.View ArticlePubMedGoogle Scholar
- Farjana T, Tuno N: Multiple blood feeding and host-seeking behavior in Aedes aegypti and Aedes albopictus (Diptera: Culicidae). J Med Entomol. 2013, 50: 838-846. 10.1603/ME12146.View ArticlePubMedGoogle Scholar
- Nasci RS: The size of emerging and host-seeking Aedes aegypti and the relation of size to blood-feeding success in the field. J Am Mosq Control Assoc. 1986, 2: 61-62.PubMedGoogle Scholar
- Takken W, Klowden MJ, Chambers GM: Effect of body size on host seeking and blood meal utilization in Anopheles gambiae sensu stricto (Diptera: Culicidae): the disadvantage of being small. J Med Entomol. 1998, 35: 639-645.View ArticlePubMedGoogle Scholar
- Takken W, Smallegange RC, Vigneau AJ, Johnston V, Brown M, Mordue-Luntz AJ, Billingsley PF: Larval nutrition differentially affects adult fitness and Plasmodium development in the malaria vectors Anopheles gambiae and Anopheles stephensi. Parasit Vectors. 2013, 6: 345-10.1186/1756-3305-6-345.PubMed CentralView ArticlePubMedGoogle Scholar
- Alto BW, Reiskind MH, Lounibos LP: Size alters susceptibility of vectors to dengue virus infection and dissemination. Am J Trop Med Hyg. 2008, 79: 688-695.PubMed CentralPubMedGoogle Scholar
- Yeap HL, Axford JK, Popovici J, Endersby NM, Iturbe-Ormaetxe I, Ritchie SA, Hoffmann AA: Assessing quality of life-shortening Wolbachia-infected Aedes aegypti mosquitoes in the field based on capture rates and morphometric assessments. Parasit Vectors. 2014, 7: 58-10.1186/1756-3305-7-58.PubMed CentralView ArticlePubMedGoogle Scholar
- Macdonald G: The Epidemiology and Control of Malaria. 1957, London: Oxford University PressGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.