Molecular characterization and geographical distribution of Zika virus worldwide from 1947 to 2022

OBJECTIVES

We conducted molecular characterization, demonstrated the geographical distribution of Zika virus (ZIKV) circulating worldwide from 1947 to 2022 and explored the potential genetic recombination site in the Thailand ZIKV genomes.

METHODS

We constructed phylogenetic trees based on ZIKV coding sequences (CDS) and determined the geographical distribution of the representative viruses by genetic relationship and timeline. We determined genetic recombination among ZIKV and between ZIKV and other flaviviruses using similarity plot and bootscan analyzes, together with the phylogeny encompassing the CDS and eight subgenomic regions.

RESULTS

The phylogenetic trees comprising 717 CDS showed two distinct African and Asian lineages. ZIKV in the African lineage formed two sublineages, and ZIKV in the Asian lineage diversified into the Asian and American sublineages. The 1966 Malaysian isolate was designated the prototype of the Asian sublineage and formed a node of only one member, while the newer viruses formed a distinct node. We detected no genetic recombination in the Thailand ZIKV.

CONCLUSION

Five Thailand isolates discovered in 2006 were the second oldest ZIKV after the Malaysian prototype. Our result suggested two independent routes of ZIKV spread from Southeast Asia to Micronesia in 2007 and French Polynesia in 2013 before further spreading to South American countries.

Zika virus isolation, propagation, and quantification using multiple methods

Zika virus (ZIKV) was isolated from the archival urine, serum, and autopsy specimens by intrathoracic inoculation of Toxorhynchitis splendens and followed by three blind sub-passaging in C6/36 mosquito cells. The virus isolates were identified using an immunofluorescence assay and real-time reverse transcription-polymerase chain reaction (real-time RT-PCR). This study analyzed 11 ZIKV isolates. One isolate (0.6%) was obtained from 171 urine samples, eight (8.7%) from 92 serum samples and two from tissues of an abortive fetus. After propagation in C6/36 cells, ZIKV was titrated by plaque and focus forming unit (FFU) assays in Vero cell monolayers, and viral genomes were determined via real-time and digital RT-PCR. Plaque and FFU assay quantitations were comparable, with the amount of infectious viruses averaging 106-107 PFU or FFU/ml. Real-time RT-PCR semi-quantified the viral genome numbers, with Ct values varying from 12 to 14. Digital RT-PCR, which precisely determines the numbers of the viral genomes, consistently averaged 10-100 times higher than the number of infectious units. There was good correlation between the results of these titration methods. Therefore, the selection of a method should be based on the objectives of each research studies.