An Evolutionary Insight into Zika Virus Strains Isolated in the Latin American Region

Epidemiology, genetic diversity, and evolutionary dynamics of Tembusu virus

Tembusu virus (TMUV) is an emerging pathogenic flavivirus associated with acute egg-drop and fatal encephalitis in domestic waterfowl. Since its initial identification in mosquitoes in 1955, TMUV has been confirmed to infect ducks, pigeons, sparrows, geese, and chickens, posing a significant threat to the poultry industry. Here, we sequenced two DTMUV strains isolated in 2019 and systematically investigated the possible origin, genetic relationships, evolutionary dynamics, and transmission patterns of TMUV based on complete virus genome sequences in the public database. We found that TMUV can be divided into four major clusters: TMUV, cluster 1, cluster 2, and cluster 3. Interestingly, we found that cluster 2.2 (within cluster 2) is the most commonly involved in interspecies transmission events, and subcluster 2.1.2 (within cluster 2.1) is currently the most prevalent cluster circulating in Asia. Notably, we also identified three positively selected sites in the E and NS1 proteins, which may be involved in virus replication, immune evasion, and host adaptation. Finally, phylogeographic analysis revealed that cluster dispersal originated in Southeast Asia and that short-distance transmission events have occurred frequently. Altogether, these data provide novel insights into the evolution and dispersal of TMUV, facilitating the development of rapid diagnostics, vaccines, and therapeutics against TMUV infection.

Analyses of Early ZIKV Genomes Are Consistent with Viral Spread from Northeast Brazil to the Americas

The Americas, particularly Brazil, were greatly impacted by the widespread Zika virus (ZIKV) outbreak in 2015 and 2016. Efforts were made to implement genomic surveillance of ZIKV as part of the public health responses. The accuracy of spatiotemporal reconstructions of the epidemic spread relies on the unbiased sampling of the transmission process. In the early stages of the outbreak, we recruited patients exhibiting clinical symptoms of arbovirus-like infection from Salvador and Campo Formoso, Bahia, in Northeast Brazil. Between May 2015 and June 2016, we identified 21 cases of acute ZIKV infection and subsequently recovered 14 near full-length sequences using the amplicon tiling multiplex approach with nanopore sequencing. We performed a time-calibrated discrete phylogeographic analysis to trace the spread and migration history of the ZIKV. Our phylogenetic analysis supports a consistent relationship between ZIKV migration from Northeast to Southeast Brazil and its subsequent dissemination beyond Brazil. Additionally, our analysis provides insights into the migration of ZIKV from Brazil to Haiti and the role Brazil played in the spread of ZIKV to other countries, such as Singapore, the USA, and the Dominican Republic. The data generated by this study enhances our understanding of ZIKV dynamics and supports the existing knowledge, which can aid in future surveillance efforts against the virus.

OUP accepted manuscript

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized ‘African’ and ‘Asian’ genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

IP-10 and CXCR3 signaling inhibit Zika virus replication in human prostate cells

Our previous studies have shown that Zika virus (ZIKV) replicates in human prostate cells, suggesting that the prostate may serve as a long-term reservoir for virus transmission. Here, we demonstrated that the innate immune responses generated to three distinct ZIKV strains (all isolated from human serum) were significantly different and dependent on their passage history (in mosquito, monkey, or human cells). In addition, some of these phenotypic differences were reduced by a single additional cell culture passage, suggesting that viruses that have been passaged more than 3 times from the patient sample will no longer reflect natural phenotypes. Two of the ZIKV strains analyzed induced high levels of the IP-10 chemokine and IFNγ in human prostate epithelial and stromal mesenchymal stem cells. To further understand the importance of these innate responses on ZIKV replication, we measured the effects of IP-10 and its downstream receptor, CXCR3, on RNA and virus production in prostate cells. Treatment with IP-10, CXCR3 agonist, or CXCR3 antagonist significantly altered ZIKV viral gene expression, depending on their passage in cells of relevant hosts (mosquito or human). We detected differences in gene expression of two primary CXCR3 isoforms (CXCR3-A and CXCR3-B) on the two cell types, possibly explaining differences in viral output. Lastly, we examined the effects of IP-10, agonist, or antagonist on cell death and proliferation under physiologically relevant infection rates, and detected no significant differences. Although we did not measure protein expression directly, our results indicate that CXCR3 signaling may be a target for therapeutics, to ultimately stop sexual transmission of this virus.

Codelivery of improved immune complex and virus-like particle vaccines containing Zika virus envelope domain III synergistically enhances immunogenicity

Zika virus (ZIKV) reemergence poses a significant health threat especially due to its risks to fetal development, necessitating safe and effective vaccines that can protect pregnant women. Zika envelope domain III (ZE3) has been identified as a safe and effective vaccine candidate, however it is poorly immunogenic. We previously showed that plant-made recombinant immune complex (RIC) vaccines are a robust platform to improve the immunogenicity of weak antigens. In this study, we altered the antigen fusion site on the RIC platform to accommodate N-terminal fusion to the IgG heavy chain (N-RIC), and thus a wider range of antigens, with a resulting 40% improvement in RIC expression over the normal C-terminal fusion (C-RIC). Both types of RICs containing ZE3 were efficiently assembled in plants and purified to >95% homogeneity with a simple one-step purification. Both ZE3 RICs strongly bound complement receptor C1q and elicited strong ZE3-specific antibody titers that correlated with ZIKV neutralization. When either N-RIC or C-RIC was codelivered with plant-produced hepatitis B core (HBc) virus-like particles (VLP) displaying ZE3, the combination elicited 5-fold greater antibody titers (>1,000,000) and more strongly neutralized ZIKV than either RICs or VLPs alone, after only two doses without adjuvant. These findings demonstrate that antigens that require a free N-terminus for optimal antigen display can now be used with the RIC system, and that plant-made RICs and VLPs are highly effective vaccines targeting ZE3. Thus, the RIC platform can be more generally applied to a wider variety of antigens.

The Rescue and Characterization of Recombinant, Microcephaly-Associated Zika Viruses as Single-Round Infectious Particles

Zika virus (ZIKV) is transmitted by Aedes mosquitoes and exhibits genetic variation with African and Asian lineages. ZIKV Natal RGN strain, an Asian-lineage virus, has been identified in brain tissues from fetal autopsy cases with microcephaly and is suggested to be a neurotropic variant. However, ZIKV Natal RGN strain has not been isolated; its biological features are not yet illustrated. This study rescued and characterized recombinant, single-round infectious particles (SRIPs) of the ZIKV Natal RGN strain using reverse genetic and synthetic biology techniques. The DNA-launched replicon of ZIKV Natal RGN was constructed and contains the EGFP reporter, lacks prM-E genes, and replicates under CMV promoter control. The peak in the ZIKV Natal RGN SRIP titer reached 6.25 × 10⁶ TCID50/mL in the supernatant of prM-E-expressing packaging cells 72 h post-transfection with a ZIKV Natal RGN replicon. The infectivity of ZIKV Natal RGN SRIPs has been demonstrated to correlate with the green florescence intensity of the EGFP reporter, the SRIP-induced cytopathic effect, and ZIKV’s non-structural protein expression. Moreover, ZIKV Natal RGN SRIPs effectively self-replicated in rhabdomyosarcoma/muscle, glioblastoma/astrocytoma, and retinal pigmented epithelial cells, displaying unique cell susceptibility with differential attachment activity. Therefore, the recombinant ZIKV Natal RGN strain was rescued as SRIPs that could be used to elucidate the biological features of a neurotropic strain regarding cell tropism and pathogenic components, apply for antiviral agent screening, and develop vaccine candidates.

Mutational landscape of Zika virus strains worldwide and its structural impact on proteins

Zika virus (ZIKV) has spread globally and has been linked to the onset of microcephaly and other brain abnormalities. The ZIKV genome consists of an ~10.7 kb positive-stranded RNA molecule that encodes three structural (C, prM and E) and seven nonstructural (5'-NS1-NS2A-NS2B-NS3- NS4A/2K-NS4B-NS5-3') proteins. In this work, we looked for genetic variants in 485 ZIKV complete genomes from GenBank (NCBI) and performed a computational systematic approach using MAESTROweb server to assess the impact of nonsynonymous mutations in ZIKV proteins (C, M, E, NS1, NS2A, NS2B-NS3 protease, NS3 helicase and NS5). Then, we merged the data and correlated it with the phenotypic reports of ZIKV circulating strains. The sensitivity profile of the proteins showed 96 mutational hotspots. We found 22 relevant mutations in proteins C (I80T), NS2A (I34M/T/V, I45V, I80T/V, L113F, A117V, I118V, L128P, V143A, T151A, M199I/V, R207K and L208I) and NS3 helicase (D436G, Y498H, R525K, Q528R and R583K) of the circulating strains. Our analysis exploited the impact of nonsynonymous mutations on ZIKV proteins, their structural and functional insights. The results presented here could advance our current understanding on ZIKV proteins functions and pathogenesis.