For best viewing of the website please use Mozilla Firefox or Google Chrome.
Citation: Tao Hu, Juan Li, Michael J. Carr, Sebastián Duchȇne, Weifeng Shi. The Asian Lineage of Zika Virus: Transmission and Evolution in Asia and the Americas .VIROLOGICA SINICA, 2019, 34(1) : 1-8.

The Asian Lineage of Zika Virus: Transmission and Evolution in Asia and the Americas

  • Corresponding author: Weifeng Shi,, ORCID:
  • Received Date: 29 October 2018
    Accepted Date: 28 November 2018
    Published Date: 25 January 2019
    Available online: 01 February 2019
  • Since first isolation in 1947 from the Zika forest in Uganda, Zika virus (ZIKV) has been principally known as a benign agent associated with sporadic human infections in a restricted number of African countries. However, during 2015–2016, an Asian lineage of ZIKV caused an unprecedentedly large outbreak in the Americas and sizeable numbers of exported cases across the globe. In this review, we critically appraise the recent advances in molecular epidemiological studies of ZIKV performed to date, and we highlight the pivotal role played by genomic surveillance in elucidating the origins, dissemination and evolution of the Asian lineage of ZIKV in Asia and in the Americas.

  • 加载中
    1. Bayer A, Lennemann NJ, Ouyang Y, Bramley JC, Morosky S, Marques ET Jr, Cherry S, Sadovsky Y, Coyne CB (2016) Type iii interferons produced by human placental trophoblasts confer protection against zika virus infection. Cell Host Microbe 19:705-712
        doi: 10.1016/j.chom.2016.03.008

    2. Besnard M, Eyrolle-Guignot D, Guillemette-Artur P, Lastere S, BostBezeaud F, Marcelis L, Abadie V, Garel C, Moutard ML, Jouannic JM, Rozenberg F, Leparc-Goffart I, Mallet HP (2016) Congenital cerebral malformations and dysfunction in fetuses and newborns following the 2013 to 2014 zika virus epidemic in french polynesia. Euro Surveill 21:30181
        doi: 10.2807/1560-7917.ES.2016.21.13.30181

    3. Bouckaert R, Heled J, Kuhnert D, Vaughan T, Wu CH, Xie D, Suchard MA, Rambaut A, Drummond AJ (2014) Beast 2: a software platform for bayesian evolutionary analysis. PLoS Comput Biol 10:e1003537
        doi: 10.1371/journal.pcbi.1003537

    4. Broutet N, Krauer F, Riesen M, Khalakdina A, Almiron M, Aldighieri S, Espinal M, Low N, Dye C (2016) Zika virus as a cause of neurologic disorders. N Engl J Med 374:1506-1509
        doi: 10.1056/NEJMp1602708

    5. Calvet GA, Santos FB, Sequeira PC (2016) Zika virus infection: epidemiology, clinical manifestations and diagnosis. Curr Opin Infect Dis 29:459-466
        doi: 10.1097/QCO.0000000000000301

    6. Cao-Lormeau VM, Roche C, Teissier A, Robin E, Berry AL, Mallet HP, Sall AA, Musso D (2014) Zika virus, french polynesia, south pacific, 2013. Emerg Infect Dis 20:1085-1086

    7. Cao-Lormeau VM, Blake A, Mons S, Lastere S, Roche C, Vanhomwegen J, Dub T, Baudouin L, Teissier A, Larre P, Vial AL, Decam C, Choumet V, Halstead SK, Willison HJ, Musset L, Manuguerra JC, Despres P, Fournier E, Mallet HP, Musso D, Fontanet A, Neil J, Ghawche F (2016) Guillain-barre syndrome outbreak associated with zika virus infection in french polynesia: a case-control study. Lancet 387:1531-1539
        doi: 10.1016/S0140-6736(16)00562-6

    8. Delatorre E, Fernandez J, Bello G (2018) Investigating the role of easter island in migration of zika virus from south pacific to americas. Emerg Infect Dis 24:2119-2121
        doi: 10.3201/eid2411.180586

    9. Driggers RW, Ho CY, Korhonen EM, Kuivanen S, Jaaskelainen AJ, Smura T, Rosenberg A, Hill DA, DeBiasi RL, Vezina G, Timofeev J, Rodriguez FJ, Levanov L, Razak J, Iyengar P, Hennenfent A, Kennedy R, Lanciotti R, du Plessis A, Vapalahti O (2016) Zika virus infection with prolonged maternal viremia and fetal brain abnormalities. N Engl J Med 374:2142-2151
        doi: 10.1056/NEJMoa1601824

    10. Drummond AJ, Suchard MA, Xie D, Rambaut A (2012) Bayesian phylogenetics with beauti and the beast 1.7. Mol Biol Evol 29:1969-1973
        doi: 10.1093/molbev/mss075

    11. Duffy MR, Chen TH, Hancock WT, Powers AM, Kool JL, Lanciotti RS, Pretrick M, Marfel M, Holzbauer S, Dubray C, Guillaumot L, Griggs A, Bel M, Lambert AJ, Laven J, Kosoy O, Panella A, Biggerstaff BJ, Fischer M, Hayes EB (2009) Zika virus outbreak on yap island, federated states of micronesia. N Engl J Med 360:2536-2543
        doi: 10.1056/NEJMoa0805715

    12. Faria NR, Azevedo R, Kraemer MUG, Souza R, Cunha MS, Hill SC, Thézé J, Bonsall MB, Bowden TA, Rissanen I et al (2016) Zika virus in the americas: early epidemiological and genetic findings. Science 352:345-349
        doi: 10.1126/science.aaf5036

    13. Faria NR, Quick J, Claro IM, Thézé J, de Jesus JG, Giovanetti M, Kraemer MUG, Hill SC, Black A, da Costa AC et al (2017) Establishment and cryptic transmission of zika virus in brazil and the americas. Nature 546:406-410
        doi: 10.1038/nature22401

    14. Gong Z, Gao Y, Han GZ (2016) Zika virus: two or three lineages? Trends Microbiol 24:521-522
        doi: 10.1016/j.tim.2016.05.002

    15. Grubaugh ND, Ladner JT, Kraemer MUG, Dudas G, Tan AL, Gangavarapu K, Wiley MR, White S, Thézeé J, Magnani DM et al (2017) Genomic epidemiology reveals multiple introductions of zika virus into the united states. Nature 546:401-405
        doi: 10.1038/nature22400

    16. Gubler DJ, Vasilakis N, Musso D (2017) History and emergence of zika virus. J Infect Dis 216:S860-S867
        doi: 10.1093/infdis/jix451

    17. Hamel R, Liegeois F, Wichit S, Pompon J, Diop F, Talignani L, Thomas F, Despres P, Yssel H, Misse D (2016) Zika virus: epidemiology, clinical features and host-virus interactions. Microbes Infect 18:441-449
        doi: 10.1016/j.micinf.2016.03.009

    18. Hancock WT, Marfel M, Bel M (2014) Zika virus, french polynesia, south pacific, 2013. Emerg Infect Dis 20:1960
        doi: 10.3201/eid2011.141380

    19. Ikejezie J, Shapiro CN, Kim J, Chiu M, Almiron M, Ugarte C, Espinal MA, Aldighieri S (2017) Zika virus transmission-region of the americas, may 15, 2015-December 15, 2016. Am J Transplant 17:1681-1686
        doi: 10.1111/ajt.2017.17.issue-6

    20. Ioos S, Mallet HP, Leparc Goffart I, Gauthier V, Cardoso T, Herida M (2014) Current zika virus epidemiology and recent epidemics. Med Mal Infect 44:302-307
        doi: 10.1016/j.medmal.2014.04.008

    21. Jagger BW, Miner JJ, Cao B, Arora N, Smith AM, Kovacs A, Mysorekar IU, Coyne CB, Diamond MS (2017) Gestational stage and ifn-lambda signaling regulate zikv infection in utero. Cell Host Microbe 22:366-376
        doi: 10.1016/j.chom.2017.08.012

    22. Kenney AD, Dowdle JA, Bozzacco L, McMichael TM, St Gelais C, Panfil AR, Sun Y, Schlesinger LS, Anderson MZ, Green PL, Lopez CB, Rosenberg BR, Wu L, Yount JS (2017) Human genetic determinants of viral diseases. Annu Rev Genet 51:241-263
        doi: 10.1146/annurev-genet-120116-023425

    23. Lanciotti RS, Kosoy OL, Laven JJ, Velez JO, Lambert AJ, Johnson AJ, Stanfield SM, Duffy MR (2008) Genetic and serologic properties of zika virus associated with an epidemic, yap state, micronesia, 2007. Emerg Infect Dis 14:1232-1239
        doi: 10.3201/eid1408.080287

    24. Lednicky J, Beau De Rochars VM, El Badry M, Loeb J, Telisma T, Chavannes S, Anilis G, Cella E, Ciccozzi M, Rashid M, Okech B, Salemi M, Morris JG Jr (2016) Zika virus outbreak in haiti in 2014: molecular and clinical data. PLoS Negl Trop Dis 10:e0004687
        doi: 10.1371/journal.pntd.0004687

    25. Lim SK, Lim JK, Yoon IK (2017) An update on zika virus in asia. Infect Chemother 49:91-100
        doi: 10.3947/ic.2017.49.2.91

    26. Liu Y, Liu J, Du S, Shan C, Nie K, Zhang R, Li XF, Zhang R, Wang T, Qin CF, Wang P, Shi PY, Cheng G (2017) Evolutionary enhancement of zika virus infectivity in aedes aegypti mosquitoes. Nature 545:482-486
        doi: 10.1038/nature22365

    27. Liu ZY, Shi WF, Qin CF (2019) The evolution of Zika virus from Asia to the Americas. Nat Rev Microbiol (in press)

    28. Marchette NJ, Garcia R, Rudnick A (1969) Isolation of zika virus from aedes aegypti mosquitoes in malaysia. Am J Trop Med Hyg 18:411-415
        doi: 10.4269/ajtmh.1969.18.411

    29. Metsky HC, Matranga CB, Wohl S, Schaffner SF, Freije CA, Winnicki SM, West K, Qu J, Baniecki ML, Gladden-Young A et al (2017) Zika virus evolution and spread in the americas. Nature 546:411-415
        doi: 10.1038/nature22402

    30. Mlakar J, Korva M, Tul N, Popovic M, Poljsak-Prijatelj M, Mraz J, Kolenc M, Resman Rus K, Vesnaver Vipotnik T, Fabjan Vodusek V, Vizjak A, Pizem J, Petrovec M, Avsic Zupanc T (2016) Zika virus associated with microcephaly. N Engl J Med 374:951-958
        doi: 10.1056/NEJMoa1600651

    31. Moi ML, Nguyen TTT, Nguyen CT, Vu TBH, Tun MMN, Pham TD, Pham NT, Tran T, Morita K, Le TQM, Dang DA, Hasebe F (2017) Zika virus infection and microcephaly in vietnam. Lancet Infect Dis 17:805-806
        doi: 10.1016/S1473-3099(17)30412-7

    32. Musso D, Gubler DJ (2016) Zika virus. Clin Microbiol Rev 29:487-524
        doi: 10.1128/CMR.00072-15

    33. Netto EM, Moreira-Soto A, Pedroso C, Hoser C, Funk S, Kucharski AJ, Rockstroh A, Kummerer BM, Sampaio GS, Luz E et al (2017) High zika virus seroprevalence in salvador, northeastern brazil limits the potential for further outbreaks. MBio 8:e01390-17

    34. O'Brien TR, Prokunina-Olsson L, Donnelly RP (2014) Ifn-lambda4: the paradoxical new member of the interferon lambda family. J Interferon Cytokine Res 34:829-838
        doi: 10.1089/jir.2013.0136

    35. Oehler E, Watrin L, Larre P, Leparc-Goffart I, Lastere S, Valour F, Baudouin L, Mallet H, Musso D, Ghawche F (2014) Zika virus infection complicated by guillain-barre syndrome-case report, french polynesia, December 2013. Euro Surveill 19:20720

    36. Olson JG, Ksiazek TG, Suhandiman Triwibowo (1981) Zika virus, a cause of fever in central java, indonesia. Trans R Soc Trop Med Hyg 75:389-393
        doi: 10.1016/0035-9203(81)90100-0

    37. Pettersson JH, Eldholm V, Seligman SJ, Lundkvist A, Falconar AK, Gaunt MW, Musso D, Nougairede A, Charrel R, Gould EA, de Lamballerie X (2016) How did zika virus emerge in the pacific islands and latin america? MBio 7:e01239-16

    38. Rossi SL, Ebel GD, Shan C, Shi PY, Vasilakis N (2018) Did zika virus mutate to cause severe outbreaks? Trends Microbiol 26:877-885
        doi: 10.1016/j.tim.2018.05.007

    39. Shen S, Shi J, Wang J, Tang S, Wang H, Hu Z, Deng F (2016) Phylogenetic analysis revealed the central roles of two african countries in the evolution and worldwide spread of zika virus. Virol Sin 31:118-130
        doi: 10.1007/s12250-016-3774-9

    40. Shi W, Zhang Z, Ling C, Carr MJ, Tong Y, Gao GF (2016) Increasing genetic diversity of zika virus in the latin american outbreak. Emerg Microbes Infect 5:e68

    41. Singapore Zika Study G (2017) Outbreak of zika virus infection in singapore: an epidemiological, entomological, virological, and clinical analysis. Lancet Infect Dis 17:813-821
        doi: 10.1016/S1473-3099(17)30249-9

    42. Siraj AS, Perkins TA (2017) Assessing the population at risk of zika virus in asia - is the emergency really over? BMJ Glob Health 2:e000309
        doi: 10.1136/bmjgh-2017-000309

    43. Song BH, Yun SI, Woolley M, Lee YM (2017) Zika virus: history, epidemiology, transmission, and clinical presentation. J Neuroimmunol 308:50-64
        doi: 10.1016/j.jneuroim.2017.03.001

    44. Stadler T, Kuhnert D, Bonhoeffer S, Drummond AJ (2013) Birthdeath skyline plot reveals temporal changes of epidemic spread in hiv and hepatitis c virus (hcv). Proc Natl Acad Sci U S A 110:228-233
        doi: 10.1073/pnas.1207965110

    45. Thézé J, Li T, du Plessis L, Bouquet J, Kraemer MUG, Somasekar S, Yu G, de Cesare M, Balmaseda A, Kuan G et al (2018) Genomic epidemiology reconstructs the introduction and spread of zika virus in central america and mexico. Cell Host Microbe 23:855-864.e7
        doi: 10.1016/j.chom.2018.04.017

    46. Wang L, Valderramos SG, Wu A, Ouyang S, Li C, Brasil P, Bonaldo M, Coates T, Nielsen-Saines K, Jiang T, Aliyari R, Cheng G (2016) From mosquitos to humans: genetic evolution of zika virus. Cell Host Microbe 19:561-565
        doi: 10.1016/j.chom.2016.04.006

    47. Wilder-Smith A, Chang CR, Leong WY (2018) Zika in travellers 1947-2017: a systematic review. J Travel Med 25:tay044.
        doi: 10.1093/jtm/tay044

    48. Wongsurawat T, Athipanyasilp N, Jenjaroenpun P, Jun SR, Kaewnapan B, Wassenaar TM, Leelahakorn N, Angkasekwinai N, Kantakamalakul W, Ussery DW, Sutthent R, Nookaew I, Horthongkham N (2018) Case of microcephaly after congenital infection with asian lineage zika virus, thailand. Emerg Infect Dis 24:1758.
        doi: 10.3201/eid2409.180416

    49. Worobey M (2017) Epidemiology: molecular mapping of zika spread. Nature 546:355-357
        doi: 10.1038/nature22495

    50. Ye Q, Liu ZY, Han JF, Jiang T, Li XF, Qin CF (2016) Genomic characterization and phylogenetic analysis of zika virus circulating in the americas. Infect Genet Evol 43:43-49
        doi: 10.1016/j.meegid.2016.05.004

    51. Yuan L, Huang XY, Liu ZY, Zhang F, Zhu XL, Yu JY, Ji X, Xu YP, Li G, Li C, Wang HJ, Deng YQ, Wu M, Cheng ML, Ye Q, Xie DY, Li XF, Wang X, Shi W, Hu B, Shi PY, Xu Z, Qin CF (2017) A single mutation in the prm protein of zika virus contributes to fetal microcephaly. Science 358:933-936
        doi: 10.1126/science.aam7120

    52. Zhang Y, Chen W, Wong G, Bi Y, Yan J, Sun Y, Chen E, Yan H, Lou X, Mao H, Xia S, Gao GF, Shi W, Chen Z (2016) Highly diversified zika viruses imported to china, 2016. Protein Cell 7:461-464
        doi: 10.1007/s13238-016-0274-5

  • 加载中


Article Metrics

Article views(7115) PDF downloads(60) Cited by()

Proportional views

    The Asian Lineage of Zika Virus: Transmission and Evolution in Asia and the Americas

      Corresponding author: Weifeng Shi,
    • 1. Key Laboratory of Etiology and Epidemiology of Emerging Infectious Diseases in Universities of Shandong, Taishan Medical College, Taian 271000, China
    • 2. Global Station for Zoonosis Control, Global Institution for Collaborative Research and Education (GI-CoRE), Hokkaido University, Sapporo 001-0020, Japan
    • 3. National Virus Reference Laboratory, School of Medicine, University College Dublin, Belfield, Dublin 4, Ireland
    • 4. Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Parkville, VIC 3020, Australia

    Abstract: Since first isolation in 1947 from the Zika forest in Uganda, Zika virus (ZIKV) has been principally known as a benign agent associated with sporadic human infections in a restricted number of African countries. However, during 2015–2016, an Asian lineage of ZIKV caused an unprecedentedly large outbreak in the Americas and sizeable numbers of exported cases across the globe. In this review, we critically appraise the recent advances in molecular epidemiological studies of ZIKV performed to date, and we highlight the pivotal role played by genomic surveillance in elucidating the origins, dissemination and evolution of the Asian lineage of ZIKV in Asia and in the Americas.

    • Zika virus (ZIKV), a mosquito-borne flavivirus, has been known for over 70 years. The positive sense, singlestranded RNA genomic structure of ZIKV and processing of the single encoded large polyprotein is characteristic of other flaviviruses, such as yellow fever virus and dengue virus. The majority of ZIKV infections are asymptomatic; however, in some cases can have a febrile presentation with mild clinical symptoms, such as, headache, skin rashes and conjunctivitis, which are typically self-limiting, resolve without consequence and do not require anti-viral treatment (Calvet et al. 2016; Hamel et al. 2016; Song et al. 2017). Notably, prior to 2007, a distinguishing feature from other closely related flaviviruses able to cause large scale outbreaks, was that ZIKV had been only reported to be associated with sporadic human infections in a small number of African countries (Calvet et al. 2016; Hamel et al. 2016; Song et al. 2017). In 2007, for the first time, ZIKV was linked to a large-scale epidemic on Yap Island in Micronesia, with a total of at least 49 laboratory-confirmed cases (Lanciotti et al. 2008; Duffy et al. 2009). ZIKV then caused an even larger outbreak in French Polynesia and surrounding islands in 2013 (Cao-Lormeau et al. 2014; Hancock et al. 2014; Ioos et al. 2014). It was in this outbreak that ZIKV was first linked to Guillain-Barrex syndrome (GBS), a severe autoimmune neuropathology (Oehler et al. 2014; Cao-Lormeau et al. 2016). Surprisingly, ZIKV subsequently disseminated across the Americas, where it caused an outbreak on an unprecedented scale for this previously considered benign agent. By the end of 2016, at least 48 countries and territories in the Americas had reported more than 175, 000 confirmed ZIKV cases (Ikejezie et al. 2017). A retrospective investigation has revealed that ZIKV-associated microcephaly was first reported in French Polynesia during the 2013-2014 epidemic (Besnard et al. 2016). More importantly, ZIKV was first definitively associated with microcephaly in the South American outbreak: a condition where a neonate has a significantly smaller head circumference than normal, or the cranium ceases normal development after birth (Broutet et al. 2016; Driggers et al. 2016; Mlakar et al. 2016). By February 2017, a total of 2654 cases of ZIKV-related fetal microcephaly and central nervous system (CNS) malfunction were reported to the World Health Organization, with approximately 90% of the cases distributed in Brazil (

      In the preceding 3 years, there have been numerous studies describing the molecular epidemiologic and genomic characterization of ZIKV identified from the 2015–2016 outbreak. The first important study that examined the ZIKV-associated American outbreak was published in Science where Faria and colleagues dated the transmission to have occurred via a single introduction event from Pacific Islands between May and December of 2013 (Faria et al. 2016). Three seminal studies were then subsequently published in the same issue of Nature in June of 2017 which delineated the transmission and evolution of ZIKV across the Americas (Faria et al. 2017; Grubaugh et al. 2017; Metsky et al. 2017; Worobey 2017). Faria et al. revealed that Bahia province in northeastern Brazil was the epicenter of the outbreak in the Americas, where ZIKV first became established and then subsequently disseminated nationally and internationally (Faria et al. 2017). Metsky and colleagues determined when the virus was transmitted to Central America and the Caribbean islands (Metsky et al. 2017). Grubaugh and co-workers further analyzed when and how often the virus was disseminated into the USA and pointed out that ZIKV had most likely been imported multiple times from the Caribbean (Grubaugh et al. 2017). Subsequently, Thézé and colleagues reported that ZIKV was imported into Central America and Mexico via multiple independent introduction events. In particular, one introduction, likely from Brazil via Honduras, led to most infections in Central America and Mexico from late 2014 onwards (Thézé et al. 2018). In this mini-review, we critically appraise the current state of knowledge with regards to the transmission and evolution of the Asian lineage of ZIKV in the unprecedentedly large epidemic in the Americas and highlight the role of Southeast Asia in the evolution of this agent.

    • Previous studies have shown that ZIKV has diverged into two principal genetic lineages: African and Asian (Faria et al. 2016; Pettersson et al. 2016; Zhang et al. 2016). Some researchers proposed a three-lineage classification nomenclature where the Asian lineage was further divided into two separable sub-lineages (Gubler et al. 2017) or, alternatively, the African strains were further divided into two sub-lineages based on phylogenetic analysis of partial genomic sequences (Shen et al. 2016). The later study also reported an African sister group to all currently known ZIKV major lineages (Shen et al. 2016) which suggested that the Asian lineage originated from the African lineage (Gong et al. 2016). This interesting finding clearly warrants further investigation and more comprehensive sampling and full genome sequencing.

      The prototypic strain of the ZIKV Asian lineage was first isolated in 1966 from mosquito pools collected in Bentong, Malaysia (Marchette et al. 1969). Although serological studies revealed that the first documented ZIKV outbreak in Asia may have occurred in central Java, Indonesia in 1977 (Olson et al. 1981) and ZIKV infections are possibly more widespread and may be circulating asymptomatically in Southeast and South Asia (Musso and Gubler 2016), only very limited partial genomic sequences are available. However, in 2007, the first documented large-scale outbreak of ZIKV was reported on a remote Pacific island, Yap in Micronesia (Lanciotti et al. 2008; Duffy et al. 2009). Complete genome sequencing showed, conclusively, that the outbreak arose from the introduction of a ZIKV belonging to the Asian lineage (Lanciotti et al. 2008). During 2013–2014, the virus caused outbreaks in four other Pacific islands with the outbreak in French Polynesia resulting in thousands of suspected human infections (Cao-Lormeau et al. 2014; Hancock et al. 2014; Ioos et al. 2014). Phylogenetic analysis revealed that the ZIKV identified from laboratory-confirmed cases during this outbreak also belonged to the Asian lineage, which potentially originated in Southeast Asia (Cao-Lormeau et al. 2014). However, these two outbreaks in the Pacific likely may have arisen from two independent introductions events from Southeast Asia (Delatorre et al. 2018), or another unsampled Asian location, because they clustered together with strains from Southeast Asia in different positions of the tree (Fig. 1). Later in 2016, Chinese travelers returning from Fiji and Samoa were confirmed to be infected by ZIKV (Zhang et al. 2016). Genomic sequencing showed that these ZIKVs did not cluster with those responsible for the outbreak in the Americas (Shi et al. 2016). Alternatively, they formed a separate cluster outside the ZIKV Asian lineage in the Americas (Shi et al. 2016), and strains of these Oceanian cases likely represent ZIKVs circulating locally within the Pacific which diverged from those responsible for the 2013–2014 outbreak.

      Figure 1.  Phylogeny and molecular characterization of the Asian lineage of ZIKV. Full-length ZIKV genomes (n = 448) derived from the Asian lineage were assembled and a phylogenetic analysis using BEAST 1.8.4 was performed (Drummond et al. 2012). The GTR + Γ nucleotide model was applied to account for rate variation among sites, with four categories for the Γ distribution. A log-normal distributed relaxed molecular clock was employed and the Bayesian skyline coalescent was designated as the tree prior. A Markov chain Monte Carlo with one hundred million steps was run twice independently and the first 10% were removed as burn-in. The posterior probability distribution of the parameters and trees from the two independent runs was combined using LogCombiner embedded in BEAST 1.8.4. Viruses identified from the ZIKV African lineage and from Indonesia cases in 1966 were also downloaded and phylogenetically analyzed using Maximum Likelihood. In the African lineage, one strain, ArD142623 (GenBank No. KF3838120), had a long branch length, which is displayed using a dashed line. In the Asian lineage, tips with circles represent confirmed ZIKV imported cases. Line color represents the collection site of the virus and the classification of the countries into different regions was based on Wikipedia ( In the right panel, different colors represent different amino acids, which have been shown in the legend in the bottom right of the panel.

      It is noteworthy that there were several ZIKVs identified from Haiti during an outbreak of chikungunya fever in the Caribbean in mid-2014 (e.g. GenBank Accession Nos. KY415986–KY415991), which also segregated within the American ZIKV cases (Lednicky et al. 2016). However, the first confirmed ZIKV infection was reported in northeastern Brazil in February of 2015, the presumed epicenter of the American outbreak (Faria et al. 2016, 2017). From the tree, we can see that shortly after ZIKV was imported into the Americas, it diverged into four major sub-lineages in a relatively short period of time between late 2013 and early 2014, probably in Brazil as strains from Brazil were always located basally (Fig. 1). The majority of viruses within sub-lineage A were from South America. Sub-lineage B was also mainly from South America; however, some viruses were transmitted to the Caribbean regions. Most of the viruses in sub-lineage C were from both Central and North America, and, finally, sub-lineage D was populated with viruses from the Caribbean and North America. Therefore, although dissemination of the virus from Brazil to Central and North America, and to the Caribbean occurred independently many times (Fig. 1), the predominant ZIKVs circulating in Central and North America and the Caribbean appear to have been transmitted via a small number of dominant importation events.

      We estimated the effective reproductive number (Re; i.e. the average number of secondary infections) of ZIKV over time using a birth–death skyline approach (Stadler et al. 2013) implemented in BEAST 2 (Bouckaert et al. 2014). There were three major phases in which ZIKV rapidly spread with Re > 1 (Fig. 2). The first phase was estimated to have occurred between March and November of 2013 when the virus caused the large-scale outbreaks on the Pacific islands (Fig. 2). Interestingly, this is also the estimated time when ZIKV was introduced into South America (Faria et al. 2016) and the time when the major South American sub-lineages diverged (Fig. 1). Therefore, the rapid increase in the average number of transmissions in the Pacific might have facilitated, or even been a prerequisite for this virus transmission event. The second phase where ZIKV population size expanded was between March 2014 and March 2015 (Fig. 2). After the virus was imported into South America, it evolved under an Re > 1 for approximately one year; then, towards the end of this phase, the first confirmed ZIKV infection case was reported in northeastern Brazil (Faria et al. 2016, 2017). More importantly, it is during this phase that many early virus dissemination events from Brazil to Central America and the Caribbean might have occurred, although evidence is scant due to undersampling (Faria et al. 2017; Metsky et al. 2017). Following phase 2, there was a short period of decrease in Re. We speculate that this might be attributable to decreases in population numbers of competent mosquito vectors for ZIKV transmission, as this time period overlaps with the lowest annual ambient temperatures in the southern hemisphere winter. However, even during this phase, ZIKV may have also circulated continuously and undetected in tropical Central America and the Caribbean (Thézé et al. 2018). The third phase began from approximately August of 2015 and ended in late 2015 (Fig. 2). During this phase, the virus continued to circulate in the Americas. It is thought that after this phase the virus was imported into the United States and led to the first autochthonous transmission events in Florida (Grubaugh et al. 2017). ZIKV laboratory-confirmed events and the cases of congenital malformations subsequently decreased in South America in 2016 and may, at least be partly attributable to the decrease in the number of susceptible individuals to maintain chains of transmission after such an unprecedently large outbreak; however, the proximate and ultimate causes are likely complex and clearly warrant further investigation.

      Figure 2.  Effective reproductive number (Re) of ZIKV of the Asian lineage through time. The effective reproductive number (Re) was estimated using a birth–death skyline approach implemented in BEAST 2 with ZIKVs identified since 2013. An uncorrelated lognormal relaxed clock model and the GTR + Γ nucleotide substitution model were used. A Markov chain Monte Carlo was run with one hundred million steps and the first 10% were removed as burn-in. In the serially sampled birth–death skyline model, the priors for the molecular clock rate, Re, the sampling proportion, and the rate at which patients recover were set according to a previous study, which were also supported by an alternative approach allowing for sampling heterogeneity (data not shown) (Thézé et al. 2018). The red line and lighter shading represented the median posterior estimate of Re and its 95% highest posterior density credible intervals, respectively. This figure was created using the R-package bdskytools (available at

    • Between 2014 and 2016, ZIKV outbreaks and sporadic autochthonous transmission events were also reported in several Southeast and South Asian countries, including Bangladesh, Cambodia, India, Indonesia, Lao People's Democratic Republic, Malaysia, Maldives, Myanmar, Philippines, Singapore, Thailand and Vietnam (Lim et al. 2017), suggesting the widespread circulation of ZIKV in Southeast and South Asia. Of particular note is the ZIKV outbreak in Singapore starting from August of 2016, with a total of 455 reported cases between August 27 and November 30, 2016 (Singapore Zika Study 2017). Phylogenetic analysis showed that ZIKV strains identified in the Singaporean outbreak were not descended from those circulating in the Americas. Alternatively, they were suggested to diverge from ZIKVs circulating locally in Southeast Asia (Singapore Zika Study 2017). In addition, Thailand reported the largest number of confirmed ZIKV cases (~ 700) in Asia and more than 200 ZIKV cases were reported in Vietnam in 2016 (Lim et al. 2017). Apart from these three countries with hundreds of ZIKV cases, few sporadic autochthonous infections were reported in the remaining countries mentioned above. In particular, though far fewer in number, at least three suspected microcephaly cases (two from Thailand and one from Vietnam) were reported (Lim et al. 2017; Moi et al. 2017; Wongsurawat et al. 2018). Precisely why congenital malformations appear more frequently in the South American cases than Asian cases is unclear but may be attributable to host genetic factors conferring differing levels of susceptibility to infection (Bayer et al. 2016; Jagger et al. 2017; Kenney et al. 2017; O'Brien et al. 2014). Genome wide association studies could potentially be employed to identify single nucleotide polymorphisms responsible for this phenomenon. Viral genomic sequencing revealed that ZIKV from the two microcephaly cases from Thailand belonged to the Asian lineage but did not cluster with the ZIKVs responsible for the outbreak in the Americas (Wongsurawat et al. 2018). Therefore, virus-specific adaptations may also play a role in the differing pathogenesis which we explore in more detail below.

    • The unprecedentedly large 2015–2016 ZIKV outbreak was mainly centered upon the Americas. Contemporaneously, however, there were also thousands of imported ZIKV cases outside the Americas. For example, more than two thousand ZIKV cases were imported into Europe by travelers returning from ZIKV-endemic regions (Wilder-Smith et al. 2018). Hypothetically, the number of imported ZIKV cases in Europe and other regions may be greatly underestimated due to a number of factors. These include the high percentage of sub-clinical infections, the overlapping clinical presentation of Zika with numerous acute-onset, travel-associated, influenza-like illnesses, the narrow window and low titer viremia precluding definitive molecular confirmation by amplification-based detection approaches, and, finally, the paucity of flavivirus-specific serological assays allowing the unambiguous discrimination from other related agents (Lim et al. 2017). Notably, many highly populous countries in Asia have also reported imported ZIKV cases, including China, India and Japan (Lim et al. 2017), and greater than one billion Asian people live in areas with competent mosquito vectors and are therefore at high risk for large Zika outbreaks (Lim et al. 2017; Siraj and Perkins 2017). Hypothetically, certain human populations may be less susceptible to ZIKV and the virus may circulate silently without detection in large population centers further adapting to more effectively transmit among humans and, thus, represent reservoirs for the origin of the founder viruses in large-scale outbreaks in immunologically naïve populations elsewhere with greater numbers of susceptible individuals. Also, whether the hyperendemicity of the four serologically distinct serotypes of dengue viruses in Southeast Asian countries impacts upon the transmission of ZIKV is unclear but may conceivably confer differing susceptibility to ZIKV which warrants further study. It should be noted that apart from the Americas, Southeast Asia has also been found, on multiple occasions, to export ZIKV to the Pacific islands and East Asia. Phylogenetic analysis revealed that these imported ZIKVs scattered across the Asian lineage, without clear geographic clustering (Fig. 1). However, the long-distance dissemination of the virus should be given special attention as evidenced by the finding that the progenitor strains responsible for the American ZIKV outbreak were likely transmitted from the Pacific via a single importation event (Faria et al. 2016).

    • It has been controversial whether the Asian or African lineages of ZIKV are intrinsically more virulent (Rossi et al. 2018), although strains of the Asian lineage are responsible for all known ZIKV outbreaks. However, different researchers have shown that there were distinct amino acid differences between ZIKV strains from the African and Asian lineages (Fig. 1) (Pettersson et al. 2016; Ye et al. 2016) and, also, between strains isolated from humans and mosquitoes (Wang et al. 2016), some of which have been experimentally verified to have important roles in ZIKV transmission and pathogenesis. For example, the S139N substitution in the prM protein has been found to significantly increase the neurovirulence of ZIKV, which might, as a consequence, have caused an excess of microcephaly cases and other congenital CNS malformations during the American outbreak (Yuan et al. 2017). In addition, the A982V substitution in the NS1 protein was reported to increase NS1 antigenemia, which was suggested to have enhanced the ZIKV infectivity and prevalence in mosquitoes and may have contributed to the widespread transmission of ZIKV in the Americas (Liu et al. 2017). The evolutionary pattern of these biologically critical amino acid positions in ZIKV has been illustrated elsewhere (Liu et al. 2019).

      Apart from the previously reported substitutions that distinguish the two main ZIKV lineages, we have also identified a panel of sub-lineage-specific amino acid substitutions (Fig. 2), including residues 109, 709, 1118 and 3144 in the Oceanian origin cases imported to China from returning travelers, and 107, 1118 and 3353 in sub-lineage A, 894, 2074 and 3045 in sub-lineage C, and 2842 and 3398 in sub-lineage D. These sub-lineage-specific amino acid substitutions may be without phenotypic consequence and thus have a neutral effect on viral fitness if they arise from a founder effect; however, they represent the further diversifying evolution of ZIKV in the Asian lineage which may contribute to a differing propensity for transmission and/or pathogenesis.

      In summary, although the ZIKV outbreak in the Americas ended suddenly (Netto et al. 2017), just like its unexpected beginning, genomic sequencing has elucidated the origins, transmission, evolution and factors underlying the differing pathogenic properties of ZIKV in the Americas and in Southeast Asia. Although recent advances have revealed the potential biological significance of several specific amino acid substitutions, these mutations should be further verified employing animal models and the biological functions of other Asian lineage-specific amino acid substitutions should also be explored. Considering the extensively distributed nature of mosquito vector species which are competent for ZIKV transmission across Southeast Asia, India and southern China, previous outbreaks caused by ZIKV and other related flaviviruses, the long-term (likely asymptomatic) circulation and adaptation to humans, the underestimated genetic diversity of ZIKV in this region due to inadequate sampling, and the ongoing ZIKV epidemic in India ( and the recent Chinese imported case from the Maldives, a comprehensive survey of ZIKV and other flaviviruses in competent vectors in the Eastern Pacific region is urgently required. This would allow a finer scale delineation of the extant genetic diversity and detailed pathogenesis studies and, critically, to facilitate a rational, evidence-based development of vaccines to stable epitope targets to mitigate further disease outbreaks.

    • This study was supported by the National Key Science and Technology Projects of China (2017ZX10104001-006) and the National Science and Technology Major Project (2018ZX10101004-002). W.S. was supported by the Taishan Scholars program of Shandong Province (ts201511056).

    • The authors declare that they have no conflict of interest.

    • The authors declare that they have no conflict of interest. This article does not contain any studies with human or animal subjects performed by any of the authors.

    • This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

    Figure (2)  Reference (52) Relative (20)



    DownLoad:  Full-Size Img  PowerPoint