Background The Asian lineage Zika virus (ZIKV) emerged as a public health emergency in 2016 causing severe neurological pathologies with no apparent historical correlate to the mild, disease-causing innocuous member of the mosquito-borne flavivirus genus that was discovered in Africa in 1947. Replication error rate of RNA viruses combined with viral protein/RNA structural plasticity can lead to evolution of virus-induced pathogenicity that is critical to identify and validate.

Methods Infection studies in cells and A129 interferon alpha/beta receptor deficient mice with ZIKV French Polynesian H/PF/2013 clinical isolate, plaque-purified isogenic clone derivatives as well as infectious cDNA clone derived wild-type and site-specific mutant viruses, were employed together with Next-Generation Sequencing (NGS) to pin-point the contributions of specific viral variants in neurovirulence recapitulated in our ZIKV mouse model.

Findings NGS analysis of the low-passage inoculum virus as well as mouse serum, brain and testis derived virus, revealed specific enrichment in the mouse brain that were not found in the other tissues. Specifically, non-structural (NS) protein 2A variant at position 117 along with changes in NS1 and NS4B were uniquely associated with the mouse brain isolate. Mutational analysis of these variants in cDNA infectious clones identified the NS2A A117V as the lethal pathogenic determinant with potential epistatic contribution of NS1 and NS4B variants in ZIKV brain penetrance.

Interpretation Our findings confirm that viral subpopulations drive ZIKV neuropathogenicity and identify specific sequence variants that expand in the mouse brain that associates with this phenotype which can serve as predictors of severe epidemics.

Funding Duke-NUS Khoo Post-doctoral Fellowship Award 2020 (KWKC) and National Medical Research Council of Singapore grants MOH-000524 (OFIRG) (SW) and MOH-OFIRG20nov-0002 (SGV).