1,652
Views
33
CrossRef citations to date
0
Altmetric
Original Articles

Re-visiting the evolution, dispersal and epidemiology of Zika virus in Asia

, ORCID Icon, , ORCID Icon, , , , , ORCID Icon, , & show all
Pages 1-8 | Received 10 Jan 2018, Accepted 03 Apr 2018, Published online: 09 May 2018

Abstract

Based on serological evidence and viral isolation, Zika virus (ZIKV) has circulated for many years relatively benignly in a sylvatic cycle in Africa and an urban cycle in South East Asia (SEA). With the recent availability of limited but novel Indian ZIKV sequences to add to the plethora of SEA sequences, we traced the phylogenetic history and spatio-temporal dispersal pattern of ZIKV in Asia prior to its explosive emergence in the Pacific region and the Americas. These analyses demonstrated that the introduction and dispersal of ZIKV on the Pacific islands were preceded by an extended period of relatively silent transmission in SEA, enabling the virus to expand geographically and evolve adaptively before its unanticipated introduction to immunologically naive populations on the Pacific islands and in the Americas. Our findings reveal new features of the evolution and dispersal of this intriguing virus and may benefit future disease control strategies.

Introduction

Zika virus (ZIKV) was first isolated in 1947 from rhesus monkey serum collected in the Zika Forest of Uganda (Africa)Citation1. Serosurveys conducted from the 1950s onwards suggested that ZIKV circulated in several countries of Africa and Asia. Less than 20 human infections were confirmed within the first 60 years following the discovery of the virus, probably because it causes mild disease in humans and is often misdiagnosed given the existence of several clinically similar diseases, for example, dengue or chikungunya feverCitation2, Citation3. ZIKV emerged in Africa (African lineage) and subsequently spread through Asia (Asian Lineage). In 2007, an Asian lineage strain of ZIKV unpredictably caused a unique explosive epidemic on the Pacific island of Yap in MicronesiaCitation4, involving 70% of the population. However, within months, the epidemic subsided, and there has only been one further report of epidemic ZIKV in Micronesia in the Kosrae State in 2016Citation5. However, in October 2013, a more recent ZIKV strain emerged in French Polynesia (FP)Citation6, probably following the accumulation of functionally important amino acid substitutionsCitation7, Citation8, which were predicted to enhance virus transmission by the primary vector species Aedes aegyptiCitation7. Subsequently, non-vector borne transmission (mother to child, sexual and blood transfusion transmission) was confirmed with the emergence of ZIKV in the AmericasCitation2, Citation9. Following its establishment on the FP islands, ZIKV radiated eastwards and westwardsCitation10 and subsequently emerged in Brazil before dispersing and causing epidemics in the AmericasCitation2, Citation11.

Many new sequences are now available from ZIKV-infected humans in Asia, including BangladeshCitation12, IndiaCitation13, IndonesiaCitation14, the PhilippinesCitation15, SingaporeCitation16, ThailandCitation17, and VietnamCitation18. These sequences may impact our interpretation of the early evolutionary history of ZIKV in Asia. We therefore re-visited the evolutionary events of ZIKV from its first emergence in Asia. Although additional sequence data will expand our knowledge, it is already possible with the new data from India and South-East Asia (SEA) to interpret the dispersal patterns of ZIKV in terms of its geographic source, mode of spread, and time of circulation.

Using updated phylogeographic and phylogenetic analyses, we illustrate and discuss how the epidemic pattern in FP and the Americas contrasts strongly with the relatively long-term and apparently benign circulation of ZIKV in SEA.

Results

Bayesian phylogenetic analysis of ZIKV sequences from Africa and Asia, including partial data from a recently available Indian ZIKV sequenceCitation13, indicate that the Indian virus shares a more recent ancestry with all of the other Asian sequences than with the Malaysian sequence, which is positioned as a sister-group to all of the other Asian sequences (Fig. , branch leading to node B, Posterior probability 0.71; Suppl. fig. S1). We also performed separate Bayesian analyses of both the partial capsid and NS2b/NS3-genes (Suppl. fig. S2 and S3, respectively). Our results demonstrated that Indian ZIKV shared a more recent ancestry with the SEA viruses compared to the Malaysian virus. However, a tree based on a partial envelope gene sequence placed the Indian virus in the most basal position of the Asian ZIKV clade (Suppl. fig. S4). Whether this truly represents the phylogenetic position of Indian ZIKV will require a complete genomic analysis of the isolate and more sequence data from geographically distant regions of India and neighbouring countries. Subsequently, and separated by a long branch, ZIKV then diverged to form a clade that harboured the majority of Asian ZIKV diversity (Fig , node C). This was followed by a series of diversification events (Fig. , nodes C–L) that predominantly occurred in SEA and were characterized by frequent introductions and localized outbreaks that probably reflected the increasing distribution and spread of ZIKV in the primary vector species, i.e., Ae aegypti, and susceptible humans. Notably, both Thailand and Vietnam experienced introductions and outbreaks of multiple ZIKV strains throughout the course of its evolution in SEA and have served as sources of ZIKV in other regions (Fig. , Suppl. fig. S1 and S5).

Fig. 1 Illustrative phylogenetic tree of ZIKV evolution and dispersal in the Asian region.

Coloured circles on the terminal nodes represent the country of infection on the map. White lines indicate the main routes of spread. Letters within the gray circles (A–M4) refer to diversification events throughout the period of ZIKV evolution in Asia and the Pacific islands, which correspond to the dates of divergence (in bold) and highest posterior density for each tMRCA. A comprehensive tree is provided in the Supplementary Materials (Suppl. fig. S5)

Fig. 1 Illustrative phylogenetic tree of ZIKV evolution and dispersal in the Asian region.Coloured circles on the terminal nodes represent the country of infection on the map. White lines indicate the main routes of spread. Letters within the gray circles (A–M4) refer to diversification events throughout the period of ZIKV evolution in Asia and the Pacific islands, which correspond to the dates of divergence (in bold) and highest posterior density for each tMRCA. A comprehensive tree is provided in the Supplementary Materials (Suppl. fig. S5)

An inclusive phylogenetic and temporal analysis, focusing on isolates positioned after the dispersal of ZIKV from Africa but prior to its emergence on the Pacific islands and the Americas, was prepared. The evolutionary rate of the Asian ZIKV lineage under a strict clock was estimated to be 7.26 × 10−4 substitutions per site per year (95% HPD: 6.28–8.19 × 10−4). Furthermore, the analysis showed that the most recent common ancestor (MRCA) of ZIKV in Asia (Fig. , node A) was present 5–12 years earlier than previously estimatedCitation7, Citation19, ~1954–04 (95% highest posterior density: 1948–08 to 1959–05). Then, ~2 years later, circa 1956–10 (95% HPD: 1949–09 to 1964–05), the MRCA of the Indian lineage and all subsequent Asian lineage ZIKV emerged (Fig. , node B), surprisingly implying the introduction of Indian ZIKV from an Asian rather than an African source. Subsequently, following a silent period of ~40 years, ZIKV appeared in SEA in early 1997 (95% HPD: 1994–09 to 1999–03; Fig. , node C). It should also be noted that the MRCA of the Micronesian and Philippine sequences was estimated to have been present in Asia ca. 1997 (95% HPD: 1995–07 to 1999–12; Suppl. fig. S5). Thus, the ancestor of these two emergent viruses was possibly present in SEA and evolved relatively silently for at least 10 years prior to the Micronesian epidemic in 2007. Following the emergence of ZIKV in SEA (node C), there was a period of circulation of ~15 years (Fig. , nodes C–L) prior to its introduction and emergence on the FP islands (Fig. , nodes M1–M4) in mid-2012 (95% HPD: 2011–12 to 2013–01). The introduction of ZIKV from FP into America followed about 1 year later, ~2013–03 (95% HPD: 2012–11 to 2013–06) (Suppl. table 1; Suppl. fig. S5). When the Indian sequence data were excluded, the results of the temporal estimates were confirmed by using the HPD values that overlapped in the two analyses, which were supported by time-regression analyses including all data (Suppl. fig. S5S7).

Discussion

Following its first isolation in April 1947 from a sentinel rhesus monkey in the canopy of the African Zika Forest of UgandaCitation1, ZIKV was isolated again in 1948 in the same forest but this time from Ae africanus, an indigenous African sylvatic mosquitoCitation1. These observations suggest that in Africa, ZIKV circulates between sylvatic Aedes species mosquitoes and monkeys that inhabit the forest canopy, which does not exclude other possible transmission cycles. Subsequent epidemiological and serological studies have revealed that ZIKV was probably also prevalent and widespread in African towns and villages during the 1950sCitation2, Citation20, Citation21. Since it is now known that both dengue viruses (DENVs) and ZIKV adapted from sylvatic to urban transmission cycles, likely due to the increasing Ae aegypti population density in urban districts of SEA, our observations provide a rational explanation for how it was possible for ZIKV to be isolated in 1966 from domestic Ae aegypti in populous areas of the Malaysian peninsulaCitation22. Importantly, attempts in the 1950s to isolate ZIKV from a wide range of mosquitoes in sylvatic environments in MalaysiaCitation22 and other Asian countries did not result in the identification of additional strains of ZIKV, implying that the Malaysian virus was not significantly established in the sylvatic Asian environment, which also appears to be the case for ZIKV in the Pacific and the Americas; that is, ZIKV is not significantly sylvatic except in Africa. Nevertheless, following the detection of ZIKV in Malaysia, subsequent serosurveys conducted between 1952 and 1989, in which 189 of 1011 tests were reported as ZIKV positiveCitation2, Citation23Citation25, support the concept of widespread distribution of ZIKV in SEACitation2, Citation19, Citation26Citation29. Importantly, these early studies were performed using the most specific serological assay, viz., neutralization against high virus challenge doses (>50 LD50 to >1000 LD50) in mice, to reveal a high incidence of neutralizing antibody responses in residents of Malaya (19% and 75%)Citation23, Citation25, Borneo (18%)Citation23 and India (16.8%)Citation24.

Why did ZIKV disperse so successfully throughout Asia for so many years without appearing to have caused major epidemicsCitation30? Clearly, the very low frequency of clinically apparent ZIKV infections and the difficulty in clinically discriminating between infections caused by DENV, chikungunya virus (CHIKV), and ZIKV were major factors in the failure to recognize the presence of ZIKV infections. Second, the inapparent spread of endemic ZIKV throughout SEA resulted in the emergence of low-level herd immunity to the virus. Third, antigenically related flaviviruses, including DENV, West Nile virus (WNV), Japanese encephalitis virus, Wesselsbron virus, and Stratford virus, circulating throughout SEA contributed a background of cross-protective reactivity to ZIKV that might have dampened ZIKV transmission. By contrast, with the exception of several islands in the Western Pacific region that have experienced outbreaks of JEV, DENV was the only human pathogenic flavivirus known to have caused outbreaks in islands further east in the Pacific prior to the arrival of epidemic ZIKV. Moreover, it is a widely held view that immunity to DENV alone does not protect humans against infection with ZIKVCitation2. Thus, the evolving ZIKV was probably introduced to populations of immunologically naive inhabitants on most of the Pacific islands, leading to the explosive epidemics that have now been reported widelyCitation7. Several human flaviviruses, including DENVs, are known to circulate in the Americas. Thus, the situation in Brazil and neighbouring South American countries should differ significantly from that described in Oceania because several human pathogenic ZIKV-related flaviviruses circulate in Brazil, including DENVs, St. Louis encephalitis virus (SLEV), yellow fever virus (YFV), Ilheus virus (ILHV), Rocio virus (ROCV), and others. Why, therefore, was ZIKV able to cause explosive epidemics in the Americas in the presence of this potential immune background? The most likely explanation is that, with the exception of YFV, which circulates primarily in the sylvatic environment in the Americas, and DENV, which is not believed to stimulate protective immunity to ZIKVCitation2, these other ZIKV-related flaviviruses are not commonly transmitted by domestic Ae aegypti and therefore do not circulate primarily in the urban environment favored by ZIKV.

Moreover, a serosurvey of Indian human sera using the mouse protection test revealed a high incidence of protective immunity against DENV-1 (40.4%: 72/196), DENV-2 (20.2%: 38/188), WNV (35.0%: 123/351), and ZIKV (16.8%: 33/196), while none of 588 sera neutralized YFV and only a few neutralized the American flaviviruses, ILHV (1.4%: 3/206) or SLEV (2.8%: 6/211)Citation24. These results suggest that the presence of an antibody-mediated cross-protection between ZIKV and either YFV or the American flaviviruses is relatively unlikely. In addition, while ROCV caused local epidemics in coastal areas of São Paulo, Brazil from 1973 to 1980, it has subsequently disappeared from this region of BrazilCitation30. Furthermore, WNV, which has only recently been detected in Brazil, does not yet appear to have become established as a significant pathogen in Latin AmericaCitation31. It should also be emphasized that even following the epidemics in the Pacific and subsequent spread of ZIKV to the Americas, Asia has experienced relatively small ZIKV outbreaks compared to the explosive and clinically apparent epidemics in the Pacific and the AmericasCitation32.

Thus, ZIKV epidemics in FP and the Americas have occurred in the presence of a high background immunity to DENV. However, as noted above, it is a widely held view that immunity to DENV does not necessarily protect humans against infection with ZIKVCitation2. How and when did African ZIKV gain access to SEA? Like DENV and CHIKV, the major transmission vector of ZIKV is Ae aegypti, the domestic descendant of the sylvatic African Ae africanusCitation2, Citation33 and the primary conduit for the introduction of DENV and CHIKV to Asia via ships trading out of East Africa during the 18th–20th centuriesCitation34Citation36. Thus, by analogy, ZIKV is also likely to have been introduced to Asia, from Africa, in infected humans and/or Ae aegypti, to which ZIKV had adapted, via the same shipping routes. This argument is also supported by contemporary evidence of frequently reported cases of ZIKV being introduced to non-tropical countries by infected individuals returning from areas in the tropics where epidemic ZIKV was known to be circulatingCitation37. Based on an earlier studyCitation7 with a similar estimated evolutionary rate (Median: 6.30 × 10−4; 95% HPD: 5.48 × 10−4 to 7.10 × 10−4) to that of the present study, the African and Asian ZIKV lineages were estimated to have diverged in ~1834, which is consistent with reports of other Ae aegypti- transmitted viruses, namely, DENV and CHIKV, being introduced via ships trading between East Africa through India to SEACitation34, Citation38.

Recently, ZIKV was isolated from a patient in northern India. A neighbour-joining analysis of the partial envelope gene of this virus implied that the Indian sequence was basal in the Asian ZIKV cladeCitation13. To investigate this further, we concatenated sequences of the partial capsid, envelope and NS2b/NS3 genes of the Indian ZIKV isolate and analyzed the data in a Bayesian phylogenetic framework. Our results suggested that the Indian sequence was not basal to the Malaysian virus isolated in 1966 (Suppl. fig. S1). By contrast, our analysis suggests that ZIKV could have been introduced to India from SEA, rather than vice versa. Although the currently available sequence data do not robustly resolve the Indian versus Malaysian ancestry, the concept of the introduction of ZIKV from SEA to India seems plausible given that ZIKV also appears to have been introduced to Bangladesh from SEA (Fig. ; Suppl. fig. S1 and S5). In this context, the highest incidence of human sera that neutralized ZIKV was collected from towns in the Ahmedabad district (north-west India) in the state from which the NIV1720741 strain used in the present analysis was isolated in 2016Citation24. This state may, therefore, represent a pocket of endemicity containing divergent ZIKV strains in that region.

Our phylogenetic and temporal analyses suggest that ZIKV was first introduced to SEA during the 1950s, and then dispersed both eastwards and westwards ~5–12 years earlier than was previously predictedCitation7, Citation19. Our estimate is also compatible with serological studies conducted in Malaysia, India, Indonesia, the Philippines, Thailand and Vietnam, indicating that ZIKV was present in these countries between 1952 and 1954Citation2 (and references therein). Based on the current phylogenies (Fig. ; Suppl. fig. S1 and S5), ZIKV was potentially first introduced from Africa to the Malaysian–Indonesian region and subsequently to India. However, although many new sequences have been published during the past few years, it is still not possible to determine precisely whether any particular country in SEA has played a major role compared with other countries in the evolution and dispersal of ZIKV in SEA.

Between 1954 and 1976, no recorded human clinical ZIKV cases in Asia were confirmed using serological tests or virus isolation procedures. Hence, it is likely that ZIKV remained undetected because it was primarily a sub-clinical disease and diagnostic tests for ZIKV were not routinely performed. Then, in 1977 and 1978 three and four febrile patients on Java, Indonesia, were diagnosed with ZIKV infection, respectivelyCitation39. Over 35 years later, ZIKV was isolated from a patient from Sumatra, IndonesiaCitation14. It is possible that ZIKV remained undetected in the Malaysian–Indonesian region for a long time, circulating mainly between Ae aegypti and human hosts, causing sporadic sub-clinical or mildly clinical human infections, while evolving and adapting to its vector species. Thus, since the late 1950s, ZIKV appears to have circulated silently in SEA according to serosurvey data. The analyses suggest that the MRCA of the Malaysian and Indian isolates is approximately 40 years earlier than that of the Micronesian/Philippine and all subsequent isolates (Fig. , node B vs C, Suppl. fig. S5). However, without further sequence data from these neglected areas of Asia, it will continue to be difficult to identify precisely where ZIKV was circulating and, potentially, causing human ZIKV infections. Nevertheless, serological studies suggest that ZIKV may have been present in the Indian–SEA region earlier than the 1950sCitation23, Citation25, Citation40.

It is clear that ZIKV has circulated without recognition for many years in Asia. This is supported by evidence from serology, virus isolationCitation2 and phylogenetic analyses, all of which have shown that the ancestral Asian virus continued to evolve and circulate (Suppl. fig. S1 and S5) prior to and after its unexpected emergence as an epidemic virus on Yap island, Micronesia, in 2007Citation4. In fact, during a period of 14 years, 1997–2011 (Fig. , nodes C–L), this ancestral lineage of Asian ZIKV appears to have dispersed between several other countries in SEA. Notably, Thailand and Vietnam were possibly the source of ZIKV in other regions of Asia. For example, Thailand appears to have been the source of ZIKV in Singapore, which experienced a large outbreak in 2016Citation16, and in Bangladesh (Suppl. fig S5). Thus, it is likely that different countries in SEA have experienced different epidemiological situations. Whereas, the Philippines seems to have had a single strain of ZIKV circulating for several years, Thailand and Vietnam appear to have experienced multiple introductions and “waves” of ZIKV that have swept across these countries (Suppl. fig S1 and S5) followed by re-occurring outbreaksCitation29, Citation41. Consequently, following the isolated Micronesian epidemic in 2007, ZIKV continued to circulate and evolve in the SEA–Indian-regionCitation32, Citation42 before causing a large-scale epidemic in FP in 2013/2014. ZIKV was subsequently introduced into the Americas, apparently via BrazilCitation7, Citation8, Citation11.

It is also believed that the ZIKV strain that emerged on Yap Island, Micronesia, did not disperse widely beyond Micronesia. Although the reasons for its failure to spread have not been identified, there are several possible explanations. First, the Micronesian lineage pre-dates the FP lineage by at least six years, during which four unique and functionally significant amino acid substitutions emerged in the FP lineage, i.e., T777M, V763M, S139N, and M/T2634VCitation7. The first three of these unique substitutions has been shown to have functions relating either to receptor attachment or virus fusion activity with cell membranesCitation43. These substitutions could therefore have an impact on the relative transmission efficiency of the FP lineage compared with earlier Asian lineage sequencesCitation7. Second, in support of this interpretation, the unique S139N substitution, first identified in 2016Citation7, emerged in FP and was detected in virtually every descendant lineage, suggesting a possible link with microcephaly. Subsequently, the unique S139N substitution was associated with increased ZIKV infectivity in humans and enhanced fetal mouse microcephalyCitation44. Third, the relatively small human population size and density on the Micronesian islands and limited connectivity with other Pacific islands might have limited the likelihood of ZIKV dispersal to other distant regions of the Pacific. Fourth, a mutation (A188V) was subsequently shown to potentially contribute to the apparently increased transmission efficiency of ZIKVCitation45. However, more than one amino acid substitution appears to determine the transmission efficiency since the A188V substitution is not consistently associated with it. For example, A188V is the ancestral state for this described substitution in Africa, but the only outbreak reported in Africa was caused by a strain belonging to the Asian lineageCitation46. Nevertheless, the data, in total, are consistent with the hypothesis that ZIKV has circulated, potentially silently in the presence of other related human pathogenic flaviviruses, in sub-Saharan Africa and SEA for many decadesCitation47Citation49.

Our analysis also indicates that the MRCA of the lineage that eventually led to the large outbreak on FP might have existed and potentially been introduced in June 2012 (95% HPD: 2011–12 to 2013–01), i.e., ~1 year before the first detection of ZIKV in FP in October 2013Citation6, Citation9, Citation50, Citation51. Based on the observed short divergence times between nodes M1 and N (Suppl. table 1; Suppl. fig. S5), a potential explanation for the rapid dispersal of ZIKV in the Pacific and the Americas could be interpreted as being due to relatively low levels of herd immunity in the populations. In addition, a similar cryptic period of about 1 year, before its appearance as an epidemic virus, was predicted for the introduction of ZIKV into BrazilCitation11. Whether this reflects the relatively mild nature of ZIKV infections and consequently low viraemic levels in the majority of the infected population is an intriguing subject for mathematical modelers.

In summary, following its emergence in Asia, after being introduced from Africa, ZIKV circulated as a relatively benign virus in Asia for several decades. ZIKV then emerged in Oceania and the Americas and, largely because of its association with central nervous system malformations in neonates, especially microcephaly, it became an international public health concern. This raises at least two important questions, (i) how did ZIKV circulate successfully throughout Asia during its apparent silent period between 1966 and 2006 and (ii) what are the major determinants of the increased epidemicity of ZIKV in the Pacific and the Americas. With regard to the former question, several possible explanations have been proposed for these differences in the epidemicity of ZIKV in Asia compared to that of ZIKV in the Pacific and the Americas: Stochastic evolutionCitation52, adaptive evolutionCitation7, NS1 protein enhancementCitation45, S139 enhancementCitation7, Citation44, and contrasting levels of background immunity due to cross-protective antibodiesCitation2. Each of these factors could contribute to the overall differences in the transmission efficiency. However, the relatively high background cross-reactive antigenicity and/or immunity, which are likely present in Asia, contrasting with the relatively low equivalent antigenicity and/or immunity in the Pacific and the Americas, in combination with the key amino acid mutations in the ZIKV genome, can potentially explain the different epidemiological patterns between Asia on the one hand and the Pacific and Americas on the other.

Materials and methods

Partial and complete genome ZIKV sequences were retrieved from the NCBI GenBank (www.ncbi.nlm.nih.gov/genbank/) focusing on Asian lineage ZIKV isolates. Initially, all available complete genomes (as of October 2017) were retrieved. Sequences collected from the Pacific islands and the Americas were subsampled to reduce the data load. Subsequently, partial phylogenetically informative sequences were collected around the Indian–SEA region and added to the dataset. These sequences were then used to construct two curated open reading-frame datasets (sequence information is available in Supplementary Table S2); (i) one open reading-frame dataset included 89 sequences, representing 5 from the African lineage and 84 of the Asian lineage; (ii) the other open reading-frame dataset included 84 sequences of Asian lineage only. The datasets were aligned with Mafft v.7.266Citation53, keeping the reading-frame consistent with amino-acid positions, and were visualized and edited in AliViewCitation54. All subsequent analyses were performed with a generalized time-reversible nucleotide substitution model with four gamma distributed rate variation categories and a proportion of invariant sites, as selected by jModeltest v.2Citation55.

To analyze the phylogenies of the 84 Asian lineage strains, including the novel Indian sequence, a Bayesian phylogenetic tree was computed using MrBayes v.3.2.6Citation56 with the dataset including the 5 African ZIKV lineage strains as an out-group. Two parallel runs with four Metropolis-coupled chains were initiated for 5 M Markov chain Monte Carlo generations using the previously determined models of nucleotide evolution with default flat Dirichlet priors, sampling every 1000 generations and discarding the first 25% as burn-in before computing a consensus tree. Since only partial sequences of the Indian isolate were available, it is likely that phylogenetic placement was also examined by constructing three separate trees based on either the capsid (453 bp), envelope (773 bp), or NS2b/NS3 (1393 bp) sequence alignments.

To estimate the evolutionary rates and time to the most recent common ancestor for the Asian lineage, BEAST v.1.8.3Citation57 was employed. Initially, the temporal structure was assessed with TempEst V.1.5.1 (http://tree.bio.ed.ac.uk/software/tempest/) plus MM-type robust regressionCitation58. Subsequently, a path-stone and stepping-stone model-test was performedCitation59, determining that a strict molecular clock with a non-informative continuous-time Markov chain (CTMC) prior and a Bayesian skyline coalescent tree prior with a piecewise-constant demographic model best suited the dataset based on Bayes factor evaluation. Based on the resultant Bayesian phylogenetic tree, which included African ZIKV as outgroup sequences, all Asian ZIKV lineages except the Malaysian 1966-sequence were annotated as a monophyletic clade. The robustness of the dating was also evaluated by excluding the Indian sequence. Using the models and parameters suggested by the model-test, two analyses were run in parallel for 100 M MCMC generations, sampling every 10,000 generations per analysis. The convergence of the two runs was assessed with Tracer 1.6 (http://tree.bio.ed.ac.uk/software/tracer/). Tree- and log-files were combined with LogCombiner (BEAST-packageCitation57), and a maximum-clade credibility tree was computed with TreeAnnotator (BEAST-packageCitation57) after discarding the first 10 M MCMC generations of each run. The resultant consensus tree was visualized and edited in FigTree v.1.4.1 (http://tree.bio.ed.ac.uk/software/figtree/). All computations were run using the CIPRES computational clusterCitation60.

Supplemental material

Supplementary figure 1

Download PDF (257.4 KB)

Supplementary figure 2

Download PDF (292.3 KB)

Supplementary figure 3

Download PDF (250.2 KB)

Supplementary figure 4

Download PDF (305.2 KB)

Supplementary figure 5

Download PDF (277.9 KB)

Supplementary figure 6

Download PDF (275.5 KB)

Supplementary figure 7

Download PDF (240.8 KB)

Supplementary table 1

Download PDF (54.2 KB)

Supplementary table 2

Download PDF (68.7 KB)

Acknowledgements

E.A.G., D.M., J.H.-O.P., M.D.-R., V.-M.C-L., V.E., and X.D.L. are members of the ZIKAlliance Programme of the European Union. M.W.G. and A.K.F. are members of the ZikaPLAN Programme of the European Union. This study was supported by the ZIKAlliance project and the ZikaPLAN project under European Union-Horizon 2020 programme (Grant Agreement Nos. 734548 and 734584, respectively). J.H.-O.P. is supported by the Swedish Research Council FORMAS (Grant Agreement 2015-710).

Conflict of interest

The authors declare that they have no conflict of interest.

Electronic supplementary material

Supplementary Information accompanies this paper at (10.1038/s41426-018-0082-5).

References

  • DickGWKitchenSHaddowAZika Virus (I). Isolations and serological specificityTrans. R. Soc. Trop. Med. Hyg.195246 509 52010.1016/0035-9203(52)90042-4
  • MussoDGublerDJZika VirusClin. Microbiol. Rev.20162948752410.1128/CMR.00072-154861986
  • GouldEAde LamballerieXZanottoPMHolmesECOrigins, evolution, and vector/host coadaptations within the genus FlavivirusAdv. Virus Res.20035927731410.1016/S0065-3527(03)59008-X
  • DuffyMRZika virus outbreak on Yap island, federated states of MicronesiaN. Engl. J. Med.20093602536254310.1056/NEJMoa0805715
  • Livinso T., Carolee M. Situation Report #17: Zika Virus & Dengue Fever, Kosrae State, Federated States of Micronesia (28 Dec 2016). ReliefWeb (2016). https://reliefweb.int/report/micronesia-federated-states/situation-report-17-zika-virus-dengue-fever-kosrae-state accessed 24 Nov 2017.
  • Cao-LormeauVMZika Virus, French Polynesia, South Pacific, 2013Emerg. Infect. Dis.2014201084108610.3201/eid2006.140138
  • PetterssonJHOHow did Zika virus emerge in the Pacific islands and Latin America?mBio20167e012391610.1128/mBio.01239-165061869
  • FariaNRZika virus in the Americas: early epidemiological and genetic findingsScience201635234534910.1126/science.aaf50364918795
  • CauchemezSAssociation between Zika virus and microcephaly in French Polynesia, 2013–15: a retrospective studyLancet20163872125213210.1016/S0140-6736(16)00651-64909533
  • Dupont-RouzeyrolMZika virus evolution on the edges of the Pacific oceanEmerg. Microbes Infect.2017610.1038/emi.2017.1025750458
  • FariaNREstablishment and cryptic transmission of Zika virus in Brazil and the AmericasNature201754640641010.1038/nature224015722632
  • Wang B., et al. The importation of the phylogenetic-transition state of Zika virus to China in 2014. JInfect. 10.1016/j.jinf.2017.09.002 (2017).
  • Sapkal G. N., et al. First laboratory confirmation on the existence of Zika virus disease in India. JInfect. 10.1016/j.jinf.2017.09.020 (2017).
  • PerkasaAIsolation of Zika virus from febrile patient, IndonesiaEmerg. Infect. Dis.20162292492510.3201/eid2205.1519154861529
  • Gu, S. H., et al. Whole-genome sequence analysis of Zika virus, amplified from urine of traveler from the Philippines. Virus Genes. 10.1007/s11262-017-1500-9 (2017).
  • HoZJMOutbreak of Zika virus infection in Singapore: an epidemiological, entomological, virological, and clinical analysisLancet Infect. Dis.20171781382110.1016/S1473-3099(17)30249-9
  • EllisonDWComplete genome sequences of Zika virus strains isolated from the blood of patients in Thailand in 2014 and the Philippines in 2012Genome Announc.20164e003591610.1128/genomeA.00359-164866850
  • TuanHMChikungunya and Zika virus cases detected against a backdrop of endemic dengue transmission in VietnamAm. J. Trop. Med. Hyg.20179714615010.4269/ajtmh.16-09795508909
  • GouldEPetterssonJHiggsSCharrelRde LamballerieXEmerging arboviruses: Why today?One Health2017411310.1016/j.onehlt.2017.06.0015501887
  • SmithburnKCNeutralizing antibodies against certain recently isolated viruses in the sera of human beings residing in East AfricaJ. Immunol.19501952223234
  • MacnamaraFNHornDWPorterfieldJSYellow fever and other arthropod-borne viruses; a consideration of two serological surveys made in South Western NigeriaTrans. R. Soc. Trop. Med. Hyg.19595320221210.1016/0035-9203(59)90072-0
  • MarchetteNJGarciaRRudnickAIsolation of Zika Virus from Aedes aegypti mosquitoes in MalaysiaAm. J. Trop. Med. Hyg.19691841141510.4269/ajtmh.1969.18.411
  • SmithburnKCNeutralizing antibodies against arthropod-borne viruses in the sera of long-time residents of Malaya and BorneoAm. J. Hyg.195459157163
  • SmithburnKCKerrJAGatnePBNeutralizing antibodies against certain viruses in the sera of residents of IndiaJ. Immunol. Baltim. Md19501954248257
  • PondWLArthropod-borne virus antibodies in sera from residents of South-East AsiaTrans. R. Soc. Trop. Med. Hyg.19635736437110.1016/0035-9203(63)90100-7
  • WeaverSCZika virus: History, emergence, biology, and prospects for controlAntivir. Res.2016130698010.1016/j.antiviral.2016.03.0104851879
  • MussoDCao-LormeauVMGublerDJZika virus: Following the path of dengue and chikungunya?Lancet201538624324410.1016/S0140-6736(15)61273-9
  • WiwanitkitVThe current status of Zika virus in Southeast AsiaEpidemiol. Health201638e201602610.4178/epih.e20160264974448
  • RutvisuttinuntWDetection of Zika virus infection in Thailand, 2012–2014Am. J. Trop. Med. Hyg.20159338038310.4269/ajtmh.15-00224530765
  • LorenzCImpact of environmental factors on neglected emerging arboviral diseasesPLoS Negl. Trop. Dis.201711e000595910.1371/journal.pntd.00059595633201
  • VeieraMACeSWest Nile Fever in Brazil: Sporadic case, silent endemic disease or epidemic in its initial stages?Rev. Inst. Med Trop. São Paulo20155727610.1590/S0036-46652015000300017
  • WikanNSmithDRZika virus from a Southeast Asian perspectiveAsian Pac. J. Trop. Med2017101510.1016/j.apjtm.2016.11.013
  • KarabatsosNSupplement to International Catalogue of Arboviruses including certain other viruses of vertebratesAm. J. Trop. Med. Hyg.19782737244010.4269/ajtmh.1978.27.372
  • CareyDEChikungunya and dengue: a case of mistaken identity?J. Hist. Med. Allied Sci.19712624326210.1093/jhmas/XXVI.3.243
  • WeaverSCArrival of chikungunya virus in the new world: prospects for spread and impact on public healthPLoS Negl. Trop. Dis.20148e292110.1371/journal.pntd.00029214072586
  • HiggsSVanlandinghamDChikungunya virus and its mosquito vectorsVector-Borne Zoonotic Dis.20151523124010.1089/vbz.2014.1745
  • MariaATZika virus infections in three travellers returning from South America and the Caribbean respectively, to Montpellier, France, December 2015 to January 2016Eur. Surveill.2016213013110.2807/1560-7917.ES.2016.21.6.30131
  • SmithCEGThe history of dengue in tropical Asia and its probable relationship to the mosquito Aedes aegyptiJ. Trop. Med. Hyg.195659243251
  • OlsonJGKsiazekTGSuhandimanTriwibowoZika virus, a cause of fever in Central Java, IndonesiaTrans. R. Soc. Trop. Med. Hyg.19817538939310.1016/0035-9203(81)90100-0
  • HammonWMSchrackWDSatherGESerological survey for a arthropod-borne virus infections in the PhilippinesAm. J. Trop. Med. Hyg.1958732332810.4269/ajtmh.1958.7.323
  • ChuDTNgocVTNTaoYZika virus infection in Vietnam: current epidemic, strain origin, spreading risk, and perspectiveEur. J. Clin. Microbiol. Infect. Dis.2017362041204210.1007/s10096-017-3030-8
  • DuongVLow circulation of Zika virus, Cambodia, 2007–2016Emerg. Infect. Dis.20172329629910.3201/eid2302.1614325324809
  • FritzRThe unique transmembrane hairpin of flavivirus fusion protein E is essential for membrane fusionJ. Virol.2011854377438510.1128/JVI.02458-103126228
  • YuanLA single mutation in the prM protein of Zika virus contributes to fetal microcephalyScience201735893393610.1126/science.aam7120
  • LiuYEvolutionary enhancement of Zika virus infectivity in Aedes aegypti mosquitoesNature201754548248610.1038/nature223655885636
  • MussoDLanteriMCEmergence of Zika virus: where does it come from and where is it going to?Lancet Infect. Dis.20171725510.1016/S1473-3099(17)30070-1
  • MacNamaraFZika virus: a report on three cases of human infection during an epidemic of jaundice in NigeriaTrans. R. Soc. Trop. Med. Hyg.19544813914510.1016/0035-9203(54)90006-1
  • GrardGZika Virus in Gabon (Central Africa) – 2007: a new threat from Aedes albopictus?PLoS Negl. Trop. Dis.20148e268110.1371/journal.pntd.00026813916288
  • Gake, B. et al. Low seroprevalence of Zika virus in Cameroonian blood donors. Braz. J. Infect. Dis. 10.1016/j.bjid.2017.03.018 (2017).
  • Cao-LormeauVMGuillain-Barré Syndrome outbreak associated with Zika virus infection in French Polynesia: a case-control studyLancet20163871531153910.1016/S0140-6736(16)00562-65444521
  • OehlerEZika virus infection complicated by Guillain-Barré syndrome – case report, French Polynesia, December 2013Eur. Surveill.2014192072010.2807/1560-7917.ES2014.19.9.20720
  • WeaverSCEmergence of epidemic Zika virus transmission and congenital Zika syndrome: are recently evolved traits to blame?mBio20178e020631610.1128/mBio.02063-165225313
  • KatohKStandleyDMMAFFT multiple sequence alignment software version 7: improvements in performance and usabilityMol. Biol. Evol.20133077278010.1093/molbev/mst0103603318
  • LarssonAAliView: a fast and lightweight alignment viewer and editor for large datasetsBioinformatics2014303276327810.1093/bioinformatics/btu5314221126
  • DarribaDTaboadaGLDoalloRPosadaDjModelTest 2: more models, new heuristics and parallel computingNat. Methods2012977210.1038/nmeth.21094594756
  • RonquistFMrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model spaceSyst. Biol.20126153954210.1093/sysbio/sys0293329765
  • DrummondAJSuchardMAXieDRambautABayesian phylogenetics with BEAUti and the BEAST 1.7Mol. Biol. Evol.2012291969197310.1093/molbev/mss0753408070
  • YohaiVJHigh breakdown-point and high efficiency robust estimates for regressionAnn. Stat.19871564265610.1214/aos/1176350366
  • BaeleGLiWLSDrummondAJSuchardMALemeyPAccurate model selection of relaxed molecular clocks in Bayesian phylogeneticsMol. Biol. Evol.20133023924310.1093/molbev/mss243
  • Miller M. A., Pfeiffer W., Schwartz T. in Gateway Computing Environments Workshop (GCE), 1–8 (2010).