On 6 June 2017, the World Health Organization (WHO) published updates to its ‘Essential Medicines List’ (EML). Read more here.

Extended deadline (from 1 July to 31 July) 2017 for call to submit papers on effectiveness and cost-effectiveness of screening and prevention of infectious diseases among newly arrived migrants in Europe. Read more here.

Eurosurveillance is on the updated list of the Directory of Open Access Journals and in the SHERPA/RoMEO database. Read more here.

Follow Eurosurveillance on Twitter: @Eurosurveillanc

In this issue

Home Eurosurveillance Edition  2016: Volume 21/ Issue 38 Article 2
Back to Table of Contents
Previous Download (pdf)

Eurosurveillance, Volume 21, Issue 38, 22 September 2016
Rapid communication
Maurer-Stroh, Mak, Ng, Phuah, Huber, Marzinek, Holdbrook, Lee, Cui, and Lin: South-east Asian Zika virus strain linked to cluster of cases in Singapore, August 2016

+ Author affiliations

Citation style for this article: Maurer-Stroh S, Mak T, Ng Y, Phuah S, Huber RG, Marzinek JK, Holdbrook DA, Lee RT, Cui L, Lin RT. South-east Asian Zika virus strain linked to cluster of cases in Singapore, August 2016. Euro Surveill. 2016;21(38):pii=30347. DOI:

Received:19 September 2016; Accepted:22 September 2016

Outbreak detection

On 22 August 2016, the Ministry of Health (MOH), Singapore, was informed by a general practitioner of a spate of cases presenting with non-specific symptoms, including rash, polyarthralgia and low grade fever. In addition, a number of cases had mild conjunctivitis. Four days later, one of the cases was referred to the Communicable Disease Centre, Tan Tock Seng Hospital. The patient’s blood and urine samples tested positive for Zika virus (ZIKV), and the blood was negative for both dengue and chikungunya viruses by polymerase chain reaction (PCR). Since the patient had no travel history within the past month, the MOH announced the first locally transmitted ZIKV infection on 27 August, after laboratory confirmation on a second set of blood and urine specimens. For the purpose of outbreak investigation and public health measures, MOH then directed the collection of clinical specimens from cases working or living in the surrounding areas, with fever (>37.8°C), rash, conjunctivitis and/or joint pain in the preceding two weeks. All clinical information and samples pertaining to the outbreak investigation were collected under the provisions of the Infectious Diseases Act in Singapore [1].

Genome sequencing of outbreak samples

A total of 153 individuals fitting the case definition were tested for ZIKV on 27 and 28 August, of which, only 56 cases were confirmed positive by real-time PCR assay [2]. Respective samples from two cases, ZKA-16–097 and ZKA-16–291, were selected based on their high viral titre (cycle threshold (Ct) values 21.7 and 24.6). Viral RNA was extracted from urine samples using the QIAamp viral RNA kit (Qiagen) and a series of overlapping reverse transcription (RT)-PCR reactions were performed using one-step Superscript III/Hi-Fidelity platinum Taq polymerase (Thermo Fisher). The primers were designed to target conserved regions among ZIKV whole genome sequences of both Asian and African lineages that were available on GenBank [3]. DNA bands of the predicted sizes were purified and sequenced by Sanger sequencing. Raw sequences were aligned and edited using CLC workbench. The consensus sequences for ZKA-16–097 and ZKA-16-291 were submitted to GenBank under the accession numbers KX813683 and KX827309, respectively.

Background on Zika virus

ZIKV is a mosquito-borne single-stranded positive-sense RNA virus belonging to the Flaviviridae family. First isolated in 1947 from a sentinel rhesus monkey in Uganda [4], ZIKV circulated enzoonotically within Africa and equatorial Asia as two distinct lineages: the African and Asian lineage [5]. Prior to 2007, only 14 sporadic human infections, confined to Africa and equatorial Asia, were documented [6]. Since then, three outbreaks of ZIKV belonging to the Asian lineage have ensued: in 2007, on Yap island within the Federated States of Micronesia, in 2013 and 2014, within the French Polynesian islands and most notably, the current large outbreak in the Americas which was first detected in Brazil in 2015. Unlike the Yap island outbreak which was characterised by cases with relatively mild dengue-like symptoms [7], the outbreaks in French Polynesia and Brazil coincided with an unusual rate of cerebral congenital anomalies, including microcephaly [8]. Since 2015, 70 countries and territories reported evidence of ZIKV transmission [9]. The magnitude of spread of the Asian lineage and its disease association prompted the World Health Organization to declare ZIKV as a Public Health Emergency of International Concern in 2016 [10].

Phylogenetic analysis

To gain further information on the ZIKV strains circulating within the current Singapore outbreak, the sequences recovered from the two cases were phylogenetically analysed. Complete ZIKV genomes available on GenBank were downloaded [3], and selected non-redundant representatives were aligned with the sequences from Singapore using Multiple Alignment using Fast Fourier Transform (MAFFT) [11]. The two Singapore sequences are 99.9% identical with only seven nucleotides different over 10,272 bases in the coding region. A maximum likelihood (ML) phylogenetic tree was created in Molecular Evolutionary Genetics Analysis (MEGA) [12] using the Tamura–Nei model with gamma distributed rate differences (5 categories, including invariant) and 1,000 bootstrap step validation. The phylogeny (Figure 1A) clearly shows the separate African and Asian lineages. Within the Asian lineage, the two sequences found in Singapore form a distinct branch, which stems from an ancestral node separating this branch from the tight cluster of strains from French Polynesia (2013), Haiti (2014) and Brazil (2015-2016) with 98% bootstrap support. The preceding shared ancestral node is with a strain circulating in Thailand in 2014 with 100% bootstrap support. Consequently, the viruses detected in Singapore evolutionarily arose in parallel to the large tight cluster of recent strains in South and Central America indicating that the two cases studied here were not infected by viruses imported from the Americas but rather by representatives of the Asian lineage already circulating in south-east Asia. The same clustering can be shown using other tree building methods and parameters (Bayesian: Figure 1B; neighbour joining, uniform rates, HKY model: data not shown).

Figure 1

Phylogenetic trees with Zika virus genomes, 1947−2016


A) Maximum likelihood phylogenetic tree of representatives from the African lineage coloured red and Asian lineage coloured blue or cyan. The tight cluster with strains from French Polynesia and the Americas is in blue and strains from south-east Asia are in cyan. The scale beneath the tree represents the number of substitutions per site. B) Time-resolved Bayesian tree of recent Asian lineage. The scale beneath the tree represents the years of divergence. C) Correlation of sampling date with phylogenetic distance for the recent Asian lineage.

Furthermore, in the tree, the long branch of the two new south-east Asian viruses found in this study suggests undetected evolution for several years. In order to estimate the time of divergence from the last common ancestor with the later cluster from French Polynesia and the Americas, we used a phylogenetic molecular clock approach through the Bayesian Evolutionary Analysis Sampling Trees (BEAST) package [13]. We first confirmed that the measured genetic distance within our set of the recent Asian lineage correlates well with the sample dates (Figure 1C) and created a time-resolved Bayesian tree with parameters previously established for ZIKV (strict clock model, Bayesian skyline, generalised time-reversible (GTR) [14]). The resulting tree (Figure 1B) is consistent with the clustering described in the ML tree and the time of the last common ancestor of the new south-east Asian strains and the French Polynesia/Americas cluster is estimated as 6.2 years ago (95% highest posterior density interval: 4.58–8.15 years) suggesting the observed clades diverged from each other around early 2010. Given this long unsampled evolution, future studies should establish possible reservoirs and circulation in the region.

Conclusions from sequence and structure analysis

This study shows that there are still multiple ZIKV strains in circulation globally in 2016 which raises the question on their antigenic diversity or similarity for vaccine development. Therefore, we systematically mapped the outer surface envelope (E) protein changes from 1947 onwards across the whole Asian lineage (Figure 2).

Figure 2

Mutations on the envelope (E) protein surface of different Zika virus strains


Mutated residues in the Asian strain of Zika virus relative to the reference strain MR766 are indicated in cyan. Amino-acid residue numbers are according to polyprotein numbering. The highly conserved glycosylation site at N444 is indicated in blue. The schematic in the bottom left illustrates how the E proteins are arranged in a set of dimeric, trimeric, and pentameric elements to form a pseudo-icosahedral structure on the viral surface. 3D homology models were created using MODELLER [15] based on PDB:5IZ7 [16] and visualised with the PyMOL Molecular Graphics System (Version 1.8 Schrödinger, LLC.). The percent identity of each depicted strain’s E protein relative to the consensus of the Brazil 2015 sequences ranges from 96.8 to 100% and is indicated in parentheses.

At least within the Asian lineage the antigenic protein E surface is highly conserved and homogeneous with, from 2007 onwards, typical identities from 99.4% to 100% (or 3 to zero mutations) relative to the consensus of recent strains from Brazil. Moreover all of the E proteins in this lineage contain the typical N-glycosylation site (N444 in polyprotein numbering which is N154 in the E protein). The lack of surface mutation drift over the past 50 years (99% identity between Malaysia 1966 and consensus of Brazil 2015) also suggests that immune pressure on the E protein has not been a dominant factor in the Asian lineages’ virus fitness and evolution so far. While the highly similar surface E proteins of the different strains in the Asian lineage should facilitate global vaccine development it is too early to judge if the strain linked to the Singapore outbreak would show any different disease characteristics from the one in French Polynesia and the Americas. Further studies are necessary and underway in Singapore to collect more data.


NPHL and its surveillance programmes are fully funded by Ministry of Health (MOH), Singapore and the other authors are funded by the Agency for Science Technology and Research (A*STAR) Singapore (HEIDI programme) and the Ministry of Education in Singapore (MOE AcRF Tier 3 Grant Number MOE2012-T3-1-008 co-held by Chandra Verma and Peter Bond of the Bioinformatics Institute (A*STAR)). The authors would like to thank Communicable Disease Division, MOH, Singapore, for its roles in field investigation and sample collection.

Conflict of interest

None declared.

Authors’ contributions

TMM, YKN, SPP, LC, RL contributed sample processing, primer preparation, sequencing, data analysis, interpretation and co-wrote the paper. SMS, RGH, JKM, DAH and RTCL contributed phylogenetic analysis, structural modelling and co-wrote the paper.


  1. Ministry of Health Singapore. Infectious Diseases Act. Singapore: Ministry of Health Singapore. [Accessed 21 Sep 2016]. Available from:

  2. Lanciotti RS, Kosoy OL, Laven JJ, Velez JO, Lambert AJ, Johnson AJ,  et al.  Genetic and serologic properties of Zika virus associated with an epidemic, Yap State, Micronesia, 2007. Emerg Infect Dis. 2008;14(8):1232-9. DOI: 10.3201/eid1408.080287 PMID: 18680646

  3. NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information.Nucleic Acids Res. 2016;44(D1):D7-19. DOI: 10.1093/nar/gkv1290 PMID: 26615191

  4. Dick GWA, Kitchen SF, Haddow AJ. Zika virus. I. Isolations and serological specificity.Trans R Soc Trop Med Hyg. 1952;46(5):509-20. DOI: 10.1016/0035-9203(52)90042-4 PMID: 12995440

  5. Haddow AD, Schuh AJ, Yasuda CY, Kasper MR, Heang V, Huy R,  et al.  Genetic characterization of Zika virus strains: geographic expansion of the Asian lineage. PLoS Negl Trop Dis. 2012;6(2):e1477. DOI: 10.1371/journal.pntd.0001477 PMID: 22389730

  6. Musso D, Gubler DJ. Zika Virus.Clin Microbiol Rev. 2016;29(3):487-524. DOI: 10.1128/CMR.00072-15 PMID: 27029595

  7. Duffy MR, Chen T-H, Hancock WT, Powers AM, Kool JL, Lanciotti RS,  et al.  Zika virus outbreak on Yap Island, Federated States of Micronesia. N Engl J Med. 2009;360(24):2536-43. DOI: 10.1056/NEJMoa0805715 PMID: 19516034

  8. European Centre for Disease Prevention and Control (ECDC). Rapid Risk Assessment. Zika virus disease epidemic: potential association with microcephaly and Guillain–Barré syndrome. Second update- 8 February 2016. Stockholm: ECDC; 2016. [Accessed 10 Sep 2016]. Available from:

  9. World Health Organization (WHO). Zika virus and complications. Geneva: WHO. [Accessed 10 Sep 2016]. Available from:

  10. World Health Organization (WHO). WHO Director-General summarizes the outcome of the Emergency Committee regarding clusters of microcephaly and Guillain-Barré syndrome. Geneva: WHO; 2016. Available from:

  11. Katoh K, Standley DM. MAFFT: iterative refinement and additional methods.Methods Mol Biol. 2014;1079:131-46. DOI: 10.1007/978-1-62703-646-7_8 PMID: 24170399

  12. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0.Mol Biol Evol. 2013;30(12):2725-9. DOI: 10.1093/molbev/mst197 PMID: 24132122

  13. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7.Mol Biol Evol. 2012;29(8):1969-73. DOI: 10.1093/molbev/mss075 PMID: 22367748

  14. Faria NR, Azevedo RS, Kraemer MU, Souza R, Cunha MS, Hill SC,  et al.  Zika virus in the Americas: Early epidemiological and genetic findings. Science. 2016;352(6283):345-9. DOI: 10.1126/science.aaf5036 PMID: 27013429

  15. Eswar N, Eramian D, Webb B, Shen MY, Sali A. Protein structure modeling with MODELLER.Methods Mol Biol. 2008;426:145-59. DOI: 10.1007/978-1-60327-058-8_8 PMID: 18542861

  16. Kostyuchenko VA, Lim EXY, Zhang S, Fibriansah G, Ng T-S, Ooi JSG,  et al.  Structure of the thermally stable Zika virus. Nature. 2016;533(7603):425-8.PMID: 27093288

Back to Table of Contents
Previous Download (pdf)

The publisher’s policy on data collection and use of cookies.

Disclaimer: The opinions expressed by authors contributing to Eurosurveillance do not necessarily reflect the opinions of the European Centre for Disease Prevention and Control (ECDC) or the editorial team or the institutions with which the authors are affiliated. Neither ECDC nor any person acting on behalf of ECDC is responsible for the use that might be made of the information in this journal. The information provided on the Eurosurveillance site is designed to support, not replace, the relationship that exists between a patient/site visitor and his/her physician. Our website does not host any form of commercial advertisement. Except where otherwise stated, all manuscripts published after 1 January 2016 will be published under the Creative Commons Attribution (CC BY) licence. You are free to share and adapt the material, but you must give appropriate credit, provide a link to the licence, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.

Eurosurveillance [ISSN] - ©2007-2016. All rights reserved

This website is certified by Health On the Net Foundation. Click to verify. This site complies with the HONcode standard for trustworthy health information:
verify here.