Endemic Chikungunya Fever in Kenyan Children
============================================

* Doris K. Nyamwaya
* Mark Otiende
* Donwilliams O. Omuoyo
* George Githinji
* Henry K. Karanja
* John N. Gitonga
* Zaydah de Laurent
* James R. Otieno
* Rosemary Sang
* Everlyn Kamau
* Stanley Cheruiyot
* Edward Otieno
* Charles N. Agoti
* Philip Bejon
* Samuel M. Thumbi
* George M. Warimwe

## ABSTRACT

**Background** Chikungunya virus (CHIKV) was first identified in Tanzania in 1952. Several epidemics including East Africa are described, but there are no descriptions of longitudinal surveillance of endemic disease. Here, we estimate the incidence of CHIKV and describe viral phylogeny in coastal Kenya.

**Methods** Over a 5-year period (2014-2018), 11,708 febrile illnesses in 5,569 children visiting two primary healthcare facilities linked to a demographic surveillance system in coastal Kenya were recorded and blood samples obtained. Reverse-transcriptase PCR was used to identify CHIKF cases in 3,500 children randomly selected from the 5,569 children.

**Results** We found CHIKF to be endemic in this setting, associated with 12.7% (95% CI 11.60, 13.80) of all febrile presentations to primary healthcare. The prevalence of CHIKV infections among asymptomatic children in a community survey was 0.7% (95% CI 0.22, 2.12). CHIKF incidence among children <1 year of age was 1703 cases/100,000-person years and 46 cases/100,000-person years among children aged ≥ 10 years. Recurrent CHIKF episodes, associated with fever and viraemia, were observed among 19 of 170 children with multiple febrile episodes during the study period and confirmed by genome sequencing. All sequenced viral genomes mapped to the ECSA genotype albeit distinct from CHIKV strains associated with the 2004 East African epidemic.

**Conclusions** CHIKF may be a substantial public health burden in primary healthcare on the East African coast outside epidemic years, and recurrent infections are common.

## BACKGROUND

Chikungunya fever (CHIKF) is a zoonotic mosquito-borne febrile illness characterised by acute, often chronic, debilitating polyarthralgia and polyarthritis that can last for months to years [1-3]. The disease is caused by chikungunya virus (CHIKV), a positive sense RNA virus in the family *Togaviridae* that was first isolated from a febrile patient in Tanzania in 1952 [4]. CHIKV transmission between humans is mainly mediated by the geographically widespread *Aedes aegypti* and *Ae. albopictus* mosquitoes, with geographic spread and spillover from sylvatic transmission cycles in monkeys thought to account for the periodic re-emergence of human disease outbreaks [5-8].

The CHIKV genome is approximately 11.8kb in length and encodes four non-structural proteins (nsP1 to nsP4) required for viral replication, three major structural proteins (the capsid protein, and envelope proteins E1 and E2), and two small polypeptides (6K/TF and E3) [9]. Exposure to CHIKV results in acquisition of protective virus neutralising antibodies that target the E2 protein, and is the basis for ongoing efforts to develop CHIKF vaccines [10, 11]. Phylogenetic studies have defined three CHIKV genotypes, namely East, Central and Southern Africa (ECSA), West Africa, and the Asian genotype [12]. In 2004, CHIKV re-emerged in coastal Kenya causing one of the largest epidemics on record, affecting millions of people as it spread along the Indian Ocean islands, India, southeast Asia and Europe [8, 13, 14]. This epidemic was associated with emergence of the Indian Ocean Lineage (IOL), which mapped within the ECSA genotype and included viruses with adaptive mutations in the E1 protein that increased their transmissibility by *Aedes albopictus* mosquitoes [15, 16].

CHIKF epidemics are now frequently reported globally [8, 17], but despite its original discovery and later re-emergence in East Africa, very little is known regarding inter-epidemic CHIKV exposure in the region. CHIKF cases have previously been detected in children in Kenya and Tanzania during inter-epidemic periods [18-20], which together with the high anti-CHIKV antibody seroprevalence observed in these settings suggest endemic CHIKV transmission [21, 22]. To address these knowledge gaps, we conducted a primary healthcare-based study linked to demographic surveillance and community survey to estimate the prevalence and incidence of CHIKV infections among children in Kilifi, coastal Kenya. Our study period, 2014 to 2018, included the most recent CHIKF epidemic year, 2016 [23], allowing an unprecedented assessment of exposure risk and disease burden before, during and after a CHIKF epidemic in this setting.

## METHODS

### Study population

We conducted a population-based observational cohort study between March 2014 and October 2018 at two dispensaries, Ngerenya and Pingilikani, based within the Kilifi Health Demographic Surveillance System (KHDSS) [24]. The KHDSS covers an area of 891km2 with approximately 290,000 residents that are enumerated every four months at the household level [24]. Ngerenya and Pingilikani have different rates of mosquito-borne pathogen exposure as inferred from previous malaria studies; Ngerenya and environs experiences low intensity malaria transmission, whereas malaria transmission around Pingilikani is moderate [25]. However, no studies on arbovirus exposure had been conducted previously. Our approach was to screen children presenting with fever (axillary temperature of ≥37.5°C) at the dispensaries for CHIKV infection by reverse-transcriptase polymerase chain reaction (RT-PCR). Children presenting at the dispensaries were linked to the KHDSS data using a unique identifier, allowing us to estimate incidence in the population, and to determine the influence of various demographic variables on the incidence estimates. Ethical approval was provided by the Kenya Medical Research Institute Scientific and Ethics Review Unit (SSC Nos. 3296, 2617 and 3149).

### CHIKV RT-PCR assays

CHIKV infection was defined as detection of CHIKV viral RNA by RT-PCR using a published primer-probe set targeting the CHIKV non-structural protein 1 (nsP1) region[26], namely CHIKV 874 (5’-AAAGGGCAAACTCAGCTTCAC-3’), CHIKV 961 (5’-GCCTGGGCTCATCGTTATTC-3’) and CHIKV 899-FAM (5’-*FAM*-CGCTGTGATACAGTGGTTTCGTGTG-*TAMRA*-3’). Briefly, 100µl of finger-prick blood samples collected at presentation to the dispensaries, and stored at −80°C, was used for total RNA isolation using TRIzol™ Reagent (ThermoFisher) as per manufacturer instructions. Presence of CHIKV viral RNA was then determined using the QuantiFast RT-PCR Kit (Qiagen) in a 25µl reaction with primers and probe at a final concentration of 100nM and 20nM, respectively, and 5µl total RNA template. RT-PCR assays were done on a 7500 Real-Time PCR System (Applied Biosystems) with cycling conditions as follows: 50°C for 20 minutes, 95°C for 15 minutes, followed by 45 cycles of 94°C for 15 seconds and 60°C for 1 minute. A second confirmatory screen was performed on some of the RT-PCR positive samples using an assay targeting the CHIKV nsP4 region as previously described[27]. For both assays, a positive result was defined as a cycle threshold (Ct) value of <40. Viral RNA from a cultured CHIKV isolate obtained from a febrile Kenyan patient was used as a positive control, while RT-PCR mastermix without template was used as a negative control.

### CHIKV genome sequencing and analysis

RT-PCR positive CHIKV samples were sequenced using the Nanopore MinION technology following PCR amplification using methods described in the PrimalSeq approach1 7. Full details of our sequencing and bioinformatics workflow can be found in the Supplementary Material. In brief, CHIKV-specific multiplex primers were designed using the Primal Scheme software and used for pre-enrichment of CHIKV viral RNA in the samples using a multiplex PCR method. Barcodes and sequencing adapters were ligated into each sample, after which all samples were pooled into a single reaction tube and the library loaded on to a MinION sequencing device for sequencing. The sequences generated in this study (N=10) were deposited in GenBank, accession numbers: [MT526796](http://medrxiv.org/lookup/external-ref?link\_type=GEN&access_num=MT526796&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom)-[MT526807](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=MT526807&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom). MUSCLE was used to align the Kilifi CHIKV genome sequences with those available in GenBank from other geographical locations (see Supplementary Material). Maximum Likelihood phylogenies were reconstructed from the alignment using RAxML[28] using the GTR substitution model with 4 gamma categories (GTR+G4) and visualized in FigTree v1.4.4.

### Sampling and statistical analyses

We selected a random sample of 3500 out of 5,669 children who visited the two dispensaries with febrile illness. Associations between available demographic variables and CHIKV RT-PCR positivity were assessed using odds ratios estimated using logistic regression models. CHIKF incidence within the dispensary catchment area, defined as radius of 5km from each dispensary, was expressed per 100,000 person-years observed and incidence rate ratios (IRRs) comparing various demographic variables determined by Poisson regression. Univariate analyses were performed using Chi2 test (categorical variables) or the Mann-Whitney U test (continuous variables). All statistical analyses were carried out in Stata™ version 15 with a two-sided p value <0.05 used for statistical significance.

## RESULTS

### CHIKV infections are common among children in coastal Kenya

Between March 2014 and October 2018, there were 29,819 visits by children aged up to 15 years to the two dispensaries (6,746 in Ngerenya, 23,073 in Pingilikani), of which 13,696 were febrile (Figure 1). Of the 13,696 fevers, 1,256 lacked a unique patient identifier and 732 samples were missing. After these exclusions, there were 11,708 febrile visits from a total of 5,569 children (median 1 febrile visit per child during the study period, IQR 1 - 2) eligible for CHIKV RT-PCR analysis (Figure 1).

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/26/2020.07.22.20116707/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2020/07/26/2020.07.22.20116707/F1)

Figure 1: Study population.
The flow of study participants, reasons for exclusion and results from CHIK RT-PCR screening are shown.

A sample of 3,500 febrile children out of the 5,569 across the two dispensaries were randomly selected and samples from their first (index) febrile episode during the study period screened for CHIKV infection (Figure 1). The median age of these 3,500 children was 3.1 years (IQR 1.3 – 6.4), with 1,701 being resident in Ngerenya and 1,799 in Pingilikani, respectively (Figure 1). Our study period, 2014 to 2018, was at least 10 years since the 2004 CHIKF epidemic occurred in Kenya, with ∼90% our study population having been born after that epidemic.

Of the 3,500 children 443 were RT-PCR positive (12.7%, 95% CI 11.60, 13.80; Table 1), with RT-PCR positivity being twice as frequent in Ngerenya (16.7%, 95% CI 15.00, 18.54) than in Pingilikani (8.8%, 95% CI 7.61, 10.24). Most of the 443 CHIKF cases had a clinical diagnosis of upper respiratory tract infection (URTI) but in no instance was a diagnosis of CHIKF assigned (Table 1). The treatments prescribed for the management of CHIKF cases were comparable to those for other causes of fever at the dispensaries, but fewer antimalarials were prescribed for CHIKF cases (Table 1). CHIKF case prevalence was highest in 2016 (23.1%), when an epidemic was reported in Kenya [23], and was significantly lower in the pre- and post-epidemic years where it ranged between 5% and 8% (Table 1). In contrast, the prevalence of asymptomatic CHIKV infections estimated by cross-sectional sampling of a different group of 435 children in the same study location in 2016 was 0.7% (95% CI 0.22, 2.12); this yielded a clinical-to-asymptomatic CHIKV infection ratio of 18:1 during the epidemic year.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/07/26/2020.07.22.20116707/T1)

Table 1. Demographic characteristics of study participants

### Incidence of CHIKF in the community

We calculated the incidence of CHIKF restricting our analysis to children aged up to 15 years whose residence was within the dispensary catchment area. This represented a total of 139,840 person-years of observation (pyo) over the study period at both locations. For the numerator (CHIKF cases), we applied the observed CHIKV RT-PCR positivity rates (Table 1) to the total number of febrile children presenting to the respective dispensary on the assumption that the case positive rate among the 3,500 we randomly selected would hold true for the remaining 2,069 samples.

The incidence of CHIKF during the study period was 340 cases/100,000 pyo (95% CI 310, 373; Table 2). However, an inverse relationship was evident between CHIKF incidence and age, consistent with acquisition of protective immunity against disease due to ongoing CHIKV exposure in this setting (Table 2). CHIKF incidence was significantly higher in Ngerenya as compared to Pingilikani; furthermore, while an increase in incidence was observed in Ngerenya during the 2016 CHIKF epidemic, such an increase was not observed in Pingilikani where CHIKF incidence had been gradually declining since 2014 (Table 2). Marked seasonality in CHIKF incidence was observed in Pingilikani but not Ngerenya (Table 2) suggesting potential differences in ecological risk factors of disease. A negative correlation was observed between CHIKF incidence and residential distance from the dispensary at both locations, which likely reflected access to care (Table 2) as has been previously reported for malaria incidence [29].

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/07/26/2020.07.22.20116707/T2)

Table 2. Incidence of CHIKF in the study population

### Recurrent CHIKF episodes

We determined whether recurrent CHIKF episodes occurred in our study population. To do this we identified all 443 children with first or only episodes of CHIKF and screened all subsequent samples collected during febrile episodes during the period 2014-2018. We defined recurrent episodes as CHIKV RT-PCR positivity during these subsequent febrile dispensary visits. Of the 443 index cases, 170 presented to the dispensary with fever on at least one other occasion during the study period, contributing a total of 320 subsequent febrile dispensary visits with varying durations since the corresponding index CHIKF case (Figure 1). Of the 170 children, 19 (11.2%, 95% CI 6.86, 16.90) were CHIKV RT-PCR positive on at least one subsequent febrile episode (Figures 1 and 2). The duration between the index and the subsequent RT-PCR positive febrile episodes ranged from 2 to 43 months (Figure 2, and Supplementary Material Table S2). Eighteen of the 19 CHIKV RT-PCR positive cases were detected during or after the 2016 CHIKF epidemic year. Neither age, sex nor geographic location were statistically associated with recurrent episodes of CHIKF (Chi2 test p>0.05 for all) and the RT-PCR Ct values and clinical diagnoses were comparable between the index case and corresponding subsequent timepoints (Supplementary Material Figure S2 and Table S2).

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/26/2020.07.22.20116707/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2020/07/26/2020.07.22.20116707/F2)

Figure 2: Recurrent CHIKF episodes.
The distribution of febrile episodes for each of 19 children who had at least two CHIKV RT-PCR positive febrile episodes is shown. The timing of the CHIKV RT-PCR positive results is shown, with the index episode shown as the first data point in the follow up timeline. Children with prefix ‘NGE’ were resident in Ngerenya and those with ‘PIN’ from Pingilikani, respectively. Further details of these children can be found in Supplementary Material Table S2.

### Viral phylogeny

We attempted to sequence all index and subsequent RT-PCR positive samples from the 19 children with recurrent CHIKF episodes (Figure 2). Complete and partial genome coding sequences were generated from 10 samples, including an additional index sample with a low Ct value (Ct=19) from a child who did not experience recurrent CHIKF episodes during the study period. All CHIKV genomes generated in this study mapped to the ECSA genotype, albeit distinct from the clade containing the IOL strains that emerged during the 2004 epidemic and the clade including genomes from CHIKF patients sampled in Mandera, northern Kenya, during the 2016 epidemic (Figure 3). The E1 A226V mutation, associated with adaptation and increased transmissibility by *Aedes albopictus* [15], was absent from all sequences sampled in Kenya [23]. Time-resolved phylogenies showed strong temporal clustering of CHIKV sequences from coastal Kenya that may suggest antigenic changes over time driven by acquisition of herd immunity (Figure 3). Index-recurrent episode case genome sequence pairs were available for four children with respective time intervals of 2.9, 3.5, 22.3 and 28.3 months between the index and recurrent episodes (Figure 3). The degree of relatedness of the sequence pairs was similarly time-dependent, with the sequence pair with the longest time interval (28.3 months) being most divergent (Figure 3).

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/26/2020.07.22.20116707/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2020/07/26/2020.07.22.20116707/F3)

Figure 3: Phylogenetic analysis of CHIKV genomes from Kenya.
BEAST MCC tree of 123 ECSA genomes collected across the world, including 24 from Kenya is shown in the left panel. On the right panel is an extraction highlighting the CHIKV sequences generated in this study. The tips are coloured according to geographic location or child identity (as in Figure 2, for those with paired index-recurrent episode sequence pairs), with posterior probability support shown for select nodes with probability >0.8. The same and more detailed MCC tree can be found in the Supplementary Material Figure S1.

## DISCUSSION

In summary, we find that CHIKF is endemic in coastal Kenya with the case burden being particularly high in young children aged <1 year. Our population-based cohort approach in primary healthcare facilities could detect an increase in CHIKF cases during the 2016 epidemic, supporting the utility of such a surveillance framework in the early detection of CHIKF epidemics. During the epidemic year an increase in CHIKF incidence was only noted in Ngerenya while that in Pingilikani (located approximately 40km south of Ngerenya) had been declining [25]; in our dataset, clinical malaria accounted for 57% of all fevers in Pingilikani but only 5% in Ngerenya, while URTI accounted for 70% of all fevers in Ngerenya but only 10% in Pingilikani. We observed significant CHIKF incidence in both locations, and an inverse relationship between CHIKF incidence and age was evident at both locations. This suggests that immunity is acquired in older children, and that CHIKF transmission is significant in both locations despite the varying malaria ecology.

Consistent with previous reports of CHIKF in children in East Africa [19, 30], none of the patients in this study had a clinical diagnosis of CHIKF or reported symptoms of joint pathology. Nineteen children were CHIKV RT-PCR positive during subsequent febrile episodes, with viral genome sequences from these recurrent episodes not showing high degrees of relatedness to the index case. Chronic relapsing and remitting CHIKF symptoms have been observed in adult patients during previous epidemics, as an inflammatory post-infectious arthropathy [31-33] with viral RNA detected in the joint tissues in one study [34]. However, we are aware of no reports of recurrent discrete CHIKF episodes associated with reinfection. This may simply be due to the rarity of longitudinal surveillance studies for CHIKF, particularly in children. In our study, the relatedness of virus genomes from first and second infections in the same child was similar to the relatedness of virus genomes from different children, suggesting that these were reinfections rather than relapses.

Neutralising antibodies that develop following CHIKV infection have been shown to correlate with decreased risk of clinical illness in a longitudinal study in the Phillipines [35, 36], though CHIKF cases in that study were all due to the Asian CHIKV genotype. The possibility that the ECSA genotype viruses, which tend to be more diverse [12] and associated with a higher clinical disease burden than Asian genotypes [37], results in a variable immune response in children in East Africa needs to be ruled out. The inverse relationship between age and incidence of CHIKF suggests that there is an acquisition of protective immunity, although the recurrent infections suggest that this immunity is not rapidly acquired in some children.

Only finger-prick whole blood samples for RT-PCR analysis were available for this study, thus precluding any immunological work. Future longitudinal studies with serial blood sampling for immunological assays will help determine the natural course of infection, including duration of viraemia, and identify host and viral factors that may underlie the occurrence of CHIKF relapses in coastal Kenya.

Consistent with other reports of CHIKV strains circulating in East Africa [12, 23], the viral genomes generated in this study all mapped to the ECSA genotype. Further, in contrast to the 2016 CHIKF epidemic viruses from Mandera (northern Kenya) that formed a single well-supported clade, sequences from coastal Kenya were highly diverse and fell within multiple well-supported clades that whose divergence increased with time. Together, these data suggest ongoing endemic CHIKV transmission at the Kenyan coast possibly sustained through frequent emergence of diverse CHIKV strains from local sylvatic transmission cycles [6].

Our data and that from others in the region [30] further reinforce the need for inclusion of CHIKF high on the list of differential diagnoses to consider when faced with a febrile child at primary healthcare facilities in coastal Kenya. Access to rapid diagnostics would be of value in antibiotic stewardship, and control of transmission in the community would substantially reduce the burden on healthcare facilities. Further studies on the associated clinical outcomes, including the incidence of severe (‘atypical’) manifestations of CHIKF, should help determine the wider public health significance of CHIKF in coastal Kenya.

## Data Availability

The replication data and analysis scripts for this manuscript shall be made available at the Harvard Dataverse: (https://dataverse.harvard.edu/dataverse/kwtrp). Some of the clinical dataset contains potentially identifying information on participants and is stored under restricted access. Requests for access to the restricted dataset should be made to the Data Governance Committee (dgc@kemri-wellcome.org) of the KEMRI-Wellcome Trust Research Programme.

## FUNDING

This work was commissioned by the National Institute for Health Research (NIHR) Global Health Research programme [16/136/33] using UK aid from the UK Government. The views expressed in this publication are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care. GMW is supported by an Oak foundation fellowship and a Wellcome Trust grant [grant number 203077\_Z\_16_Z].

## CONFLICT OF INTEREST

All authors declare no conflicts of interest.

## ACKNOWLEDGEMENTS

This manuscript was submitted for publication with permission from the Director of the Kenya Medical Research Institute.

*   Received July 22, 2020.
*   Revision received July 22, 2020.
*   Accepted July 26, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## REFERENCES

1.  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.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0035-9203(55)90080-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14373834&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

2.  2.Lumsden WH. An epidemic of virus disease in Southern Province, Tanganyika Territory, in 1952-53. II. General description and epidemiology. Trans R Soc Trop Med Hyg 1955; 49(1): 33–57.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0035-9203(55)90081-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14373835&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

3.  3.Tritsch SR, Encinales L, Pacheco N, et al. Chronic joint pain 3 years after chikungunya virus infection largely characterized by relapsing-remitting symptoms. J Rheumatol 2019.
    
    

4.  4.Ross RW. The Newala epidemic. III. The virus: isolation, pathogenic properties and relationship to the epidemic. J Hyg (Lond) 1956; 54(2): 177–91.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0022172400044442&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=13346078&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

5.  5.Althouse BM, Guerbois M, Cummings DAT, et al. Role of monkeys in the sylvatic cycle of chikungunya virus in Senegal. Nat Commun 2018; 9(1): 1046.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-03332-7&link_type=DOI) 

6.  6.Eastwood G, Sang RC, Guerbois M, Taracha ELN, Weaver SC. Enzootic Circulation of Chikungunya Virus in East Africa: Serological Evidence in Non-human Kenyan Primates. Am J Trop Med Hyg 2017; 97(5): 1399–404.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4269/ajtmh.17-0126&link_type=DOI) 

7.  7.Weaver SC, Lecuit M. Chikungunya virus and the global spread of a mosquito-borne disease. N Engl J Med 2015; 372(13): 1231–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMra1406035&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25806915&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

8.  8.Wahid B, Ali A, Rafique S, Idrees M. Global expansion of chikungunya virus: mapping the 64-year history. Int J Infect Dis 2017; 58: 69–76.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2017.03.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

9.  9.Voss JE, Vaney MC, Duquerroy S, et al. Glycoprotein organization of Chikungunya virus particles revealed by X-ray crystallography. Nature 2010; 468(7324): 709–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature09555&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21124458&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000284836700044&link_type=ISI) 

10. 10.Rezza G, Weaver SC. Chikungunya as a paradigm for emerging viral diseases: Evaluating disease impact and hurdles to vaccine development. PLoS Negl Trop Dis 2019; 13(1): e0006919.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0006919&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30653504&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

11. 11.Gouglas D,  Thanh Le T, Henderson K, et al. Estimating the cost of vaccine development against epidemic infectious diseases: a cost minimisation study. Lancet Glob Health 2018; 6(12): e1386–e96.
    
    

12. 12.Volk SM, Chen R, Tsetsarkin KA, et al. Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates. J Virol 2010; 84(13): 6497–504.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjEwOiI4NC8xMy82NDk3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMjYvMjAyMC4wNy4yMi4yMDExNjcwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

13. 13.Chretien JP, Anyamba A, Bedno SA, et al. Drought-associated chikungunya emergence along coastal East Africa. Am J Trop Med Hyg 2007; 76(3): 405–7.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI3Ni8zLzQwNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzI2LzIwMjAuMDcuMjIuMjAxMTY3MDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

14. 14.Sergon K, Njuguna C, Kalani R, et al. Seroprevalence of Chikungunya virus (CHIKV) infection on Lamu Island, Kenya, October 2004. Am J Trop Med Hyg 2008; 78(2): 333–7.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI3OC8yLzMzMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzI2LzIwMjAuMDcuMjIuMjAxMTY3MDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

15. 15.Tsetsarkin KA, Vanlandingham DL, McGee CE, Higgs S. A single mutation in chikungunya virus affects vector specificity and epidemic potential. PLoS Pathog 2007; 3(12): e201.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.0030201&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18069894&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

16. 16.Weaver SC, Forrester NL. Chikungunya: Evolutionary history and recent epidemic spread. Antiviral Res 2015; 120: 32–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.antiviral.2015.04.016&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25979669&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

17. 17.Nsoesie EO, Kraemer MU, Golding N, et al. Global distribution and environmental suitability for chikungunya virus, 1952 to 2015. Euro Surveill 2016; 21(20).
    
    

18. 18.Kajeguka DC, Kaaya RD, Mwakalinga S, et al. Prevalence of dengue and chikungunya virus infections in north-eastern Tanzania: a cross sectional study among participants presenting with malaria-like symptoms. BMC Infect Dis 2016; 16: 183.
    
    

19. 19.Chipwaza B, Mugasa JP, Selemani M, et al. Dengue and Chikungunya fever among viral diseases in outpatient febrile children in Kilosa district hospital, Tanzania. PLoS Negl Trop Dis 2014; 8(11): e3335.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0003335&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25412076&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

20. 20.Hertz JT, Munishi OM, Ooi EE, et al. Chikungunya and dengue fever among hospitalized febrile patients in northern Tanzania. Am J Trop Med Hyg 2012; 86(1): 171–7.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI4Ni8xLzE3MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzI2LzIwMjAuMDcuMjIuMjAxMTY3MDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

21. 21.Weller N, Clowes P, Dobler G, et al. Seroprevalence of alphavirus antibodies in a cross-sectional study in southwestern Tanzania suggests endemic circulation of chikungunya. PLoS Negl Trop Dis 2014; 8(7): e2979.
    
    

22. 22.LaBeaud AD, Banda T, Brichard J, et al. High rates of o’nyong nyong and Chikungunya virus transmission in coastal Kenya. PLoS Negl Trop Dis 2015; 9(2): e0003436.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0003436&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25658762&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

23. 23.Maljkovic Berry I, Eyase F, Pollett S, et al. Global Outbreaks and Origins of a Chikungunya Virus Variant Carrying Mutations Which May Increase Fitness for Aedes aegypti: Revelations from the 2016 Mandera, Kenya Outbreak. Am J Trop Med Hyg 2019; 100(5): 1249–57.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4269/ajtmh.18-0980&link_type=DOI) 

24. 24.Scott JA, Bauni E, Moisi JC, et al. Profile: The Kilifi Health and Demographic Surveillance System (KHDSS). Int J Epidemiol 2012; 41(3): 650–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys062&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22544844&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306417300014&link_type=ISI) 

25. 25.Mogeni P, Williams TN, Fegan G, et al. Age, Spatial, and Temporal Variations in Hospital Admissions with Malaria in Kilifi County, Kenya: A 25-Year Longitudinal Observational Study. PLoS Med 2016; 13(6): e1002047.
    
    

26. 26.Lanciotti RS, Kosoy OL, Laven JJ, et al. Chikungunya virus in US travelers returning from India, 2006. Emerg Infect Dis 2007; 13(5): 764–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1305.070015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17553261&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246212400018&link_type=ISI) 

27. 27.Pabbaraju K, Wong S, Gill K, Fonseca K, Tipples GA, Tellier R. Simultaneous detection of Zika, Chikungunya and Dengue viruses by a multiplex real-time RT-PCR assay. J Clin Virol 2016; 83: 66–71.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jcv.2016.09.001&link_type=DOI) 

28. 28.Kozlov AM, Darriba D, Flouri T, Morel B, Stamatakis A. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 2019; 35(21): 4453–5.
    
    

29. 29.Bejon P, Williams TN, Nyundo C, et al. A micro-epidemiological analysis of febrile malaria in Coastal Kenya showing hotspots within hotspots. Elife 2014; 3: e02130.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.02130&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24843017&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

30. 30.Waggoner J, Brichard J, Mutuku F, et al. Malaria and Chikungunya Detected Using Molecular Diagnostics Among Febrile Kenyan Children. Open Forum Infect Dis 2017; 4(3): ofx110.
    
    

31. 31.Couturier E, Guillemin F, Mura M, et al. Impaired quality of life after chikungunya virus infection: a 2-year follow-up study. Rheumatology (Oxford) 2012; 51(7): 1315–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/rheumatology/kes015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22427407&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305631500024&link_type=ISI) 

32. 32.Staikowsky F, Le Roux K, Schuffenecker I, et al. Retrospective survey of Chikungunya disease in Reunion Island hospital staff. Epidemiol Infect 2008; 136(2): 196–206.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/s0950268807008424&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17433130&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259337000006&link_type=ISI) 

33. 33.Chang AY, Encinales L, Porras A, et al. Frequency of Chronic Joint Pain Following Chikungunya Virus Infection: A Colombian Cohort Study. Arthritis Rheumatol 2018; 70(4): 578–84.
    
    

34. 34.Hoarau JJ, Jaffar Bandjee MC, Krejbich Trotot P, et al. Persistent chronic inflammation and infection by Chikungunya arthritogenic alphavirus in spite of a robust host immune response. J Immunol 2010; 184(10): 5914–27.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE4NC8xMC81OTE0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMjYvMjAyMC4wNy4yMi4yMDExNjcwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

35. 35.Yoon IK, Alera MT, Lago CB, et al. High rate of subclinical chikungunya virus infection and association of neutralizing antibody with protection in a prospective cohort in the Philippines. PLoS Negl Trop Dis 2015; 9(5): e0003764.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0003764&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25951202&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F26%2F2020.07.22.20116707.atom) 

36. 36.Yoon IK, Srikiatkhachorn A, Alera MT, Fernandez S, Cummings DAT, Salje H. Pre-existing chikungunya virus neutralizing antibodies correlate with risk of symptomatic infection and subclinical seroconversion in a Philippine cohort. Int J Infect Dis 2020.
    
    

37. 37.Bustos Carrillo F, Collado D, Sanchez N, et al. Epidemiological Evidence for Lineage-Specific Differences in the Risk of Inapparent Chikungunya Virus Infection. J Virol 2019; 93(4).