Genetic diversity and phylogenetic analyses of Asian lineage Zika virus whole genome sequences derived from Culex quinquefasciatus mosquitoes and urine of patients during the 2020 epidemic in Thailand

Abstract

Zika virus (ZIKV), a mosquito-borne flavivirus, has been continually emerging and re-emerging since 2010, with sporadic cases reported annually in Thailand, peaking at over 1000 confirmed positive cases in 2016. Leveraging high-throughput sequencing technologies, specifically whole genome sequencing (WGS), has facilitated rapid pathogen genome sequencing. In this study, we used multiplex amplicon sequencing on the Illumina Miseq instrument to describe ZIKV WGS. Six ZIKV WGS were derived from three samples of field-caught Culex quinquefasciatus mosquitoes (two males and one female) and three urine samples collected from patients in three different provinces of Thailand. Additionally, successful isolation of a ZIKV isolate occurred from a female Cx. quinquefasciatus. The WGS analysis revealed a correlation between the 2020 outbreak and the acquisition of five amino acid changes in the Asian lineage ZIKV strains from Thailand (2006), Cambodia (2010 and 2019), and the Philippines (2012). These changes, including C-T106A, prM-V1A, E-V473M, NS1-A188V, and NS5-M872V, were identified in all seven WGS, previously linked to significantly higher mortality rates. Furthermore, phylogenetic analysis indicated that the seven ZIKV sequences belonged to the Asian lineage. Notably, the genomic region of the E gene showed the highest nucleotide diversity (0.7–1.3%). This data holds significance in informing the development of molecular tools that enhance our understanding of virus patterns and evolution. Moreover, it may identify targets for improved methods to prevent and control future ZIKV outbreaks.