Comparative genomic analysis of pre-epidemic and epidemic Zika virus strains for virological factors potentially associated with the rapidly expanding epidemic

Abstract

Less than 20 sporadic cases of human Zika virus (ZIKV) infection were reported in Africa and Asia before 2007, but large outbreaks involving up to 73% of the populations on the Pacific islands have started since 2007, and spread to the Americas in 2014. Moreover, the clinical manifestation of ZIKV infection has apparently changed, as evident by increasing reports of neurological complications, such as Guillain-Barré syndrome in adults and congenital anomalies in neonates. We comprehensively compared the genome sequences of pre-epidemic and epidemic ZIKV strains with complete genome or complete polyprotein sequences available in GenBank. Besides the reported phylogenetic clustering of the epidemic strains with the Asian lineage, we found that the topology of phylogenetic tree of all coding regions is the same except that of the non-structural 2B (NS2B) coding region. This finding was confirmed by bootscan analysis and multiple sequence alignment, which suggested the presence of a fragment of genetic recombination at NS2B with that of Spondweni virus. Moreover, the representative epidemic strain possesses one large bulge of nine bases instead of an external loop on the first stem-loop structure at the 3'-untranslated region just distal to the stop codon of the NS5 in the 1947 pre-epidemic prototype strain. Fifteen amino acid substitutions are found in the epidemic strains when compared with the pre-epidemic strains. As mutations in other flaviviruses can be associated with changes in virulence, replication efficiency, antigenic epitopes and host tropism, further studies would be important to ascertain the biological significance of these genomic changes.