Skip to main content

Probable contribution of Culex quinquefasciatus mosquitoes to the circulation of chikungunya virus during an outbreak in Mombasa County, Kenya, 2017–2018

Abstract

Background

Chikungunya virus is an alphavirus, primarily transmitted by Aedes aegypti and Ae. albopictus. In late 2017–2018, an outbreak of chikungunya occurred in Mombasa county, Kenya, and investigations were conducted to establish associated entomological risk factors.

Methods

Homes were stratified and water-filled containers inspected for immature Ae. aegypti, and larval indices were calculated. Adult mosquitoes were collected in the same homesteads using BG-Sentinel and CDC light traps and screened for chikungunya virus. Experiments were also conducted to determine the ability of Culex quinquefasciatus to transmit chikungunya virus.

Results

One hundred thirty-one houses and 1637 containers were inspected; 48 and 128 of them, respectively, were positive for immature Ae. aegypti, with the house index (36.60), container index (7.82) and Breteau index (97.71) recorded. Jerry cans (n = 1232; 72.26%) and clay pots (n = 2; 0.12%) were the most and least inspected containers, respectively, while drums, the second most commonly sampled (n = 249; 15.21%), were highly positive (65.63%) and productive (60%). Tires and jerry cans demonstrated the highest and lowest breeding preference ratios, 11.36 and 0.2, respectively. Over 6900 adult mosquitoes were collected and identified into 15 species comprising Cx. quinquefasciatus (= 4492; 65.04%), Aedes vittatus (n = 1137; 16.46%) and Ae. aegypti (n = 911; 13.19%) and 2 species groups. Simpson’s dominance and Shannon-Wiener diversity indices of 0.4388 and 1.1942 were recorded, respectively. Chikungunya virus was isolated from pools of Ae. aegypti (1) and Cx. quinquefasciatus (4), two of which were males. Minimum infection rates of 3.0 and 0.8 were observed for female Ae. aegypti and Cx. quinquefasciatus, respectively. Between 25 and 31.3% of exposed mosquitoes became infected with CHIKV 7, 14 and 21 days post-exposure. For the experimentally infected Cx. quinquefasciatus mosquitoes, between 13 and 40% had the virus disseminated, with 100% transmission being observed among those with disseminated infection.

Conclusions

These results demonstrated high risk of chikungunya transmission for residents in the sampled areas of Mombasa. Transmission data confirmed the probable role played by Cx. quinquefasciatus in the outbreak while the role of Ae. vittatus in the transmission of chikungunya virus remains unknown.

Graphical Abstract

Introduction

Chikungunya virus (CHIKV) is a global re-emerging mosquito-borne alphavirus, which was first detected in 1952 along the Tanzania-Mozambique border [1]. In Kenya, a chikungunya re-mergence was reported during a major outbreak in 2004–2005 in the coastal region starting in Lamu and spreading to the Indian Ocean islands [2]. This was followed by a huge outbreak in Mandera City, after a lull of 12 years, which affected 10.31% of the population [3], and a year later, another outbreak was reported in Mombasa in 2017/2018.

Historically, CHIKV has been endemic in tropical and subtropical regions of sub-Saharan Africa and Southeast Asia, where two distinct transmission cycles, sylvatic and urban, exist. In Africa, the mosquito species that are traditionally involved in the urban transmission cycle include Aedes aegypti and, more recently, Ae. albopictus [4]. Ae. aegypti is widely distributed in Kenya [5] while Ae. albopictus has not been reported. In the sylvatic transmission cycle, a wider range of species is involved including Ae. aegypti, Ae. africanus, Ae. luteocephalus, Ae. furcifer, Ae. taylori, Ae. dalzieli, Ae. metallicus, Ae. neoafricanus, Ae. centropunctatus, Ae. hirsutus, Anopheles domicola, An. funestus, An. coustani, Mansonia uniformis and Culex poicilipes [6,7,8,9,10,11]. Because of these known traditional vectors involved in the urban transmission cycle of CHIKV, entomological surveillance and outbreak response activities are usually biased toward Ae. aegypti or Ae. albopictus as the main vectors, leading to neglect of other species even where their densities and possible involvement in the transmission of this virus should be explored. For instance, during chikungunya outbreaks in Reunion and Comoros Islands, CHIKV was isolated from Culex P. quinquefasciatus and a pool of Culex spp. mosquitoes, respectively, although this was attributed to possible midgut infection and not investigated further [12, 13]. This was probably because species other than Ae. aegypti and Ae. albopictus have rarely been associated with CHIKV transmission.

Infection with CHIKV is characterized by a spectrum of clinical manifestations ranging from asymptomatic to a mild flu-like syndrome including fever, headache, fatigue, nausea, chills; severe arthralgia and, recently, mortality [14]. These symptoms start 4 to 7 days post-exposure, and most resolve within 2 weeks of the acute phase. However, joint pain can persist for months or years following the initial infection [15].

After laboratory confirmation of human clinical cases during the 2017–2018 outbreak in Mombasa, we investigated entomological risk factors for CHIKV transmission as a means of identifying opportunities for targeted control. We evaluated adult vector occurrence and abundance as well as larval indices and screened the mosquitoes for virus presence. Following the unexpected isolation of CHIKV from Cx. quinquefasciatus, we investigated further the ability of this species to amplify and transmit the virus.

Materials and methods

Study site description

Mombasa, the study site, is well described by Lutomiah et al. [5]. It is administratively divided into six sub-counties: Changamwe (mainland), Likoni (South coast), Kisauni, Jomvu and Nyali (north coast) and Mvita (Mombasa Island). While chikungunya cases were reported in all six sub-counties, we only conducted sampling in Changamwe (Miritini ward), Mvita (Tononoka ward) and Jomvu (Mikindani ward), which recorded the most cases (Fig. 1).

Fig. 1
figure 1

The map of Mombasa county showing outbreak sites where vector sampling was conducted

Weather data

The prevailing weather data for Mombasa during the outbreak were obtained from the Kenya Meteorological Department. Most of December 2017, except for the first week, and January and February 2018 were relatively dry with the lowest mean humidity (67.8%) and mean temperature of 27.8 °C. Subsequently, March to June 2018 were wet months with May receiving the highest amount of rainfall, > 260 mm, and recording the highest mean humidity of 79% (Fig. 2).

Fig. 2
figure 2

(Source: Kenya Meteorological Department, 2017–2018)

Rainfall, temperature and relative humidity for Mombasa county during the chikungunya outbreak

Epidemiological data

An outbreak of chikungunya fever was first reported in Mombasa county, endemic to both dengue and chikungunya, in mid-December 2017. This overlapped with the dengue outbreak that had persisted in the area since April of the same year. The first suspected human clinical cases were reported on 6 December 2017 and confirmed as CHIKV infection on 14 December 2017 at the Arbovirus and Viral Hemorrhagic Fevers (VHF) laboratory at KEMRI’s Centre for Virus Research (CVR). Apparently three of the initial six clinical samples were positive for both DENV and CHIKV IgM, thus indicating co-infections. In total, 19 CHIKV-positive cases were reported in the month of December 2017. The peak of the chikungunya outbreak was in January 2018, and between this month and May 2018, a further 140 positive samples were recorded, mostly in Mvita sub-county. Two of these were co-infected with CHIKV and DENV.

Collection and rearing of immature mosquitoes from containers in the environment

Entomological surveillance was conducted during the month of January 2018. This coincided with the period when the highest number (n = 100) of CHIKV-positive cases was reported. The immatures (larvae and pupae) of Ae. aegypti mosquitoes were collected indoors and outdoors over a 10-day sampling period in the areas where chikungunya cases were confirmed. Sampling followed the verbal/oral consent of the household head in the presence of the area headman or sub-chief. All indoor sampling was conducted by field workers accompanied by a household member while outdoor sampling was in the presence of a headman or sub-chief. Individual houses were georeferenced and the demographic data (number of households that spent the night in each house prior to the sampling) recorded. All water-holding containers occurring indoors and outdoors were inspected. Wet containers were scored as either negative (with no Ae. aegypti immatures) or positive (with at least one immature Ae. aegypti). All immatures from each positive container were collected using white ladles or Pasteur pipettes depending on the container type. For jerry cans, the water was first screened with the help of a flashlight, and if found positive, the water was poured through a fine sieve into a clean white basin. All the immatures were then picked from the sieve immersed in water in the basin using a Pasteur pipette. The type of breeding habitats and their locations were recorded. Live immatures from each container type were put in a whirl pack with information on container type, area/site, whether indoor or outdoor and date of collection and transported in a cool box to the field laboratory. Pupae were separated from larvae, put in pupal cups and then placed in mosquito-rearing cages to develop to adults. Larvae were placed in white enamel trays and fed on fish food (Tetramin®) until pupation. The emerging adults were cryo-preserved in liquid nitrogen (LN2) for transportation to the VHF laboratory for identification using standard taxonomic keys [16].

Collection and identification of adult mosquitoes

Adult mosquitoes were collected in the areas where chikungunya cases were confirmed using CDC light and BG-Sentinel (Biogent) traps (John W. Hock). Ten of each of the two trap types were randomly placed in pairs outdoors per day/night, and each was baited with 0.5 kg of dry ice held in igloos next to it [17]. Each pair of CDC and BG-Sentinel traps, approximately 5 m apart, was positioned between 40 and 60 m away from the next pair. In total, 200 traps were set over the 10-day sampling period from 7 to 16 January 2018. The CDC light traps were set at 17:00 and retrieved at 06:00 for 10 consecutive sampling days. The BG-Sentinel traps were set from 17:00 to 17:00 the following day after a 24-h cycle. However, dry ice and collection bags were replaced every morning to target day-biting mosquitoes. All trapped mosquitoes were taken to the site laboratory for sorting and preservation and later transported to the VHF laboratory for identification [16] and pooling (maximum 25/pool) by area/site, sex, species, blood feeding status and date of collection.

Mosquito processing and virus isolation

Mosquito pools were homogenized using minimum essential medium (MEM) homogenization media supplemented with 15% heat-inactivated fetal bovine serum (FBS) (Sigma-Aldrich), 2% l-glutamine (Sigma-Aldrich) and 2% antibiotic/antimycotic solution (Sigma-Aldrich) with 10,000 U penicillin, 10 mg streptomycin and 25 μg amphotericin B per milliliter. The homogenates were centrifuged at 12,000 rpm for 10 min and the supernatants and pellets stored separately. Fifty microliters each of the mosquito homogenates was inoculated per well in 24-well plates containing confluent monolayers of Vero cells (CCL-81™) (ATCC®), grown in MEM (Sigma-Aldrich) growth media, with Earle’s salts and reduced NaHCO3, supplemented with 10% FBS, 2% l-glutamine and 2% antibiotic/antimycotic solution. The inoculated plates were incubated at 37 °C for 1 h to allow for virus adsorption and maintained in MEM with 2% FBS, 2% l-glutamine and 2% antibiotic/antimycotic solution at 37 °C. The cells were observed for cytopathic effect (CPE) over a 14-day period and the virus harvested when 75% CPE was observed. Cultures presenting with CPE were passaged to confirm infection and harvested for downstream PCR analysis and sequencing.

Total RNA extraction and virus identification by reverse transcription PCR

Total RNA was extracted from the supernatants of cell cultures exhibiting CPE using the Qiagen kit following the manufacturer’s protocol. The RNA extract was transcribed into cDNA using the First Strand cDNA synthesis kit (Invitrogen) and random hexamers followed by PCR amplification with Amplitaq Gold PCR master mix (Applied Biosystems) using a panel of primers targeting flavivirus, orthobunyavirus and alphavirus genera. The samples which tested positive with the alphavirus primers were further analyzed using specific primers for CHIKV, ndumu (NDUV), babanki (BBKV) and sindbis (SINV) (Table 1). Amplification products were resolved in 1.5% agarose gel in Tris-borate-EDTA buffer stained with ethidium bromide.

Table 1 Sequences and target regions of the primers used in the identification of CHIKV isolates

Sequencing and phylogenetic analysis

PCR products of the five isolates were purified using the DNA Clean & Concentrator Kit (ZymoResearch, US), and the clean products were sequenced using the Sanger method (Macrogen, Korea). Chromatogram files of the three successfully sequenced isolates were imported into Chromas v2.6, and low-quality regions for each of the forward and reverse reads were trimmed independently and a consensus sequence generated based on the forward and reverse sequences for each of the samples. Sequences generated were compared to those available in Genbank using Blast and phylogenetic analysis was performed using a set of reference sequences obtained from Genbank. Sequences generated were submitted to Genbank under the following ccession numbers: MT992066, MT992067 and MT992068.

Vector competence studies

Isolation of virus from mosquitoes is not indicative of the mosquitoes being a vector. Therefore, following the isolation of CHIKV from Cx. quinquefasciatus mosquitoes, we conducted studies to determine the susceptibility and transmission potential of this species for the virus using an isolate from human serum during the outbreak for all the infection assays. The virus was passaged once in a confluent monolayer of Vero cells in a T-25 cell culture flask, grown in MEM growth media and quantified by plaque assay. The virus was tenfold serial diluted, and the first five dilutions inoculated in 6-well plates, 100 μl/well, containing confluent Vero monolayers as described [22]. The infected cells were maintained with 2.5% methylcellulose mixed with 2× MEM (GIBCO), further incubated at 37 °C for 4 days and then fixed for 1 h with 10% formalin, stained for 2 h with 0.5% crystal violet, washed and the plaques counted to determine plaque-forming units per milliliter (PFU/ml), using the following formula:

$$\frac{{\text{Number of plaques}}}{d \times V} = {\text{ PFU/ml}}$$

where d is the dilution factor and V is the volume of diluted virus added to the wells [23].

Mosquito rearing and identification

Immature Culex mosquitoes collected from Mombasa were reared in the insectary under standard conditions of 28 °C temperature and 80% relative humidity, with a 12:12-h (light:dark) photoperiod. Adult mosquitoes were inactivated at – 20 °C for 45 s and morphologically identified using the identification key [16] under a dissecting microscope to select Cx. quinquefasciatus for use in the study. A leg was also detached from each for subsequent molecular confirmation of species as described by Smith and Fonseca [24]. The mosquitoes were then fed on anaesthetized clean laboratory mice for 45 min to stimulate egg production and provided with an oviposition cup, and the eggs collected were hatched to first filial generation (F1) for experimental studies.

Oral infection of Cx. quinquefasciatus mosquitoes

Three-to-four-day-old Cx. quinquefasciatus females (F1) were deprived of glucose for 24 h. The mosquitoes were allowed to feed for 1 h in batches of 50–100 on CHIKV-infectious bloodmeal (log108.6 PFU/ml, determined after 1 h of feeding) in a Hemotek artificial feeding system maintained at 37 °C. Those engorged were transferred to 1-l plastic cages (15–30/cage) with screened tops and maintained on 10% glucose at standard insectary conditions and a 12:12 h light:dark (L:D) photoperiod. The non-engorged mosquitoes were killed by freezing and subsequently incinerated.

Manipulation of the infected mosquitoes for virus analysis

At 7, 14 and 21 days post-exposure (dpe), 30% of the exposed mosquitoes were picked and sucrose-starved for 16 h. The mosquitoes were cold anesthetized, the legs and wings of each individual mosquito removed using different forceps and pins to avoid cross contamination and placed in the same 1.5-ml microfuge tubes (Eppendorf). The bodies were placed on a sticky tape and proboscises individually inserted into capillary tubes each containing 10–20 µl homogenization media; they were allowed to expectorate for 30 min. The media containing salivary expectorate was expelled into a cryovial with 200 µl of MEM homogenization media. Bodies of mosquitoes, with heads still attached, were each placed in a labeled 1.5-ml microcentrifuge tube containing 1000 μl of MEM homogenization media and homogenized using a mini bead beater (BioSpec Products Inc., Bartlesville, OK, USA) with the aid of copper beads (BB-caliber airgun shot). The homogenates were centrifuged at 12,000 rpm for 10 min at 4 °C and the supernatants stored in 1.5-ml tubes containing 500 µl of MEM homogenization media.

Determination of infection and dissemination rates and transmission potential

To determine the infection rates, mosquito body homogenates were inoculated in 24-well plates containing confluent monolayers of Vero cells and incubated at 37 °C in 5% CO2 for 1 h to allow for virus adsorption. One milliliter of MEM maintenance media was added and the cultures incubated further for 7 days while observing for CPE. Abdominal supernatant of each positive body was serial diluted (tenfold), and 100 µl of each dilution was inoculated to each of ten wells of the 12-well plate with Vero cells grown in MEM. The remaining two wells were inoculated with homogenized male Cx. quinquefasciatus mosquito supernatants as negative controls. The plates were incubated at 37 °C for 1 h with frequent agitation every 15 min to allow for virus adsorption. The infected cell monolayers were then overlaid with 2.5% methylcellulose supplemented with 2% FBS, 2% l-glutamine and 2% antibiotic/antimycotic and further incubated at 37 °C. On day 4, plates were fixed for 1 h with 10% formalin, stained for 2 h with 0.5% crystal violet, washed in running tap water and dried overnight. The plates were observed on a light box, and plaques counted and used to calculate the abdominal viral titer and infection rates. For each positive abdomen, corresponding legs/wings were homogenized and their infection status and viral titers determined as described for the abdomens. This process was also repeated for salivary expectorates corresponding to positive legs to determine the transmission rate and the virus titer in the saliva.

Statistical analysis

Data on adult mosquito collection were analyzed using Simpson’s dominance index and the Shannon-Wiener diversity index. The data on the Ae. aegypti survey were analyzed and calculated in terms of four larval indices, which include the house index (HI), container index (CI), Breteau index (BI) and pupal index (PI) or pupae-per-household index (PHI) using the following formulas [25]:

$${\text{HI}} = \frac{{\text{Number of houses with immature mosquitoes}}}{{\text{Number of inspected houses}}} \times 100$$
$${\text{CI}} = \frac{{\text{Number of containers with immature mosquitoes}}}{{\text{Number of wet containers}}} \times 100$$
$${\text{BI}} = \frac{{\text{Number of containers with immature mosquitoes}}}{{\text{Number of inspected houses}}} \times {1}00$$
$${\text{PHI}}/{\text{PI}} = \frac{{\text{Number of pupae collected}}}{{\text{Total number of inspected houses}}} \times {1}00$$

The pupae-per-person index (PPI) was also calculated as the total number of pupae sampled divided by the total population of the inspected households [26]. The risk of transmission of CHIKV was estimated using the WHO criteria. In this case, in an area where the HI, CI and BI exceed 35, 20 and 50, respectively, the risk of Aedes-borne viruses is considered to be high; at BI between 5 and 50, the density of Ae. aegypti is considered to be sufficient to promote an outbreak; at HI, CI and BI of 4, 3 and 5, respectively, it is considered unlikely for Aedes-borne virus transmission to occur [27]. Descriptive analysis was used for the distribution of wet containers and Ae. aegypti immatures. Wet containers with any number of larvae or pupae were considered “positive containers” while houses with positive containers were considered as “positive houses.” Prevalence of containers was calculated by dividing the total number of container types by the total number of all containers sampled. Container productivity was calculated as a percentage using the formula: number of immatures that emerged from a particular container type/total immatures that emerged × 100. However, to avoid bias arising from one type of container being more prevalent and therefore being sampled than others, mean container productivity was calculated by dividing the total number of Ae. aegypti immatures from a particular container type by the total number of that container type that were sampled (number of immatures from a container type/total number of that container type). The 95% confidence interval for container productivity for each category of container was also calculated [28]. The container preferences of Ae. aegypti breeding were assessed by determining the breeding preference ratio (BPR) as described by Kumar et al. [29]. This was calculated by dividing the percentage of positive containers by the percentage of containers sampled. Minimum infection rate (MIR) for each species was calculated as the [number of positive pools/total specimens tested] × 1000. All data collected were analyzed using the statistical software package STATA 13.1 (StataCorp LP, TX, USA). The maximum likelihood (ML) tree of Kenyan CHIKV was generated using the Tamura-Nei substitution model and tested using 10,000 bootstraps in MEGA version 7.0.26.

Results

Aedes aegypti larval infestation levels in different domestic container types

In total, 1636 containers and 7 container types were inspected indoors and outdoors in Tononoka, with jerry cans being the most abundant (n = 1232; 75.17%), followed by drums (n = 249; 15.19%), buckets (n = 96; 5.86%) and basins (n = 41; 2.5%), while clay pots were the least abundant. A total of 128 containers (7.81%) were infested with at least one immature Ae. aegypti mosquito. Of these, 65.63% were drums followed by jerry cans (14.84%) while clay pots registered 0.78% of all positive containers. All the inspected containers produced a combined total of 528 immature Ae. aegypti mosquitoes, with drums being the most productive cumulatively (n = 317; 60%) followed by jerry cans (n = 74; 14%) and water tanks (n = 62; 11.7%) while clay pots were the least productive (n = 1; 0.2%). However, mean container productivity showed that water tanks were the most productive with a mean infestation of 6.2 Ae. aegypti mosquitoes per container, tire (3.1), drum (1.3) and basin (1.05). Among the sampled container types, BPR was highest for tires (11.36) followed by clay pots (6.50), water tanks (5.13) and drums (4.31) while least for buckets (0.93) (Table 2).

Table 2 Container types, positivity, productivity and BPR for immature Ae. aegypti mosquitoes

Level of household risk for CHIK transmission

A total of 547 people spent the night in the 131 houses inspected on the eves of the sampling, representing a mean population of 4.18/household. Forty-eight of the houses and 128 out of the 1637 inspected containers were positive for immature Ae. aegypti; thus, HI = 36.6, CI = 7.82 and BI = 97.71 were recorded during the outbreak. Seventy pupae were also collected, resulting in PHI/PI of 53.44%, while PPI was 12.80%. According to the WHO, thresholds for risk of outbreak/transmission for CI is > 3, HI > 4 and BI > 5 [27]. Therefore, based on our findings, the risk of chikungunya transmission in Mombasa was very high.

Diversity of adult mosquito species sampled during the outbreak

A total of 6899 adult mosquitoes were collected by BG-Sentinel and CDC light traps from the three sampled sites, Tononoka (n = 4733), Mikindani (n = 1638) and Miritini (= 528), in decreasing order of abundance. These were identified into 15 mosquito species and 2 species groups, including Cx. quinquefasciatus (n = 4492; 65.04%), which was the predominant species followed by Ae. vittatus (= 1137; 16.46%), Ae. aegypti (n = 911; 13.19%) and 14 other species (n = 367; 5.31%). The least sampled species included Ae. simpsoni, Anopheles gambiae, Cx. poicilipes, Eretmapodites chrysogaster and Mansonia uniformis. As expected, overall, BG-Sentinel traps collected more Ae. aegyti mosquitoes (n = 838) compared to the CDC light traps (n = 73) as well as Ae. vittatus (n = 1106 and n = 31, respectively). However, it was interesting to observe higher numbers of Cx. quinquefasciatus caught by BG-Sentinel traps (n = 4060) compared to CDC light traps (n = 432) (Table 3). Simpson’s dominance index and Shannon-Wiener diversity index of 0.4388 and 1.1942 were also respectively recorded for all the adult mosquitoes collected.

Table 3 Mosquito species and total collections by BG-Sentinel and CDC light traps

Chikungunya virus detection and isolation from sampled adult mosquitoes

From 524 pools of adult mosquitoes processed, CHIKV was isolated from 5, 2 each of female (MSA/S24/3013 and MSA/S24/3031) and male (MSA/S23/2481 and MSA/S24/3010) Cx. quinquefasciatus mosquitoes and 1 pool of female Ae. aegypti (MSA/S23/2444), all collected in Tononoka. These isolations represented minimum infection rates (MIRs) of 3.0 for female Ae. aegypti, 0.8 for female Cx. quinquefasciatus and 1.0 for male Cx. quinquefasciatus mosquitoes (Table 4).

Table 4 Mosquito species collected, CHIKV-positive pools and the MIR

Sequencing and phylogenetic analysis

Only three (MSA/S23/2481, MSA/S24/3010 and MSA/S24/3031) out of the five isolates sequenced returned successful sequences. A comparison of the obtained partial sequences showed great similarity with each other. Importantly, they clustered together and formed a clade with isolates previously collected between 2014 and 2018 from Kilifi, a county located in the Southeast Coast region and neighboring Mombasa. These isolates aligned within the ECSA lineage of CHIKV, a clade that also included one isolate from South Sudan (Fig. 3), which neighbors Kenya.

Fig. 3
figure 3

Maximum likelihood phylogenetic tree generated based on sequences belonging to the ECSA, Asian and West African lineages. Sequences obtained in this study are in red fonts

Susceptibility and transmission rates from vector competence studies

Forty-seven Cx. quinquefasciatus mosquitoes were exposed to infectious bloodmeals with viremia of 108.6 PFU/ml. Twenty-five percent of those sampled on day 7 post exposure (dpe) were infected, while 31.3% and 26.7% on 14 and 21 dpe respectively were also infected. The difference in infection rates between 7 and 21 dpe was not statistically significant, Fisher’s exact test = 1.000. Dissemination rates were 25%, 40% and 13.3% on 7, 14 and 21 dpe respectively, with dissemination titers ranging between 104.1 PFU/ml and 105.17 PFU/ml. All the mosquitoes with disseminated infection on 7 and 14 dpe transmitted the virus while there was no transmission on 21 dpe. Transmission titers ranged from 105.13 to 105.78 PFU/ml (Table 5).

Table 5 Infection, dissemination and transmission rates for Cx. quinquefasciatus mosquitoes after exposure to infectious bloodmeals with viremia of log 108.6 PFU/ml

Discussion

The global reemergence and spread of CHIKV is driven by the widely distributed primary vector Ae. aegypti [30]. There is also evidence that mutations adapt the virus for transmission by new mosquito vectors such as Ae. albopictus [31, 32]. Although whether mutations are involved in the transmission of CHIKV by Cx. quinquefasciatus has not been determined, control of CHIKV would be best facilitated through identification of breeding habitats of all vectors involved during outbreaks.

Fewer container types were inspected during this outbreak (7) than during the dengue outbreak (17) some years before in the same sites [5]. We attribute this to the dry conditions prevailing during the 2017–2018 sampling period compared to that of 2013–2014, which was during moderate rainfall. Most wet containers were found indoors, and jerry cans were the most frequently sampled (75.26%). However, drums were the most productive containers with 60% of Ae. aegypti immatures although they had a moderate BPR of 4.31. While tires had the highest BPR (11.36) and were the second most sampled, jerry cans had the lowest (0.2). The high percentage productivity for drums may be because of their wide openings, which allow easy access for gravid mosquitoes seeking oviposition sites. Most importantly, they also store water for longer, therefore allowing continuous breeding of mosquitoes. Like drums, tires and tanks also remain wet for longer, although their individual productivity during the sampling period was low. For the tires, this was because of the dry spell; hence, only a few of them had water that could support mosquito breeding. Although buckets and basins are used for temporary storage of water in kitchens and/or bathrooms, the water usage practice results in frequent depletion and replenishment, hence not providing sufficient time for mosquitoes to complete the breeding cycle. Traditionally, water tanks are difficult to sample from because of their depth [5], so the observed low numbers of Ae. aegypti immatures may not reflect their true productivity given that they can retain water for a long time, thus allowing for continuous breeding.

A community’s water storage practice impacts arboviral disease occurrence. Mombasa City authorities continue to face challenges in supplying reliable piped water to its residents. This has influenced the “container-water-storage” practice such that even during the dry season there are wet containers which allow for continuous breeding of Ae. aegypti. Hence, the observed larval indices of HI = 36.6%, CI = 7.82% and BI = 97.71% above the minimum thresholds (HI > 4, CI > 3 and BI > 5) for disease transmission [27] may explain the occurrence of chikungunya outbreaks in this region [3]. However, these indices are inadequate to measure transmission risk [33] due to lack of clear correlation with disease transmission [34]. Hence, currently, focus is on the pupal indices as an alternative for assessing the risk of transmission. During this investigation, the PHI/PI = 53.44% and PPI = 12.80% were comparatively higher than those observed in Honduras (0.25) and Brazil (0.15) [35, 36]. Therefore, it is becoming increasingly clear that transmission thresholds are dynamic and differ by geographic regions based on the complexity of factors that influence disease risk and must be determined independently for each region [37,38,39,40]. These thresholds have not been determined specifically for Kenya with respect to dengue and chikungunya transmissions. Now with the entry of Cx. quinquefasciatus as a competent vector of CHIKV, these indices will play a minimal role in determining risk of disease transmission, so other better means of predicting outbreaks are needed for guiding intervention measures.

Approximately 6900 adult mosquitoes were collected, 911 Ae. aegypti and more than 4400 Cx. quinquefasciatus. The Simpson’s dominance index of 0.4388 suggested that, at the time, Mombasa County had moderate diversity of mosquito species predominated by Cx. quinquefasciatus. The collection of significantly more adult Ae. vittatus than Ae. aegypti (P = 0.0001) was also interesting. During previous dengue outbreak investigations in Mombasa and routine surveillance activities in neighboring rural villages, only insignificant numbers were recorded [5, 41] suggesting that this is a recent invasion. This observation, and the sampling of four larvae of this species from drums indoors and tires outdoors, further confirms the changing patterns in their breeding behavior [42] and suggests that the main breeding habitat in Mombasa is unidentified. Traditionally, Ae. vittatus breeds in rock pools and tree holes [43] but also in diverse habitats including discarded containers and cement tanks, among others [42, 44, 45].

No CHIKV was isolated from Ae. vittatus; hence, its role in the outbreak remains unclear. However, previous experimental studies have shown that they are competent vectors [42, 46, 47]. This species also has an anthropophilic host selection pattern [41, 48, 49]. A combination of these factors would increase its vectorial capacity for CHIKV in urban cycles.

However, CHIKV was isolated from a pool of Ae. aegypti and four pools of male and female Cx. quinquefasciatus mosquitoes. Traditionally, investigations of chikungunya outbreaks result in few isolates from the primary vector, Ae. aegypti. For instance, during the Comoros outbreak investigation, CHIKV was isolated from just two pools of Ae. aegypti, despite the magnitude of the epidemic, and one pool of Culex spp., although this was attributed to midgut infection and not investigated further [13]. Before this, the virus had also been isolated from Cx. quinquefasciatus in Reunion Island but none from Ae. aegypti [12]. Some studies have attributed this to inefficient sampling methods for Ae. aegypti [50, 51], which are required in large numbers for arbovirus detection [52]. Whereas human landing catch is the most effective, it raises ethical concerns especially during outbreaks due to increased risk of exposure [53]. BG-Sentinel traps, though also considered effective [54], collected only a few (n = 333) female Ae. aegypti mosquitoes.

Isolation of CHIKV from Cx. quinquefasciatus males is the first observation. The main mode of transmitting arboviruses to vectors is through a horizontal process by feeding on viremic hosts. Hence, this observation significantly adds to existing knowledge about vertical transmission playing a role in the natural maintenance of CHIKV in the environment during inter-epidemic periods [55], besides the sylvatic cycle involving mosquitoes, primates and probably rodents. Vertical transmission of CHIKV was previously documented in field-collected and experimentally infected Ae. aegypti and Ae. albopictus mosquitoes [56,57,58,59,60,61].

The similarity of the current CHIKV isolates to those found circulating in this region from as early as 2014 [62] suggests that this strain was acquired by Cx. quinquefasciatus while feeding on viremic hosts.

Vector competence data incriminated this species as being able to transmit CHIKV. The infectious bloodmeal used had a titer of 108.6 PFU/ml, which is well above the 103.6–106.1 PFU/ml that has been used before to infect Ae. aegypti and Ae. albopictus mosquitoes [56] but within range (104–109 PFU/ml) of what has been observed in viremic patients [57]. Infection and dissemination rates of Cx. quinquefasciatus on 7, 14 and 21 dpe ranged from 25–31% and 13–40%, respectively, while transmission rates on day 7 and 14 dpe were 6.3% and 12.5%, respectively, based on all exposed mosquitoes regardless of infection status, and 100% for those with disseminated infection. This suggests moderate midgut infection (MIB) and escape barriers (MEB) but a weak salivary gland escape barrier (SGEB). However, experimental demonstration of transmission by the capillary method is less efficient than animal models, which are currently unavailable for CHIKV. These experiments confirmed that Cx. quinquefasciatus may be contributing to the perpetuation of CHIKV outbreaks that remain unevaluated in the coastal region. Previous experimental studies with Culex pipiens showed no evidence of CHIKV transmission [63]. The first laboratory evaluation of vector competence of Cx. quinquefasciatus yielded no infections either [8] while subsequent studies showed only midgut infection [64]. Therefore, data from this study add to the increasing spectrum of potential vectors of CHIKV that should be considered in entomological investigations and control strategies during outbreaks.

Other than predominantly feeding on avians [65, 66], Cx. quinquefasciatus are also opportunistic feeders and readily feed on humans in urban and peri-urban areas [67]. Therefore, in cities such as Mombasa where this competent species is abundant year-round, Cx. quinquefasciatus can contribute to sustained chikungunya outbreaks. Its coexistence with Ae. aegypti in the same geographic locales means that when recommending control measures, consideration should incorporate those that target Cx. quinquefasciatus.

Conclusion

We observed high larval indices, which may have significantly contributed to the chikungunya outbreak. The unexpected isolation of CHIKV from Cx. quinquefasciatus and the experimental confirmation of transmission of CHIKV point to its possible role in virus circulation in the coastal region. This implies that each outbreak should be approached with an open mind concerning possible incrimination of new vectors involved in the transmission for appropriate vector control measures to be instituted. The coexistence of Cx. quinquefasciatus and Ae. aegypti in this region will also likely further complicate the vector control processes considering that both species require different approaches. Therefore, these findings can inform future considerations of both species when planning vector control programs during chikungunya outbreaks and routine surveillance activities.

Availability of data and materials

All data generated or analysed during this study are included in this published article.

Abbreviations

BBKV:

Babanki virus

BG:

Biogent

BI:

Breteau index

BPR:

Breeding preference ratio

CDC:

Centers for Disease Control

CHIKV:

Chikungunya virus

CI:

Container index

CVR:

Centre for Virus Research

CPE:

Cytopathic effect

DENV:

Dengue virus

ECSA:

East, Central and Southern Africa

FBS:

Fetal bovine serum

HI:

House index

LN:

Liquid nitrogen

MEM:

Minimum essential medium

MIR:

Minimum infection rate

NDUV:

Ndumu virus

PCR:

Polymerase chain reaction

PFU:

Plaque-forming units

PHI:

Pupae-per-household index

PI:

Pupal index

PPI:

Pupae-per-person index

RH:

Relative humidity

SINV:

Sindbis virus

VHF:

Arbovirus and viral hemorrhagic fevers

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.

    Article  CAS  PubMed  Google Scholar 

  2. Sergon K, Njuguna C, Kalani R, Ofula V, Onyango C, Konongoi LS, et al. Seroprevalence of chikungunya virus (CHIKV) infection on Lamu Island, Kenya, October 2004. Am J Trop Med Hyg. 2008;78(2):333–7.

    Article  PubMed  Google Scholar 

  3. Konongoi SL, Nyunja A, Ofula V, Owaka S, Koka H, Koskei E, et al. Human and entomologic investigations of chikungunya outbreak in Mandera, Northeastern Kenya, 2016. PLoS ONE. 2018;13(10):e0205058. https://doi.org/10.1371/journal.pone.0205058.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Higgs S, Vanlandingham D. Chikungunya virus and its mosquito vectors. Vector-Borne Zoonotic Dis. 2015;15:231–40. https://doi.org/10.1089/vbz.2014.1745.

    Article  PubMed  Google Scholar 

  5. Lutomiah J, Barrera R, Makio A, Mutisya J, Koka H, Owaka S, et al. Dengue outbreak in Mombasa City, Kenya, 2013–2014: entomologic investigations. PLoS Negl Trop Dis. 2016;10(10):e0004981. https://doi.org/10.1371/journal.pntd.0004981.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Diallo M, Thonnon J, Traore-lamizana M, Fontenille D. Vectors of chikungunya virus in Senegal: current data and transmission cycles. Am J Trop Med Hyg. 1999;60:281–6.

    Article  CAS  PubMed  Google Scholar 

  7. Jupp PG, McIntosh BM. Aedes furcifer and other mosquitoes as vectors of chikungunya virus at Mica, northeastern Transvaal, South Africa. J Am Mosq Control Assoc. 1990;6:415–20.

    CAS  PubMed  Google Scholar 

  8. Jupp PG, McIntosh BM, Dos Santos I, DeMoor P. Laboratory vector studies on six mosquito and one tick species with chikungunya virus. Trans R Soc Trop Med Hyg. 1981;75:15–9. https://doi.org/10.1016/0035-9203(81)90005-5.

    Article  CAS  PubMed  Google Scholar 

  9. McIntosh BM, Paterson HE, McGillivray G, Desousa J. Further studies on the chikungunya outbreak in Rhodesia. I. Mosquitoes, wild primates and birds in relation to the epidemic. Ann Trop Med Parasitol. 1964;58:45–51. https://doi.org/10.1080/00034983.1964.11686213.

    Article  CAS  PubMed  Google Scholar 

  10. Weinbren MP, Haddow AJ, Williams MC. The occurrence of chikungunya virus in Uganda. I. Isolation from mosquitoes. Trans R Soc Trop Med Hyg. 1958;52:253–7. https://doi.org/10.1016/0035-9203(58)90084-1.

    Article  CAS  PubMed  Google Scholar 

  11. Diallo D, Sall AA, Buenemann M, Chen R, Faye O, Diagne CT, et al. Landscape ecology of sylvatic chikungunya virus and mosquito vectors in Southeastern Senegal. PLoS Negl Trop Dis. 2012;6(6):e1649. https://doi.org/10.1371/journal.pntd.0001649.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Bessaud M, Peyrefitte CN, Pastorino BAM, Tock F, Merle O, Colpart JJ, Dehecq JS, Girod R, Jaffar-Bandjee MC, Glass PJ, et al. Chikungunya virus strains, Reunion Island outbreak. Emerg Infect Dis. 2006;12:1604–6.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Sang RC, Ahmed O, Faye O, Kelly CLH, Yahaya AA, Mmadi I, et al. Entomologic investigations of a chikungunya virus epidemic in the Union of the Comoros, 2005. Am J Trop Med Hyg. 2008;78(1):77–82.

    Article  CAS  PubMed  Google Scholar 

  14. Josseran L, Paquet C, Zehgnoun A, Caillere N, Le Tertre A, Solet JL, et al. Chikungunya disease outbreak, Reunion Island. Emerg Infect Dis. 2006;12:1994–5.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Schwartz O, Albert ML. Biology and pathogenesis of chikungunya virus. Nat Rev Microbiol. 2010;8:491–500.

    Article  CAS  PubMed  Google Scholar 

  16. Edwards FW. Mosquitoes of the Ethiopian region. Iii.-culicine adults and pupae. Mosquitoes of the Ethiopian region. III.-Culicine Adults and Pupae. 1941.

  17. McNelly JR. The CDC trap as a special monitoring tool. In: Proceedings of the seventy-sixth annual meeting of the New Jersey Mosquito Control Association, Inc. 1989; pp. 26–33.

  18. Bryant JE, Crabtree MB, Nam VSN, Yen T, Duc HM, Miller BR. Isolation of arboviruses from mosquitoes collected in northern Vietnam. Am J Trop Med Hyg. 2005;73:47–53.

    Article  Google Scholar 

  19. Kuno G, Mitchell CJ, Chang GJJ, Smith GC. Detecting bunyaviruses of the Bunyamwera and California serogroups by a PCR technique. J Clin Microbiol. 1996;34:1184–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Eshoo MW, Whitehouse CA, Zoll ST, Massire C, Pennella TD, Blyn LB, et al. Direct broad-range detection of alphaviruses in mosquito extracts. Virology. 2007;368:286–95.

    Article  CAS  PubMed  Google Scholar 

  21. GenBank Accession AF369024.

  22. Gargan TP II, Bailey CL, Higbee GA, Gad A, El Said S. The effect of laboratory colonization on the vector-pathogen interactions of Egyptian Culex pipiens and Rift Valley fever virus. Am J Trop Med Hyg. 1983;32(5):1154–63.

    Article  PubMed  Google Scholar 

  23. Reiskind MH, Pesko K, Westbrook CJ, Mores CN. Susceptibility of Florida mosquitoes to infection with chikungunya virus. Am J Trop Med Hyg. 2008;78(3):422–5.

    Article  PubMed  Google Scholar 

  24. Smith JL, Fonseca DM. Rapid assays for identification of members of the Culex (Culex) pipiens complex, their hybrids, and other sibling species (Diptera: Culicidae). Am J Trop Med Hyg. 2004;70(4):339–45.

    Article  CAS  PubMed  Google Scholar 

  25. Pan American Health Organization. Dengue and dengue haemorrhagic fever in the Americas: guidelines for prevention and control. Washington DC: Pan American Health Organization; 1994.

    Google Scholar 

  26. Jiménez-Alejo A, Morales-Pérez A, Nava-Aguilera E, Flores-Moreno M, Apreza-Aguilar S, Carranza-Alcaraz W, et al. Pupal productivity in rainy and dry seasons: findings from the impact survey of a randomised controlled trial of dengue prevention in Guerrero, Mexico. BMC Public Health. 2017;2017(17 Suppl 1):428. https://doi.org/10.1186/s12889-017-4294-8.

    Article  Google Scholar 

  27. Connor ME, Monroe WM. Stegomyia indices and their value in yellow fever control. Am J Trop Med Hyg. 1923;3:9–19.

    Article  Google Scholar 

  28. Hammond SN, Gordon AL, Lugo EC, Moreno G, Kuan GM, Lopez MM, et al. Characterization of Aedes aegypti (Diptera: Culcidae) production sites in urban Nicaragua. J Med Entomol. 2007;44:851–60.

    Article  PubMed  Google Scholar 

  29. Kumar R, Kamal S, Patnaik SK, Sharma RC. Breeding habitats and larval indices of Aedes aegypti (L) in residential areas of Rajahmundry town Andhra Pradesh. J Commun Dis. 2002;34:50–8.

    PubMed  Google Scholar 

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

    Google Scholar 

  31. Tsetsarkin KA, Vanlandingham DL, McGee CE, Higgs S. A single mutation in chikungunya virus affects vector specificity and epidemic potential. PLoS Pathog. 2007;3:e201.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Vazeille M, Moutailler S, Coudrier D, Rousseaux C, Khun H, 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. https://doi.org/10.1371/journal.pone.0001168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Focks DA, Chadee DD. Pupal survey: an epidemiologically significant surveillance method for Aedes aegypti: an example using data from Trinidad. Am J Trop Med Hyg. 1997;56:159–67.

    Article  CAS  PubMed  Google Scholar 

  34. Weaver SC, Morrison AC. Aedes aegypti density and the risk of dengue virus transmission. In: Takken W, Scott TW, editors. Ecological aspects for application of genetically modified mosquitoes. Dordrecht: FRONTIS; 2003. p. 187–206.

    Google Scholar 

  35. 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.

    Article  CAS  PubMed  Google Scholar 

  36. Pilger D, Lenhart A, Manrique-Saide P, Siqueira JB, Da Rocha WT, et al. Is routine dengue vector surveillance in central Brazil able to accurately monitor the Aedesaegypti population? Results from a pupal productivity survey. Trop Med Int Health. 2011;16:1143–50.

    Article  PubMed  Google Scholar 

  37. Kuno GG. Review of the factors modulating dengue transmission. Epidemiol Rev. 1995;17:321–35.

    Article  CAS  PubMed  Google Scholar 

  38. Reich NG, Shrestha S, King AA, Rohani P, Lessler J, et al. Interactions between serotypes of dengue highlight epidemiological impact of cross-immunity. J R Soc Interface. 2013;10:86. https://doi.org/10.3201/eid1302.060539.

    Article  Google Scholar 

  39. Rabaa MA, Simmons CP, Fox A, Le MQ, Nguyen TTT, et al. Dengue virus in sub-tropical northern and central Viet Nam: population immunity and climate shape patterns of viral invasion and maintenance. PLoS Negl Trop Dis. 2013;7(12):e2581. https://doi.org/10.1371/journal.pntd.0002581.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Sommerfeld J, Kroeger A. Eco-bio-social research on dengue in Asia: a multicountry study on ecosystem and community-based approaches for the control of dengue vectors in urban and peri-urban Asia. Pathog Glob Health. 2012;106:428–35.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Lutomiah J, Bast J, Clark J, Richardson J, Yalwala S, Oullo D, et al. Abundance, diversity, and distribution of mosquito vectors in selected ecological regions of Kenya: public health implications. J Vector Ecol. 2013;38:134–42.

    Article  PubMed  Google Scholar 

  42. Diallo D, Diagne CT, Hanley KA, Sall AA, Buenemann M, Ba Y, Dia I, Weaver SC, Diallo M. Larval ecology of mosquitoes in sylvatic arbovirus foci in southeastern Senegal. Parasites Vectors. 2012;5:286. https://doi.org/10.1186/1756-3305-5-286.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Service MW. Studies on the biology and taxonomy of Aedes (Stegomya) vittatus (Bigot) (Diptera: Culicidae) in northern Nigeria. Trans R Soc Entomol Soc Lond. 1970;122:101–43.

    Article  Google Scholar 

  44. Adebote DA, Oniye SJ, Muhammed YA. Studies on mosquitoes breeding in rock pools on inselbergs around Zaria, northern Nigeria. J Vector Borne Dis. 2008;45(1):21–8.

    PubMed  Google Scholar 

  45. Tewari SC, Thenmozhi V, Katholi CR, Manavalan R, Munirathinam A, Gajanana A. Dengue vector prevalence and virus infection in a rural area in south India. Trop Med Int Health. 2004;9(4):499–507.

    Article  CAS  PubMed  Google Scholar 

  46. Mulwa F, Lutomiah J, Chepkorir E, Okello S, Eyase F, Tigoi C, et al. Vector competence of Aedes bromeliae and Aedes vitattus mosquito populations from Kenya for chikungunya virus. PLoS Negl Trop Dis. 2018;12(10):e0006746. https://doi.org/10.1371/journal.pntd.0006746.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Diagne CT, Faye O, Guerbois M, Knight R, Diallo D, Ba Y, Dia I, Faye O, Weaver SC, Sall AA, Diallo M. Vector competence of Aedes aegypti and Aedes vittatus (Diptera: Culicidae) from Senegal and Cape Verde archipelago for West African lineages of chikungunya virus. Am J Trop Med Hyg. 2014;91(3):635–41.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Wilson JJ, Sevarkodiyone SP. Host preference of blood feeding mosquitoes in rural areas of Southern Tamil Nadu, India. Acad J Entomol. 2015;8(2):80–3.

    Google Scholar 

  49. Lee VH, Moore DL. Vectors of the 1969 yellow fever epidemic on the Jos Plateau, Nigeria. Bull World Health Organ. 1972;46(5):669–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Diallo M, Dia I, Diallo D, Diagne CT, Ba Y, Yactayo S. Perspectives and challenges in entomological risk assessment and vector control of chikungunya. J Infect Dis. 2016;214(suppl 5):S459–65. https://doi.org/10.1093/infdis/jiw397.

    Article  PubMed  Google Scholar 

  51. Sivagnaname N, Gunasekaran K. Need for an efficient adult trap for the surveillance of dengue vectors. Indian J Med Res. 2012;136(5):739–49.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Gu W, Novak RJ. Short report: detection probability of arbovirus infection in mosquito populations. Am J Trop Med Hyg. 2004;71(5):636–8.

    Article  PubMed  Google Scholar 

  53. Williams CR, Long SA, Russell RC, Ritchie SA. Field efficacy of the BG-Sentinel compared with CDC Backpack Aspirators and CO2-baited EVS traps for collection of adult Aedes aegypti in Cairns, Queensland, Australia. J Am Mosq Control Assoc. 2006;22(2):296–300.

    Article  PubMed  Google Scholar 

  54. Harwood JF, Arimoto H, Nunn P, Richardson AG, Obenauer PJ. Assessing carbon dioxide and synthetic lure-baited traps for dengue and chikungunya vector surveillance (3). J Am Mosq Control Assoc. 2015;31(3):242–7.

    Article  PubMed  Google Scholar 

  55. Hardy JL, Houk EJ, Kramer LD, Reeves WC. Intrinsic factors affecting vector competence of mosquitoes for arboviruses. Annu Rev Entomol. 1983;28:229–62. https://doi.org/10.1146/annurev.en.28.010183.001305.

    Article  CAS  PubMed  Google Scholar 

  56. Thavara U, Tawatsin A, Pengsakul T, Bhakdeenuan P, Chanama S, Anantapreecha S, et al. Outbreak of Chikungunya fever in Thailand and virus detection in field population of vector mosquitoes, Aedes aegypti (L.) and Aedes albopictus Skuse (Diptera: Culicidae). Southeast Asian J Trop Med Public Health. 2009;40:951–62.

    PubMed  Google Scholar 

  57. Delatte H, Paupy C, Dehecq JS, Thiria J, Failloux AB, Fontenille D. Aedes albopictus, vector of chikungunya and dengue viruses in Reunion Island: biology and control. Parasite. 2008;15:3–13.

    Article  CAS  PubMed  Google Scholar 

  58. Ratsitorahina M, Harisoa J, Ratovonjato J, Biacabe S, Reynes JM, Zeller H, et al. Outbreak of dengue and chikungunya fevers, Toamasina, Madagascar, 2006. Emerg Infect Dis. 2008;14:1135–7.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Honório AN, Wiggins K, Eastmond B, Câmara CPD, Alto BW. Experimental vertical transmission of chikungunya virus by Brazilian and Florida Aedes albopictus populations. Viruses. 2019;11:353. https://doi.org/10.3390/v11040353.

    Article  CAS  PubMed Central  Google Scholar 

  60. Chompoosri J, Thavara U, Tawatsin A, Boonserm R, Phumee A, Sangkitporn S, Siriyasatien P. Vertical transmission of Indian Ocean Lineage of chikungunya virus in Aedes aegypti and Aedes albopictus mosquitoes. Parasit Vectors. 2016;9:227.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. Agarwal A, Dash PK, Singh AK, Sharma S, Gopalan N, Rao PV, Parida MM, Reiter P. Evidence of experimental vertical transmission of emerging novel ECSA genotype of chikungunya virus in Aedes aegypti. PLoS Negl Trop Dis. 2014;8:e2990. https://doi.org/10.1371/journal.pntd.0002990 (eCollection 2014).

    Article  PubMed  PubMed Central  Google Scholar 

  62. Nyamwaya DK, Otiende M, Omuoyo DO, Githinji G, Karanja HK, Gitonga JN, et al. Endemic chikungunya fever in Kenyan. medRxiv preprint. 2020. https://doi.org/10.1101/2020.07.22.20116707.

  63. Brugman A, Hernández-Triana LM, Medlock JM, Fooks AR, Carpenter S, Johnson N. The role of Culexpipiens L. (Diptera: Culicidae) in virus transmission in Europe. Int J Environ Res Public Health. 2018;15(2):389. https://doi.org/10.3390/ijerph15020389.

    Article  PubMed Central  Google Scholar 

  64. Richards SL, Anderson SL, Smartt CT. Vector competence of Florida mosquitoes for chikungunya virus. J Vector Ecol. 2010;35(2):439–43. https://doi.org/10.1111/j.1948-7134.2010.00105.x.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Kay BH, Boreham PF, Fanning ID. Host-feeding patterns of Culexannulirostris and other mosquitoes (Diptera: Culicidae) at Charleville, southwestern Queensland, Australia. J Med Entomol. 1985;22:529–35.

    Article  CAS  PubMed  Google Scholar 

  66. Jansen CC, Webb CE, Graham GC, Craig CB, 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:849–57.

    Article  PubMed  Google Scholar 

  67. Janssen N, Fernandez-Salas I, Gonzalez EED, Gaytan-Burns A, Medina-delaGarza CE, Sanchez-Casas RM, et al. Mammalophilic feeding behaviour of Culex quinquefasciatus mosquitoes collected in the cities of Chetumal and Cancun, Yucatan Peninsula, Mexico. Trop Med Int Health. 2015;20:1488–91.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dunstone Beti, John Gachoya and Reuben Lugalia (KEMRI) for assisting in field collection and identification of the mosquitoes.

Disclaimer

The views expressed are those of the authors and should not be construed to represent the position of KEMRI.

Funding

This study was made possible through the financial support by the Government of Kenya to the Kenya Medical Research Institute (KEMRI) to respond to and investigate the chikungunya outbreak in Mombasa County.

Author information

Authors and Affiliations

Authors

Contributions

JL, FE and RS designed and conceptualized the study. JL, FM, JM, SO, VO and SK coordinated the fieldwork and conducted mosquito sampling. HK, EK, EC, VO and SK conducted cell culture work and isolated the viruses. AN and SL carried out the molecular laboratory work and characterized the virus isolates. JL, FM and JM conducted vector competence studies. JL and RS performed the analysis and interpreted the results. JL drafted the manuscript. RS and FE coordinated the entire study. All authors reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Joel Lutomiah.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lutomiah, J., Mulwa, F., Mutisya, J. et al. Probable contribution of Culex quinquefasciatus mosquitoes to the circulation of chikungunya virus during an outbreak in Mombasa County, Kenya, 2017–2018. Parasites Vectors 14, 138 (2021). https://doi.org/10.1186/s13071-021-04632-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-021-04632-6

Keywords