The Usual Suspects: Comparison of the Relative Roles of Potential Urban Chikungunya Virus Vectors in Australia

Abstract

The global re-emergence of chikungunya virus (CHIKV) over the last decade presents a serious public health risk to Australia. An increasing number of imported cases further underline the potential for local transmission to occur if local mosquitoes bite an infected traveller. Laboratory experiments have identified a number of competent Australian mosquito species, including the primary vectors of CHIKV abroad, Aedes aegypti and Aedes albopictus, and local endemic species Aedes vigilax and Aedes notoscriptus. The implication of these additional endemic species as potential vectors has generated much uncertainty amongst public health professionals regarding their actual role in CHIKV transmission in the field. Using data estimated from or documented in the literature, we parameterise a simple vectorial capacity model to evaluate the relative roles of Australian mosquito species in potential CHIKV transmission. The model takes into account a number of key biological and ecological variables which influence the role of a species in field transmission, including population density, human feeding rates, mosquito survival rates and vector competence. We confirm the relative importance of Ae. aegypti and Ae. albopictus in sustaining potential CHIKV transmission in Australia. Even at maximum estimated densities and human feeding rates, Ae. vigilax and Ae. notoscriptus are likely to play a relatively minor role in CHIKV transmission, when compared with either Ae. aegypti or Ae. albopictus. This relatively straightforward analysis has application for any region where mosquito species have been incriminated in vector competence experiments, but where their actual role in CHIKV transmission has not been established.

About the Authors:

Cassie C. Jansen

* E-mail: [email protected]

Affiliation: Metro North Public Health Unit, Queensland Health, Windsor, Queensland, Australia

Craig R. Williams

Affiliation: Sansom Institute for Health Research, University of South Australia, Adelaide, South Australia, Australia

Andrew F. van den Hurk

Affiliation: Forensic and Scientific Services, Department of Health, Queensland Government, Coopers Plains, Queensland, Australia

Introduction

Chikungunya virus (CHIKV) is a mosquito-borne pathogen that has re-emerged globally in the last decade. First isolated in 1952–53 in Tanzania, the name ‘chikungunya’ describes the painful signs associated with the polyarthritic disease syndrome caused by infection with this virus [1]. The virus is transmitted between humans by mosquitoes, particularly Aedes aegypti and Aedes albopictus in urban settings. Since 2004, CHIKV has spread from East Africa, through the islands of the western Indian Ocean, into the Indian sub-continent and Southeast Asia [2, 3]. In 2007, a traveller from India introduced the virus into north-eastern Italy, causing the first European outbreak of the virus [4]. Emerging for the first time in the Western Hemisphere in late 2013, autochthonous cases have been reported throughout countries in the Caribbean, and in Central and South America [5–8] and, most recently, the first locally transmitted cases in the USA have been recorded in Florida [8, 9]. Outbreaks of CHIKV are often explosive, as evidenced on Reunion Island, where an estimated 244,000 cases represented a prevalence of 35% of the population [10].

Recent outbreaks in nearby locations in the Pacific region including, but not limited to, Papua New Guinea (PNG), New Caledonia, Tonga, Samoa and the Federated states of Micronesia [11–13] have caused considerable concern amongst public health professionals in Australia. Not surprisingly, there have been an increasing number of imported cases into Australia. Between 2002 and 2012, 168 imported cases were reported over the ten-year period [14], but in 2013 alone, 133 cases were notified to the National Notifiable Diseases Surveillance System, followed by 107 in 2014 [15]. If appropriate mosquito vectors are present, local transmission in Australia arising from an imported case is certainly possible, in a scenario similar to what has already occurred in Italy and the Caribbean.

Given this threat, laboratory experiments were undertaken to determine which Australian mosquito species could become infected with and transmit CHIKV [16, 17]. Such physiological competence demonstrates compatibility between a vector species and the virus in question, and is a requisite for a mosquito to be implicated in the transmission of a virus such as CHIKV. These studies demonstrated that a number of endemic species are physiologically competent at becoming infected with, and transmitting CHIKV. Not surprisingly, Australian populations of both Ae. aegypti and Ae. albopictus were highly competent. However, the current distribution of these two species in Australia is relatively limited, with Ae. aegypti occurring predominantly in north Queensland—with isolated populations in central and southern Queensland [18, 19] and Ae. albopictus present only in the Torres Strait Islands to the north of Cape York [20, 21].

In addition to confirming the vector competence of Australian populations of both Ae. aegypti and Ae. albopictus, van den Hurk and others [16] also identified native Australian species, including Aedes vigilax and Aedes notoscriptus as competent vectors. This has important implications, as these species are geographically widespread, and are abundant in urban areas in both tropical and temperate Australia. However, their potential role in field transmission is uncertain, because the virus has never been isolated from them in nature. Further, species-specific biological traits, such as blood feeding behaviour and relative abundance, provide an additional source of uncertainty when predicting their potential roles in field transmission, and have not been considered in detail. Thus, the potential role of various vector species is prone to misinterpretation when assessing risk and during the development of guidelines for CHIKV risk management.

Herein, we apply a simple vectorial capacity model first developed for malaria vectors [22, 23] to evaluate the relative roles of Australian mosquito species in potential CHIKV transmission. The model takes into account a number of key variables such as density in relation to humans, human feeding rates, mosquito survival rates, the incubation period of the virus in the mosquito and the transmission rate.

Materials and Methods

Model

We employed a simple vectorial capacity model presented by Black & Moore [24] and as described by Kramer & Ebel [25]. This equation is based on the “basic reproductive rate” model described and developed by Macdonald [22, 23] and later modified by Garrett-Jones [26]. The model we employed estimates the mean number of infective bites arising per unit time from one infective host and is represented as:

Whereby VC = vectorial capacity; m = vector density in relation to the host; a = probability a vector feeds on a host in one day; p = probability of vectors surviving through one day; n = duration of the extrinsic incubation period (EIP) in days; b = vector competence (the proportion of vectors ingesting an infective meal that successfully become infective to another host), and 1/(-ln(p)) = duration of vector’s life in days after surviving the EIP (where ln denotes the natural logarithm).

Estimation of parameters

Four species were chosen for consideration, based on implication in CHIKV transmission abroad and/or their association with humans and previous demonstration of vector competence in the laboratory [16]. The parameters used to estimate the vectorial capacity (VC) for each species were derived from previously published studies and unpublished data, with particular emphasis placed on observations from Australian populations, when such information was available.

Values for m (density) were identified as biting rates using an extrapolation from human landing collections which measure the number of attempted mosquito bites per human host over a defined period which, in our study, was per hour. The probability that a mosquito feeds on a host in one day (a) was calculated as the product of human host preference, feeding frequency and the number of multiple feeds comprising one replete bloodmeal. Host preference was identified from retrospective host feeding studies of bloodfed field-collected mosquitoes and defined as the proportion of bloodmeals identified from a human host (as opposed to meals from other animals). As no reservoir hosts, other than primates, have been identified for CHIKV, only the feeding rate on humans is applicable to the current study.

As per the Black and Moore definition of vectorial capacity adopted in this study [24], the feeding frequency is approximately equal to one divided by the length of the female reproductive cycle and was estimated from the duration of the gonotrophic cycle (in this case, the period of time between obtaining a bloodmeal and egg laying), after which a female mosquito would seek a subsequent bloodmeal. Vector competence (b) was based on transmission rates identified from laboratory infection studies using a Reunion Island isolate containing the alanine to valine amino acid substitution in the E1 envelope glycoprotein (E1:A226V) [16, 17]. For consistency, we chose day 14 transmission estimates for Ae. albopictus as these were the only values available for the other species. The vectorial capacity was calculated using an EIP of 10 days for all species because a recent study using Ae. albopictus demonstrated that the transmission rate peaked on this day [17]. Notwithstanding the likely variability in EIP due to variation in CHIKV strains and temperature, we selected 10 days as the EIP for the purposes of comparison between vector species as this likely represents the average EIP rather than the shortest observed EIP which may bias comparison between species [27]. Finally, all survival estimates (p) were estimated from field mark-release-recapture studies of wild populations, and represent the probability that a mosquito would survive through one day.

Where considerable differences in estimates of density and/or human feeding rates were found for some species, we chose minimum and maximum values to parameterise the model for comparison within those species.

For the purposes of comparison, it may be assumed that the vectorial capacity values reported here are valid for conditions of 28°C. Vector competence experiments for CHIKV [16, 17] were conducted at this temperature, and many field observations of host landing rates and daily survival are made during the summer months when temperatures are typically in the high 20s°C in Australian cities.

Results

Aedes aegypti, Ae. albopictus, Ae. notoscriptus and Ae. vigilax were chosen for calculation of VC based on their association with humans, and demonstrated vector competence in the laboratory. The values identified from the literature and used to parameterise the model for each species are listed in Table 1, along with their sources. All parameter estimates pertaining to endemic species Ae. notoscriptus and Ae. vigilax were sourced from Australian studies, whereas estimates of parameters comprising a and p for Ae. aegypti and Ae. albopictus were, where necessary, sourced from studies conducted on mosquito populations from locations other than Australia (Table 1). The length of the gonotrophic cycle was designated as three days for all species, as these data were not available for all species, but was consistently identified as three days in most cases. Importantly, we acknowledge that these parameters are, to some extent, temperature-dependent, and the values used here will be subject to the laboratory and field conditions at which data were collected in source studies.

[Figure omitted. See PDF.]

Table 1. Parameters identified from the literature and used to estimate vectorial capacity of Australian mosquito species for chikungunya virus.

https://doi.org/10.1371/journal.pone.0134975.t001

The influence of multiple feeds was only considered for Ae. aegypti as this behaviour has not been described for either Ae. vigilax or Ae. notoscriptus. While Ae. albopictus may occasionally feed on more than one host to complete a blood meal [33], it is anticipated that this behaviour occurs in less than 20% of bloodmeals [40] and, even when it does occur, we consider the likelihood of resuming a blood meal on a human to be low due to more eclectic feeding behaviour (a propensity to feed on a range of host types other than humans) when compared with Ae. aegypti.

Minimum and maximum estimates of host feeding rates (host preference) are based on the lowest and highest records identified in the literature. For Ae. albopictus, the maximum human host feeding rate identified was 0.957 [34] and the minimum was 0.2 [33], while the maximum identified for Ae. aegypti was 0.95 [32] and minimum was 0.75 [31]. Importantly, no Australian records of blood feeding habits or host preference were available for Ae. albopictus. The maximum and minimum densities identified from the literature for Ae. vigilax were 550.07 [30] and 8.2 [29], respectively, and for Ae. albopictus 48 and 3.6 (Queensland Health, state government data), respectively.

Overall Results

When maximum and minimum human host feeding rates were employed in the model, Ae. aegypti demonstrated a VC of 3.00 and 4.85, respectively, which was considerably greater than that of Ae. notoscriptus and Ae. vigilax under all conditions, and that of Ae. albopictus under three of the four tested conditions (Fig 1). Even with minimum host feeding rates employed in the model, Ae. aegypti consistently demonstrated a much higher VC than all of the other species, under all conditions. When entered into the model with maximum density and human host feeding rate estimates, Ae. albopictus demonstrated a VC of 1.44, which was the highest estimate of the remaining species. However, when minimum density and/or host feeding rates were used, VC values for Ae. albopictus were consistently less than 0.15. Values of VC calculated for both Ae. vigilax and Ae. notoscriptus were less than 0.22 under all conditions. Even with a maximum density estimate of over 500 bites per person per day (the greatest estimate for any of the species considered), the VC of Ae. vigilax was 75% less than that of Ae. aegypti. Similarly, even at the maximum human host feeding rate, Ae. notoscriptus demonstrated a VC of < 0.02.

[Figure omitted. See PDF.]

Fig 1. Comparison of vectorial capacity for four species, given maximum and minimum estimates of density and/or host feeding preferences where data available from the literature (refer Table 1).

https://doi.org/10.1371/journal.pone.0134975.g001

Discussion

Preparedness for potential urban transmission of CHIKV in Australia is critical. A key element of a pre-emptive management plan is the incrimination of potential local vectors, with consideration of their vector competence together with biological and ecological attributes, which determine their relative roles in field transmission. This activity can enhance preparedness measures, in addition to reducing confusion when interpreting results from laboratory studies in a public health setting. Using a simple vectorial capacity model, we compared the relative potential roles of urban mosquito species in field transmission of CHIKV. We chose four species for consideration, based on laboratory vector competence experiments [16], widespread geographical distributions, association with urban habitats, and the relative abundance of these species in such locations.

Unsurprisingly, the model shows that Ae. aegypti is a key candidate vector of CHIKV in Australia regardless of whether minimum or maximum estimates of human feeding rates were employed, and suggests that Ae. aegypti would most likely be the primary potential Australian vector across most conditions. This species was the first species to be implicated in urban CHIKV transmission and, despite the incrimination of Ae. albopictus in a number of recent outbreaks (particularly since the virus mutation described by Schuffenecker et al. [49] was associated with enhanced infectivity in Ae. albopictus;[50]), Ae. aegypti remains an important vector in a number of locations [51, 52], and has been confirmed to be equally competent for the recently emerged strain of CHIKV [53]. The important role of Ae. aegypti in the transmission of pathogens that are maintained in a human-mosquito-human transmission cycle has been repeatedly demonstrated using similar models to that described herein [43, 54, 55]. Notably, as the probability of the vector feeding on a host in a day (a) appears as a squared element in the model, the VC is very sensitive to changes in the values for a. Thus, consistent anthropophilly throughout its urban distribution, combined with the propensity for Ae. aegypti to obtain multiple bloodmeals in a replete meal (two unique features amongst the species considered herein), can explain why Ae. aegypti is consistently implicated as the most important candidate vector in the current study under all conditions.

Currently, the geographical distribution of Ae. aegypti in Australia is restricted to northern coastal areas and some townships of central and southern Queensland. In Northern Queensland Ae. aegypti periodically transmits dengue viruses, following the importation of virus by travellers that have acquired the infection overseas. In a similar way, it is certainly feasible that Ae. aegypti could be a likely candidate vector of imported CHIKV throughout its geographic distribution, and these regions should consider appropriate control strategies.

Vectorial capacity changes are most sensitive to the daily survival probability of the vector (p), and the extrinsic incubation rate of the virus (n) [24]. Furthermore, changes to CHIKV characteristics that impact EIP will influence VC of Ae. aegypti. Christofferson and others [27] showed the effect of EIP on transmission potential was much greater than the effect of vector competence and that divergent strains of CHIKV vary in terms of the observed EIP between and within vector species. Similarly, even slight increases in daily mosquito survival can lead to disproportionately large increases in VC. Thus, if CHIKV were to become established in Australia, it would immediately be critical to examine both the EIP in local populations and confirm the typical survival rates of candidate species under local conditions as a matter of priority.

Using estimates of maximum density and maximum human host feeding rate, Ae. albopictus was also shown to be a potential Australian candidate vector for CHIKV, with VC estimates around half that calculated for Ae. aegypti under its conservative (minimum estimate) conditions. This suggests that Ae. albopictus could present a significant risk in urban areas where alternative hosts may be scarce, and if domestic habitats are abundant and productive. Currently, the geographical distribution of Ae. albopictus is restricted to the Torres Strait region, but there is considerable concern that this species may establish on mainland Australia [21, 56, 57]. The Torres Strait should be considered a high risk location for potential CHIKV transmission due to the presence of both Ae. albopictus and Ae. aegypti, in addition to its close proximity to PNG which has recently experienced explosive CHIKV activity [12]. Indeed, PNG was the source of the majority of imported chikungunya cases in Queensland during the outbreak (11 of 14 cases reported in Queensland in 2013; D. Francis, Queensland Health, pers. comm.). Importantly, it must also be recognised that management of vector borne diseases in Torres Strait region is particularly challenging due to the relative remoteness of the region and a lack of data describing behavioural characteristics and population dynamics of local Ae. albopictus populations, as this species was only detected in the Torres Strait relatively recently [21]. Our analysis suggests that when at high densities, Ae. albopictus may be an important vector of CHIKV, provided that humans are the predominant source of bloodmeals. As observed for Ae. aegypti, the human host feeding rate can dramatically influence their projected VC, and it is noteworthy that in the current study estimates of human feeding rates of Ae. albopictus were sourced from studies conducted abroad, as host feeding patterns of Australian populations have not yet been examined. Furthermore, should multiple feeding on different human hosts by Ae. albopictus occur, this would increase VC strongly. If field data demonstrating this were available, recalculated VC for this species would likely approach that calculated for Ae. aegypti. Overall, the implication of Ae. albopictus as a potential vector of CHIKV in Australia underscores the urgent need to obtain host feeding data describing the behaviour of Australian populations of Ae. albopictus to further resolve the potential role of this species in CHIKV transmission.

In addition to human feeding rates, density was also shown to be a key determinant of estimated VC for Ae. albopictus by the current model. Importantly, the density of Ae. albopictus varies considerably between islands of the Torres Strait and the current study underscores the importance of understanding such variation across islands to support decision making on vector control with respect to CHIKV transmission risk. Further, should Ae. albopictus manage to establish in southern regions on the Australian mainland [58, 59], the VC of Ae. albopictus for CHIKV transmission in temperate Australia would need to be re-evaluated cognisant of the temperature-dependent nature of CHIKV EIP and mosquito biology and ecology.

Confirmation of the relative importance of Ae. aegypti and Ae. albopictus in the potential transmission of CHIKV in Australia demonstrated here highlights the value of efforts to prevent these species from expanding their geographical ranges and underscores the importance of targeted surveillance activities to provide early detection of the incursion of these vectors in novel locations. This is of particular importance to urban population centres, particularly those which experience high numbers of imported case notifications. Should populations of Ae. aegypti or Ae. albopictus become established in such urban centres on the Australian mainland, it would represent a significant risk of explosive CHIKV transmission, potentially placing a considerable proportion of the human population at risk of infection.

The remaining species considered herein are likely to play a relatively minor role in potential CHIKV transmission, when compared with either Ae. aegypti or Ae. albopictus. For example, despite the ability of Ae. vigilax to reach very high population densities particularly after high tides during humid conditions, this species would likely play a limited role in transmission, according to the current model. Nevertheless, at the highest population density considered, this species would likely play a larger role in urban transmission of CHIKV than Ae. notoscriptus despite the tendency for Ae. notoscriptus to exploit backyard habitats. Indeed, the model suggests that Ae. notoscriptus would have a minimal role in the potential transmission of CHIKV in Australia. However, recent genetic data suggests that Ae. notoscriptus may comprise at least three species in Australia [60], so geographical differences in behaviour, ecology and biology should be further considered. However, the current model suggests that neither Ae. vigilax nor Ae. notoscriptus alone would likely maintain large and explosive outbreaks similar to those observed elsewhere in recent years. This may be viewed as propitious, given the widespread geographical distribution of these two species across all states of Australia, including prevalence in capital cities, where the majority of imported cases of CHIKV are reported.

Notwithstanding, a number of other factors would likely contribute to the likelihood and magnitude of a CHIKV outbreak in Australia. Of these, the number of imported cases in a given region is arguably among the most important. Imported cases represent an opportunity for local transmission if the case is viraemic and suitable vectors are present. Imported cases of CHIKV have been recorded in all Australian states and mainland territories (with the exception of the Australian Capital Territory where CHIKV is not considered a notifiable condition), with the majority of cases reported in the most populous states of New South Wales, Victoria, Queensland and Western Australia [15]. Of these, Queensland is the only state with established populations of Ae. aegypti and/or Ae. albopictus. Importantly, the Cairns local government area (the location of the majority of local dengue transmission) has recorded the second highest number of total CHIKV importations between 2010–2014 in Queensland, outnumbered only by the Brisbane metropolitan region, highlighting the potential vulnerability of northern Queensland where these candidate vectors are present [61].

Importantly, the strain of virus can determine the likelihood of local transmission, and the potential magnitude of ensuing local transmission. Molecular characterization of an Indian Ocean CHIKV strain revealed that the E1:A226V mutation led to increased replication and transmission efficiency in Ae. albopictus which likely contributed to the scale of the Indian Ocean epidemic [50, 53, 62]. Subsequent substitutions in the E2 envelope glycoprotein in some virus populations have further enhanced CHIKV fitness in Ae. albopictus [63, 64]. Interestingly, these substitutions did not confer a fitness advantage in Ae. aegypti. In another study, Vega-Rua et al. [65] demonstrated considerable variation within and between North and South American populations of Ae. aegypti and Ae. albopictus in their ability to become infected with and transmit strains of CHIKV with and without the E1:A226V mutation. The data we used to estimate transmission rates was sourced from vector competence experiments conducted using a Reunion Island isolate containing the E1:A226V mutation [16, 17], and is the only strain that has been characterized in Australian mosquitoes. Given the emergence of the Asian lineage of CHIKV in the Americas and some locations in the Pacific, vector competence experiments are urgently required to assess the ability of Australian mosquito populations to transmit strains of CHIKV which do not contain the E1:A226V mutation. The recalculated VC of incriminated species for other strains of CHIKV should be used to assess the potential for outbreaks caused by different CHIKV lineages.

Australia possesses a number of unique attributes which may impact the likelihood of CHIKV causing explosive transmission similar to that which has occurred in the Indian Ocean region and the Americas. Importantly, CHIKV is antigenically related to endemic Ross River virus (RRV), which circulates in an enzootic cycle in all mainland states and is responsible for the majority of arbovirus disease notifications in Australia [15, 66]. Given that marsupials are the primary hosts of RRV [66], it is possible that marsupials could develop sufficient CHIKV viraemia to infect mosquitoes. Anti-RRV antibodies can reduce the viraemia and severity of disease associated with CHIKV infection in mice [67], so it is plausible that humans with prior infection with RRV may be protected against CHIKV infection and associated disease. Concomitantly, a reduction in viraemia levels could reduce the pool of available hosts capable of infecting recipient mosquitoes and perpetuating transmission. However, the 2012–2013 epidemic in PNG (where RRV is likewise endemic) suggests that whilst prior exposure to RRV may modulate CHIKV infection, it will not prevent transmission.

Public health professionals need to exercise caution when interpreting the results from laboratory vector competence studies, particularly when inferring the potential role of endemic mosquito species in potential field transmission of exotic viruses and developing pre-emptive management strategies. The vectorial capacity model is a tool to aid the understanding of the relative potential importance of vector species in the transmission of a given pathogen, and provides a relative quantitative index of a mosquito population’s capacity to transmit the pathogen to a susceptible host population. A vectorial capacity model can also highlight key biological and ecological factors which require further characterisation to underpin informed risk assessment. Importantly, this model can be applied to assess the role of endemic species in the transmission of other viruses of public health concern to Australia, including the recently emergent Zika virus [68, 69].

When applied to candidate Australian vectors of CHIKV, the VC model suggests that Ae. aegypti and Ae. albopictus are, overwhelmingly, the most likely species to sustain local transmission of CHIKV. Secondary vector species are likely to be able to maintain local transmission, but unlikely to sustain outbreaks of an explosive nature as observed in recent outbreaks abroad. Nevertheless, their potential role in urban transmission of CHIKV cannot be discounted and the biological and behavioural characteristics of local populations must be considered when developing control strategies or management plans. Although our study focussed on the Australian context, the relatively simple analysis presented herein has application to other regions of the world where vector competence experiments have implicated potential vectors, but for which there is limited information on their actual role in CHIKV transmission cycles.

Acknowledgments

Authors thank Stephen Frances, Odwell Muzari and Joe Davis for providing biting rate data. Authors also acknowledge Leon Hugo for advice on calculating survival rates, and Cameron Webb, Richard Russell, John Mackenzie and Phil Lounibos for general advice. We also thank Scott Ritchie for commenting on an early draft of the manuscript.

Author Contributions

Conceived and designed the experiments: CCJ AFvdH. Performed the experiments: CCJ AFvdH CRW. Analyzed the data: CCJ AFvdH CRW. Contributed reagents/materials/analysis tools: CCJ AFvdH CRW. Wrote the paper: CCJ AFvdH CRW.

Citation: Jansen CC, Williams CR, van den Hurk AF (2015) The Usual Suspects: Comparison of the Relative Roles of Potential Urban Chikungunya Virus Vectors in Australia. PLoS ONE 10(8): e0134975. https://doi.org/10.1371/journal.pone.0134975

References

1. Robinson MC. An epidemic of virus disease in Southern Province, Tanganyika Territory, in 1952–53. I. Clinical features. Trans R Soc Trop Med Hyg. 1955;49(1): 28–32. pmid:14373834.

2. Burt FJ, Rolph MS, Rulli NE, Mahalingam S, Heise MT. Chikungunya: a re-emerging virus. Lancet. 2012;379(9816): 662–71. pmid:22100854.

3. Schwartz O, Albert ML. Biology and pathogenesis of chikungunya virus. Nat Rev Microbiol. 2010;8(7): 491–500. pmid:20551973.

4. Rezza G, Nicoletti L, Angelini R, Romi R, Finarelli AC, Panning M, et al. Infection with chikungunya virus in Italy: an outbreak in a temperate region. Lancet. 2007;370(9602): 1840–6. pmid:18061059.

5. Leparc-Goffart I, Nougairede A, Cassadou S, Prat C, de Lamballerie X. Chikungunya in the Americas. Lancet. 2014;383(9916): 514. pmid:24506907.

6. Enserink M. Infectious diseases. Crippling virus set to conquer Western Hemisphere. Science. 2014;344(6185): 678–9. pmid:24833366.

7. Lindsey NP, Prince HE, Kosoy O, Laven J, Messenger S, Staples JE, et al. Chikungunya virus infections among travelers- United States, 2010–2013. Am J Trop Med Hyg. 2015;92(1): 82–7. pmid:25349374

8. Higgs S, Vanlandingham D. Chikungunya Virus and Its Mosquito Vectors. Vector Borne Zoonotic Dis. 2015. pmid:25674945.

9. Staples JE, Fischer M. Chikungunya virus in the Americas—what a vectorborne pathogen can do. N Engl J Med. 2014;371(10): 887–9. pmid:25184860.

10. Renault P, Solet JL, Sissoko D, Balleydier E, Larrieu S, Filleul L, et al. A major epidemic of chikungunya virus infection on Reunion Island, France, 2005–2006. Am J Trop Med Hyg. 2007;77(4): 727–31. pmid:17978079.

11. Horwood P, Bande G, Dagina R, Guillaumot L, Aaskov J, Paviln B. The threat of chikungunya in Oceania. WPSAR. 2013;4(2).

12. Horwood PF, Reimer LJ, Dagina R, Susapu M, Bande G, Katusele M, et al. Outbreak of chikungunya virus infection, Vanimo, Papua New Guinea. Emerg Infect Dis. 2013;19(9): 1535–8. pmid:23965757; PubMed Central PMCID: PMC3810919.

13. Roth A, Mercier A, Lepers C, Hoy D, Duituturaga S, Benyon E, et al. Concurrent outbreaks of dengue, chikungunya, and Zika virus infections—an unprecedented epidemic wave of mosquito-borne viruses in the Pacific 2012–2014. Euro Surveill. 2014;19(41): 20929. pmid:25345518

14. Viennet E, Knope K, Faddy HM, Williams C, Harley D. Assessing the threat of chikungunya virus emergence in Australia. Commun Dis Intell. 2013;37(2): E136–E43.

15. Australian Government Department of Health. National Notifiable Diseases Surveillance System 2014 [cited 2015 7 March 2015]. Available from: http://www9.health.gov.au/cda/source/rpt_2_sel.cfm.

16. van den Hurk AF, Hall-Mendelin S, Pyke AT, Smith GA, Mackenzie JS. Vector competence of Australian mosquitoes for chikungunya virus. Vector Borne Zoonotic Dis. 2010;10(5): 489–95. pmid:19877822.

17. Nicholson J, Ritchie SA, van den Hurk AF. Aedes albopictus (Diptera: Culicidae) as a potential vector of endemic and exotic arboviruses in Australia. J Med Entomol. 2014;51(3): 661–9. pmid:24897860

18. Mottram P, Gould D, Florian P. Aedes aegypti distribution in southern Queensland 1990–2009, are they really marching toward Brisbane? Mosquito Bites in the Asia Pacific Region. 2009;4(1): 45–50.

19. Beebe NW, Cooper RD, Mottram P, Sweeney AW. Australia's dengue risk driven by human adaptation to climate change. PLoS Negl Trop Dis. 2009;3(5): e429. pmid:19415109; PubMed Central PMCID: PMC2671609.

20. Ritchie SA, Moore P, Carruthers M, Williams C, Montgomery B, Foley P, et al. Discovery of a widespread infestation of Aedes albopictus in the Torres Strait, Australia. J Am Mosq Control Assoc. 2006;22(3): 358–65. Epub 2006/10/28. pmid:17067032.

21. Beebe NW, Ambrose L, Hill LA, Davis JB, Hapgood G, Cooper RD, et al. Tracing the tiger: population genetics provides valuable insights into the Aedes (Stegomyia) albopictus invasion of the Australasian Region. PLoS Negl Trop Dis. 2013;7(8): e2361. pmid:23951380; PubMed Central PMCID: PMC3738475.

22. Macdonald G. The Epidemiology and Control of Malaria. London: Oxford University Press; 1957.

23. Macdonald G. Epidemiologic models in studies of vectorborne diseases. Public health reports. 1961;76: 753–64. pmid:13764730; PubMed Central PMCID: PMC1929763.

24. Black WI, Moore CG. Population biology as a tool for studying vector-borne diseases. In: Beaty BJ, Marquards WC, editors. The Biology of Disease Vectors. Niwot: University Press of Colorado; 1996.

25. Kramer LD, Ebel GD. Dynamics of flavivirus infection in mosquitoes. Adv Virus Res. 2003;60: 187–232. pmid:14689695.

26. Garrett-Jones C. Prognosis for interruption of malaria transmission through assessment of the mosquito's vectorial capacity. Nature. 1964;204: 1173–5. pmid:14268587.

27. Christofferson RC, Chisenhall DM, Wearing HJ, Mores CN. Chikungunya viral fitness measures within the vector and subsequent transmission potential. PLoS One. 2014;9(10): e110538. pmid:25310016; PubMed Central PMCID: PMC4195746.

28. Kay BH, Watson TM, Ryan PA. Definition of productive Aedes notoscriptus (Diptera: Culicidae) habitats in western Brisbane, and a strategy for their control. Aust J Entomol. 2008;47: 142–8. pmid:ISI:000255418900010.

29. Jones RE, Barker-Hudson P, Kay BH. Comparison of dry ice baited light traps with human bait collections for surveillance of mosquitoes in northern Queensland, Australia. J Am Mosq Control Assoc. 1991;7(3): 387–94. pmid:1686444.

30. Frances SP, MacKenzie DO, Rowcliffe KL, Corcoran SK. Comparative field evaluation of repellent formulations containing deet and IR3535 against mosquitoes in Queensland, Australia. J Am Mosq Control Assoc. 2009;25(4): 511–3. pmid:20099600.

31. Jansen CC, Webb CE, Graham GC, Craig SB, Zborowski P, Ritchie SA, et al. Blood sources of mosquitoes collected from urban and peri-urban environments in eastern Australia with species-specific molecular analysis of avian blood meals. Am J Trop Med Hyg. 2009;81(5): 849–57. Epub 2009/10/29. pmid:19861621.

32. Scott TW, Chow E, Strickman D, Kittayapong P, Wirtz RA, Lorenz LH, et al. Blood-Feeding Patterns of Aedes aegypti (Diptera, Culicidae) Collected in a Rural Thai Village. J Med Entomol. 1993;30(5): 922–7. pmid:ISI:A1993LV97700014.

33. Richards SL, Ponnusamy L, Unnasch TR, Hassan HK, Apperson CS. Host-feeding patterns of Aedes albopictus (Diptera: Culicidae) in relation to availability of human and domestic animals in suburban landscapes of central North Carolina. J Med Entomol. 2006;43(3): 543–51. pmid:16739414; PubMed Central PMCID: PMC2577020.

34. Valerio L, Marini F, Bongiorno G, Facchinelli L, Pombi M, Caputo B, et al. Host-feeding patterns of Aedes albopictus (Diptera: Culicidae) in urban and rural contexts within Rome province, Italy. Vector Borne Zoonotic Dis. 2010;10(3): 291–4. Epub 2009/06/03. pmid:19485771.

35. Kay BH, Boyd AM, Ryan PA, Hall RA. Mosquito feeding patterns and natural infection of vertebrates with Ross River and Barmah Forest viruses in Brisbane, Australia. Am J Trop Med Hyg. 2007;76(3): 417–23. pmid:17360861.

36. Watson TM. Ecology and behaviour of Aedes notoscriptus (Skuse): implications for arbovirus transmission in southeast Queensland [Ph.D. Thesis]. St Lucia, Queensland: University of Queensland; 1998.

37. Pant CP, Yasuno M. Field studies on the gonotrophic cycle of Aedes aegypti in Bangkok, Thailand. J Med Entomol. 1973;10(2): 219–23. pmid:4707760.

38. Hien DS. Biology of Aedes aegypti (L.,1762) and Aedes albopictus (Skuse, 1895) (Diptera, Culicidae). V. The gonotrophic cycle and oviposition. Acta Parasitologica Polonica. 1976;24(6): 37–55.

39. Christophers SR. Aëdes aegypti (L.), the yellow fever mosquito; its life history, bionomics, and structure. Cambridge Eng.: University Press; 1960. xii, 738 p.

40. Hawley WA. The biology of Aedes albopictus. J Am Mosq Control Assoc Suppl. 1988;1: 1–39. pmid:3068349.

41. Watson TM, Marshall KL, Kay BH. Colonization and laboratory biology of Aedes notoscriptus from Brisbane, Australia. J Am Mosq Control Assoc. 2000;16(2): 138–42. Epub 2000/07/20. pmid:10901637.

42. Lee DJ, Hicks MM, Griffiths M, Russell RC. The Culicidae of the Australasian region. Volume 3. Lee DJ, Commonwealth Institute of Health, editors. Canberra: Australian Government Publishing Service; 1984.

43. Jetten TH, Focks DA. Potential changes in the distribution of dengue transmission under climate warming. Am J Trop Med Hyg. 1997;57(3): 285–97. pmid:9311638.

44. Muir LE, Kay BH. Aedes aegypti survival and dispersal estimated by mark-release-recapture in northern Australia. Am J Trop Med Hyg. 1998;58(3): 277–82. Epub 1998/04/18. pmid:9546403.

45. Maciel-de-Freitas R, Codeco CT, Lourenco-de-Oliveira R. Daily survival rates and dispersal of Aedes aegypti females in Rio de Janeiro, Brazil. Am J Trop Med Hyg. 2007;76(4): 659–65. Epub 2007/04/12. doi: 76/4/659 [pii]. pmid:17426166.

46. Watson TM, Saul A, Kay BH. Aedes notoscriptus (Diptera: Culicidae) survival and dispersal estimated by mark-release-recapture in Brisbane, Queensland, Australia. J Med Entomol. 2000;37(3): 380–4. Epub 2004/11/13. pmid:15535581.

47. Russell RC. The mosquito fauna of Conjola state forest on the South coast of New South Wales. Part 3. Female survival and population age composition. Gen Appl Ent. 1987;19: 25–30.

48. Cousineau MM. Gene flow in the salt-marsh mosquito Aedes vigilax (Skuse) [Ph.D. Thesis]: University of Queensland; 1994.

49. Schuffenecker I, Iteman I, Michault A, Murri S, Frangeul L, Vaney MC, et al. Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak. PLoS Med. 2006;3(7): e263. pmid:16700631; PubMed Central PMCID: PMC1463904.

50. Tsetsarkin KA, Chen R, Sherman MB, Weaver SC. Chikungunya virus: evolution and genetic determinants of emergence. Curr Opin Virol. 2011;1(4): 310–7. pmid:21966353; PubMed Central PMCID: PMC3182774.

51. Jupp PG, McIntosh BM. Chikungunya virus disease. In: Monath TP, editor. The arboviruses: epidemiology and ecology. II. Boca Raton, FL: CRC Press; 1988. p. 137–57.

52. Weaver SC. Arrival of chikungunya virus in the new world: prospects for spread and impact on public health. PLoS Negl Trop Dis. 2014;8(6): e2921. pmid:24967777; PubMed Central PMCID: PMC4072586.

53. Dubrulle M, Mousson L, Moutailler S, Vazeille M, Failloux AB. Chikungunya virus and Aedes mosquitoes: saliva is infectious as soon as two days after oral infection. PLoS One. 2009;4(6): e5895. pmid:19521520; PubMed Central PMCID: PMC2690823.

54. Anderson JR, Rico-Hesse R. Aedes aegypti vectorial capacity is determined by the infecting genotype of dengue virus. Am J Trop Med Hyg. 2006;75(5): 886–92. pmid:17123982; PubMed Central PMCID: PMC1993907.

55. 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; PubMed Central PMCID: PMC3088608.

56. Russell RC, Williams CR, Sutherst RW, Ritchie SA. Aedes (Stegomyia) albopictus—a dengue threat for southern Australia? Commun Dis Intell. 2005;29(3): 296–8. Epub 2005/10/14. pmid:16220868.

57. van den Hurk AF, Craig SB, Tulsiani SM, Jansen CC. Emerging tropical diseases in Australia. Part 4. Mosquitoborne diseases. Ann Trop Med Parasitol. 2010;104(8): 623–40. pmid:ISI:000285292100002.

58. Hill MP, Axford J, Hoffmann AA. Predicting the spread of Aedes albopictus in Australia under current and future climates: Multiple approaches and datasets to incorporate potential evolutionary divergence. Austral Ecol. 2014(39): 469–78.

59. Nicholson J, Ritchie SA, Russell RC, Zalucki MP, Hurk AFVD. Ability for Aedes albopictus (Diptera: Culicidae) to survive at the climatic limits of its potential range in eastern Australia. J Med Entomol. 2014;51(5): 948–57. pmid:25276922

60. Endersby NM, White VL, Chan J, Hurst T, Rasic G, Miller A, et al. Evidence of cryptic genetic lineages within Aedes notoscriptus (Skuse). Infect Genet Evol. 2013;18: 191–201. pmid:23681021.

61. Queensland Government. Queensland Notifiable Conditions System (NOCS) 2015 [cited 2015 31 March].

62. Vazeille M, Moutailler S, Coudrier D, Rousseaux C, Khun H, Huerre M, et al. Two Chikungunya isolates from the outbreak of La Reunion (Indian Ocean) exhibit different patterns of infection in the mosquito, Aedes albopictus. PLoS One. 2007;2(11): e1168. pmid:18000540; PubMed Central PMCID: PMC2064959.

63. Tsetsarkin KA, Weaver SC. Sequential adaptive mutations enhance efficient vector switching by Chikungunya virus and its epidemic emergence. PLoS Pathog. 2011;7(12): e1002412. pmid:22174678; PubMed Central PMCID: PMC3234230.

64. Tsetsarkin KA, Chen R, Yun R, Rossi SL, Plante KS, Guerbois M, et al. Multi-peaked adaptive landscape for chikungunya virus evolution predicts continued fitness optimization in Aedes albopictus mosquitoes. Nature communications. 2014;5: 4084. pmid:24933611.

65. Vega-Rua A, Zouache K, Girod R, Failloux AB, Lourenco-de-Oliveira R. High level of vector competence of Aedes aegypti and Aedes albopictus from ten American countries as a crucial factor in the spread of Chikungunya virus. Journal of virology. 2014;88(11): 6294–306. pmid:24672026.

66. Russell RC. Ross River virus: ecology and distribution. Annu Rev Entomol. 2002;47: 1–31. pmid:11729067.

67. Gardner J, Anraku I, Le TT, Larcher T, Major L, Roques P, et al. Chikungunya virus arthritis in adult wild-type mice. Journal of virology. 2010;84(16): 8021–32. pmid:20519386; PubMed Central PMCID: PMC2916516.

68. Pyke AT, Daly MT, Cameron JN, Moore PR, Taylor CT, Hewitson GR, et al. Imported zika virus infection from the Cook Islands into Australia, 2014. PLoS currents. 2014;6. pmid:24944843; PubMed Central PMCID: PMC4055592.

69. Musso D, Nilles EJ, Cao-Lormeau VM. Rapid spreading of emerging Zika virus in the Pacific area. Clin Microbiol Infect. 2014;20: 595–6. pmid:24909208.