Genome-Scale Phylogenetic Analyses of Chikungunya Virus Reveal Independent Emergences of Recent Epidemics and Various Evolutionary Rates

Chikungunya virus (CHIKV), a mosquito-borne alphavirus, has traditionally circulated in Africa and Asia, causing human febrile illness accompanied by severe, chronic joint pain. In Africa, epidemic emergence of CHIKV involves the transition from an enzootic, sylvatic cycle involving arboreal mosquito vectors and nonhuman primates, into an urban cycle where peridomestic mosquitoes transmit among humans. In Asia, however, CHIKV appears to circulate only in the endemic, urban cycle. Recently, CHIKV emerged into the Indian Ocean and the Indian subcontinent to cause major epidemics. To examine patterns of CHIKV evolution and the origins of these outbreaks, as well as to examine whether evolutionary rates that vary between enzootic and epidemic transmission, we sequenced the genomes of 40 CHIKV strains and performed a phylogenetic analysis representing the most comprehensive study of its kind to date. We inferred that extant CHIKV strains evolved from an ancestor that existed within the last 500 years and that some geographic overlap exists between two main enzootic lineages previously thought to be geographically separated within Africa. We estimated that CHIKV was introduced from Africa into Asia 70 to 90 years ago. The recent Indian Ocean and Indian subcontinent epidemics appear to have emerged independently from the mainland of East Africa. This finding underscores the importance of surveillance to rapidly detect and control African outbreaks before exportation can occur. Significantly higher rates of nucleotide substitution appear to occur during urban than during enzootic transmission. These results suggest fundamental differences in transmission modes and/or dynamics in these two transmission cycles.Chikungunya virus (CHIKV; Togaviridae: Alphavirus) is an arbovirus (arthropod-borne virus) vectored by Aedes mosquitoes to humans in tropical and subtropical regions of Africa and Asia (Fig. ​(Fig.1;1; reviewed in references 26 and 46). CHIKV has a single-stranded, positive-sense RNA genome of ∼12 kb and causes chikungunya fever (CHIK), a febrile illness associated with severe arthralgia and rash (2, 15, 31, 35); the name is derived from a Bantu language word describing the severe arthritic signs (32), which can persist for years. Thus, CHIK has enormous economic costs in addition to its public health impact (9). Because the signs and symptoms of CHIK overlap with those of dengue and because CHIKV is transmitted sympatrically in urban areas by the same mosquito vectors, it is grossly underreported in the absence of laboratory diagnostics (10, 37).Open in a separate windowFIG. 1.Distribution of the CHIKV strains used in this study. The map, based on a world map template from http://www.presentationmagazine.com, was edited with permission.CHIKV was first isolated during a 1953 outbreak in present-day Tanzania by Ross (48, 49). Since then, outbreaks have been documented in Africa and Asia, including the Indian subcontinent (Fig. ​(Fig.1)1) (1, 4). In 2005, CHIKV emerged from East Africa to cause an explosive urban epidemic in popular tourist island destinations in the Indian Ocean (Fig. ​(Fig.1;1; reviewed in reference 31). In late 2005, CHIKV spread into the Indian subcontinent, where millions of people have been affected (5). However, the geographic source of spread into India, from the mainland of Africa or from the Indian Ocean Islands, has not been delineated. India had seen large epidemics of CHIK in the past (reviewed in reference 30), but CHIKV apparently disappeared during the 1970s (5). Since 2006, CHIKV has been imported into Europe and the western hemisphere (including the United States) via many viremic travelers, and an epidemic was initiated in Italy by a traveler from India (4, 11, 47). The dramatic spread since 1980 of dengue viruses (DENV) throughout tropical America, via the same vectors, portends the severity of the public health problem if CHIKV becomes established in the western hemisphere.The first phylogenetic analysis of CHIKV (45) identified three geographically associated genotypes: the West African (WAf), East/Central/South African (ECSA), and Asian genotypes. More recent analyses indicate that the recent Indian Ocean and Indian strains form a monophyletic group within the ECSA lineage (5, 12, 14, 27, 40, 51, 52). However, most CHIKV phylogenetic studies (1, 14, 28, 29, 38, 40, 41, 47, 52) have utilized only partial sequences from the envelope glycoprotein E1 gene, preventing a robust assessment of some of the relationships among strains and of their evolutionary dynamics.The CHIKV strains represented in different geographic lineages apparently circulate in different ecological cycles. In Asia, CHIKV appears to circulate primarily in an urban transmission cycle involving the peridomestic mosquitoes Aedes aegypti and A. albopictus, as well as humans (25, 45). Asian epidemics typically infect thousands-to-millions of people over the course of several years (46). In contrast, African CHIKV circulates primarily in a sylvatic/enzootic cycle, transmitted by arboreal primatophilic Aedes mosquitoes (e.g., A. furcifer and A. africanus) and probably relies on nonhuman primates as reservoir hosts (reviewed in reference 16). Epidemics in rural Africa usually occur on a much smaller scale than in Asia, likely a result of the lower human population densities, and possibly more stable herd immunity. Although the assignments of “urban” and “sylvatic/enzootic” are based on the most common mode of transmission, CHIKV strains of African origin are capable of urban transmission by A. aegypti and A. albopictus, as evidenced by outbreaks in the Democratic Republic of the Congo (41), Nigeria (36), Kenya (27), and Gabon (42). The ecological differences between the sylvatic/enzootic (henceforth called enzootic) and urban/endemic/epidemic transmission cycles (henceforth called epidemic) such as seasonality of vector larval habitats, vertebrate host abundance and herd immunity, and vector host preferences, prompted us to hypothesize that the evolutionary dynamics of CHIKV may differ between the two transmission cycles. To test this hypothesis, to provide more robust estimates of the evolutionary relationships among the CHIKV strains including the sources of the recent epidemics, and to elucidate the temporal and spatial history of CHIKV evolution, we performed an extensive, genome-scale phylogenetic analysis, utilizing complete open reading frame (ORF) sequences of a large collection of 80 isolates with broad temporal, spatial, and host coverage.