Inter- and intra-lineage genetic diversity of wild-type Zika viruses reveals both common and distinctive nucleotide variants and clusters of genomic diversity

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

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.

African ZIKV lineage fails to sustain infectivity in an in vitro mimetic urban cycle

Zika virus (ZIKV) is an arbovirus maintained in nature in two distinct cycles of transmission: urban and sylvatic. Each cycle includes specific vertebrate and invertebrate hosts, and through alternate infections, a conserved consensus sequence is maintained that might vary depending on the cycle. The current study aimed to investigate the ability of ZIKVAF and ZIKVBR to maintain an infectious cycle by alternating passages in cells mimicking the urban (UC) and semi-sylvatic (SC) cycles. The complete genome of the original inoculum and the last passages for each cycle were sequenced by Sanger. Ten passages were performed, as planned, for ZIKVBR UC, ZIKVAF SC, and ZIKVBR SC. ZIKVBR SC showed significant variation in viral titers along the passages, suggesting that the virus is not well adapted to the non-human primate host. ZIKVAF passage in UC was abrogated in the third passage, showing the inability of the African lineage to sustain cycles in human cells, suggesting a low capacity to establish an urban cycle. Several mutations were found in both strains along the passages, but not occurring at equivalent positions. Further studies are needed to elucidate whether any of these specific mutations affect viral fitness. ZIKV strains behave differently in artificial transmission cycles in vitro: Brazilian ZIKV was able to establish urban and semi-sylvatic cycles in vitro. African ZIKV proved unable to cycle among human and mosquito cells and is compatible only with the semi-sylvatic cycle. The main mutations arose in the NS2A region after artificial transmission cycles for both ZIKV strains but not at equivalent positions.

Tissue-specific expansion of Zika virus isogenic variants drive disease pathogenesis

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).

High-throughput sequencing approaches applied to SARS-CoV-2

High-throughput sequencing is crucial for surveillance and control of viral outbreaks. During the ongoing coronavirus disease 2019 (COVID-19) pandemic, advances in the high-throughput sequencing technology resources have enhanced diagnosis, surveillance, and vaccine discovery. From the onset of the pandemic in December 2019, several genome-sequencing approaches have been developed and supported across the major sequencing platforms such as Illumina, Oxford Nanopore, PacBio, MGI DNBSEQTM and Ion Torrent. Here, we share insights from the sequencing approaches developed for sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) between December 2019 and October 2022.

Unravelling the epidemiological diversity of Zika virus by analyzing key protein variations

The consequences of Zika virus (ZIKV) infections were limited to sporadic mild diseases until almost a decade ago, when epidemic outbreaks took place, with quick spread into the Americas. Simultaneously, novel severe neurological manifestations of ZIKV infections were identified, including congenital microcephaly. However, why the epidemic strains behave differently is not yet completely understood, and many questions remain about the actual significance of genetic variations in the epidemiology and biology of ZIKV. In this study, we analysed a large number of viral sequences to identify genes with different levels of variability and patterns of genomic variations that could be associated with ZIKV diversity. We compared numerous epidemic strains with pre-epidemic strains, using the BWA-mem algorithm, and we also examined specific variations among the epidemic ZIKV strains derived from microcephaly cases. We identified several viral genes with dissimilar mutation rates among the ZIKV strain groups and novel protein variation profiles that might be associated with epidemiological particularities. Finally, we assessed the impact of the detected changes on the structure and stability of the NS1, NS5, and E proteins using the I-TASSER, trRosetta, and RaptorX modelling algorithms, and we found some interesting variations that might help to explain the heterogeneous features of the diverse ZIKA strains. This work contributes to the identification of genetic differences in the ZIKV genome that might have a phenotypic impact, providing a basis for future experimental analysis to elucidate the genetic causes of the recent ZIKV emergency.

A safe replication-defective Zika virus vaccine protects mice from viral infection and vertical transmission

With the explosive emergence of Zika virus (ZIKV) and the consequent devastating fetal malformations in infected expectant women, a safe and effective vaccine is urgently needed. Here, using our established NS1 trans-complementation system, we generated high titer of replication-defective ZIKV with NS1 deletion (ZIKV-ΔNS1) in the BHK-21 cell line stably expressing NS1 (BHK_NS1). NS1 deletion of ZIKV-ΔNS1 was stably maintained as no replicative virus was found in naïve BHK-21 cells after continuous passaging of ZIKV-ΔNS1 in BHK_NS1 cells. The safety of ZIKV-ΔNS1 was demonstrated when a high dose of ZIKV-ΔNS1 (10⁷ IU) was used to infect the highly susceptible type I and type II interferon (IFN) receptor-deficient mice. ZIKV-ΔNS1 could induce antibody responses in both immunocompetent (BALB/c) and immunodeficient mice and a single dose of ZIKV-ΔNS1 vaccine protected the immunodeficient mice from a highly lethal dosage of challenge with WT ZIKV. ZIKV-ΔNS1 immunization also attenuated vertical transmission during pregnancy of type I IFN receptor-deficient IFNAR^-/- mice and protected fetuses from ZIKV infection. Our data reported here not only provide a promising ZIKV vaccine candidate with a satisfied balance between safety and efficacy, but also demonstrate the potential of the NS1 trans-complementation system as a platform for flavivirus vaccine development, especially for highly pathogenic flaviviruses.

An overview of Zika virus genotypes and their infectivity

Zika virus (ZIKV) is an enveloped, single-stranded RNA arbovirus belonging to the genus Flavivirus. It was first isolated from a sentinel monkey in Uganda in 1947. More recently, ZIKV has undergone rapid geographic expansion and has been responsible for outbreaks in Southeast Asia, the Pacific Islands, and America. In this review, we have highlighted the influence of viral genetic variants on ZIKV pathogenesis. Two major ZIKV genotypes (African and Asian) have been identified. The Asian genotype is subdivided into Southwest Asia, Pacific Island, and American strains, and is responsible for most outbreaks. Non-synonymous mutations in ZIKV proteins C, prM, E, NS1, NS2A, NS2B, NS3, and NS4B were found to have a higher prevalence and association with virulent strains of the Asian genotype. Consequently, the Asian genotype appears to have acquired higher cellular permissiveness, tissue persistence, and viral tropism in human neural cells. Therefore, mutations in specific coding regions of the Asian genotype may enhance ZIKV infectivity. Considering that mutations in the genomes of emerging viruses may lead to new virulent variants in humans, there is a potential for the re-emergence of new ZIKV cases in the future.

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.

Neuroinvasiveness of the MR766 strain of Zika virus in IFNAR-/- mice maps to prM residues conserved amongst African genotype viruses

Zika virus (ZIKV) strains are classified into the African and Asian genotypes. The higher virulence of the African MR766 strain, which has been used extensively in ZIKV research, in adult IFNα/β receptor knockout (IFNAR-/-) mice is widely viewed as an artifact associated with mouse adaptation due to at least 146 passages in wild-type suckling mouse brains. To gain insights into the molecular determinants of MR766's virulence, a series of genes from MR766 were swapped with those from the Asian genotype PRVABC59 isolate, which is less virulent in IFNAR-/- mice. MR766 causes 100% lethal infection in IFNAR-/- mice, but when the prM gene of MR766 was replaced with that of PRVABC59, the chimera MR/PR(prM) showed 0% lethal infection. The reduced virulence was associated with reduced neuroinvasiveness, with MR766 brain titers ≈3 logs higher than those of MR/PR(prM) after subcutaneous infection, but was not significantly different in brain titers of MR766 and MR/PR(prM) after intracranial inoculation. MR/PR(prM) also showed reduced transcytosis when compared with MR766 in vitro. The high neuroinvasiveness of MR766 in IFNAR-/- mice could be linked to the 10 amino acids that differ between the prM proteins of MR766 and PRVABC59, with 5 of these changes affecting positive charge and hydrophobicity on the exposed surface of the prM protein. These 10 amino acids are highly conserved amongst African ZIKV isolates, irrespective of suckling mouse passage, arguing that the high virulence of MR766 in adult IFNAR-/- mice is not the result of mouse adaptation.

Systemic inflammation, innate immunity and pathogenesis after Zika virus infection in cynomolgus macaques are modulated by strain-specificity within the Asian lineage

Zika virus (ZIKV) is an emerging arbovirus with recent global expansion. Historically, ZIKV infections with Asian lineages have been associated with mild disease such as rash and fever. However, recent Asian sub-lineages have caused outbreaks in the South Pacific and Latin America with increased prevalence of neurological disorders in infants and adults. Asian sub-lineage differences may partially explain the range of disease severity observed. However, the effect of Asian sub-lineage differences on pathogenesis remains poorly characterized. Current study conducts a head-to-head comparison of three Asian sub-lineages that are representative of the circulating ancestral mild Asian strain (ZIKV-SG), the 2007 epidemic French Polynesian strain (ZIKV-FP), and the 2013 epidemic Brazil strain (ZIKV-Brazil) in adult Cynomolgus macaques. Animals infected intervenously or subcutaneously with either of the three clinical isolates showed sub-lineage-specific differences in viral pathogenesis, early innate immune responses and systemic inflammation. Despite the lack of neurological symptoms in infected animals, the epidemiologically neurotropic ZIKV sub-lineages (ZIKV-Brazil and/or ZIKV-FP) were associated with more sustained viral replication, higher systemic inflammation (i.e. higher levels of TNFα, MCP-1, IL15 and G-CSF) and greater percentage of CD14+ monocytes and dendritic cells in blood. Multidimensional analysis showed clustering of ZIKV-SG away from ZIKV-Brazil and ZIKV-FP, further confirming sub-lineage differences in the measured parameters. These findings highlight greater systemic inflammation and monocyte recruitment as possible risk factors of adult ZIKV disease observed during the 2007 FP and 2013 Brazil epidemics. Future studies should explore the use of anti-inflammatory therapeutics as early treatment to prevent ZIKV-associated disease in adults.

Attenuated dengue viruses are genetically more diverse than their respective wild-type parents

Dengue poses a significant burden of individual health, health systems and the economy in dengue endemic regions. As such, dengue vaccine development has been an active area of research. Previous studies selected attenuated vaccine candidates based on plaque size. However, these candidates led to mixed safety outcome in clinical trials, suggesting it is insufficiently informative as an indicator of dengue virus (DENV) attenuation. In this study, we examined the genome diversity of wild-type DENVs and their attenuated derivatives developed by Mahidol University and tested in phase 1 clinical trials. We found that the attenuated DENVs, in particular the strain under clinical development by Takeda Vaccines, DENV2 PDK53, showed significantly higher genome diversity than its wild-type parent, DENV2 16681. The determinant of genomic diversity was intrinsic to the PDK53 genome as infectious clone of PDK53 showed greater genomic diversity after a single in vitro passage compared to 16681 infectious clone. Similar trends were observed with attenuated DENV1 and DENV4, both of which were shown to be attenuated clinically, but not DENV3 that was not adequately attenuated clinically. Taken together, evidence presented here suggests that genome diversity could be developed into a marker of DENV attenuation.

Zika Virus Pathogenesis: A Battle for Immune Evasion

Zika virus (ZIKV) infection and its associated congenital and other neurological disorders, particularly microcephaly and other fetal developmental abnormalities, constitute a World Health Organization (WHO) Zika Virus Research Agenda within the WHO’s R&D Blueprint for Action to Prevent Epidemics, and continue to be a Public Health Emergency of International Concern (PHEIC) today. ZIKV pathogenicity is initiated by viral infection and propagation across multiple placental and fetal tissue barriers, and is critically strengthened by subverting host immunity. ZIKV immune evasion involves viral non-structural proteins, genomic and non-coding RNA and microRNA (miRNA) to modulate interferon (IFN) signaling and production, interfering with intracellular signal pathways and autophagy, and promoting cellular environment changes together with secretion of cellular components to escape innate and adaptive immunity and further infect privileged immune organs/tissues such as the placenta and eyes. This review includes a description of recent advances in the understanding of the mechanisms underlying ZIKV immune modulation and evasion that strongly condition viral pathogenesis, which would certainly contribute to the development of anti-ZIKV strategies, drugs, and vaccines.

A genetically stable Zika virus vaccine candidate protects mice against virus infection and vertical transmission

Although live attenuated vaccines (LAVs) have been effective in the control of flavivirus infections, to date they have been excluded from Zika virus (ZIKV) vaccine trials due to safety concerns. We have previously reported two ZIKV mutants, each of which has a single substitution in either envelope (E) glycosylation or nonstructural (NS) 4B P36 and displays a modest reduction in mouse neurovirulence and neuroinvasiveness, respectively. Here, we generated a ZIKV mutant, ZE4B-36, which combines mutations in both E glycosylation and NS4B P36. The ZE4B-36 mutant is stable and attenuated in viral replication. Next-generation sequence analysis showed that the attenuating mutations in the E and NS4B proteins are retained during serial cell culture passages. The mutant exhibits a significant reduction in neuroinvasiveness and neurovirulence and low infectivity in mosquitoes. It induces robust ZIKV-specific memory B cell, antibody, and T cell-mediated immune responses in type I interferon receptor (IFNR) deficient mice. ZIKV-specific T cell immunity remains strong months post-vaccination in wild-type C57BL/6 (B6) mice. Vaccination with ZE4B-36 protects mice from ZIKV-induced diseases and vertical transmission. Our results suggest that combination mutations in E glycosylation and NS4B P36 contribute to a candidate LAV with significantly increased safety but retain strong immunogenicity for prevention and control of ZIKV infection.

Mapping the evolutionary landscape of Zika virus infection in immunocompromised mice

The fundamental basis of how arboviruses evolve in nature and what regulates the adaptive process remain unclear. To address this problem, we established a Zika virus (ZIKV) vector-borne transmission system in immunocompromised mice to study the evolutionary characteristics of ZIKV infection. Using this system, we defined factors that influence the evolutionary landscape of ZIKV infection and show that transmission route and specific organ microenvironments impact viral diversity and defective viral genome production. In addition, we identified in mice the emergence of ZIKV mutants previously seen in natural infections, including variants present in currently circulating Asian and American strains, as well as mutations unique to the mouse infections. With these studies, we have established an insect-to-mouse transmission model to study ZIKV evolution in vivo. We also defined how organ microenvironments and infection route impact the ZIKV evolutionary landscape, providing a deeper understanding of the factors that regulate arbovirus evolution and emergence.

Using Next Generation Sequencing to Study the Genetic Diversity of Candidate Live Attenuated Zika Vaccines

Zika virus (ZIKV) is a mosquito-transmitted positive-sense RNA virus in the family Flaviviridae. Candidate live-attenuated vaccine (LAV) viruses with engineered deletions in the 3’ untranslated region (UTR) provide immunity and protection in animal models of ZIKV infection, and phenotypic studies show that LAVs retain protective abilities following in vitro passage. The present study investigated the genetic diversity of wild-type (WT) parent ZIKV and its candidate LAVs using next generation sequencing analysis of five sequential in vitro passages. The results show that genomic entropy of WT ZIKV steadily increases during in vitro passage, whereas that of LAVs also increased by passage number five but was variable throughout passaging. Additionally, clusters of single nucleotide variants (SNVs) were found to be present in the pre-membrane/membrane (prM), envelope (E), nonstructural protein NS1 (NS1), and other nonstructural protein genes, depending on the specific deletion, whereas in the parent WT ZIKV, they are more abundant in prM and NS1. Ultimately, both the parental WT and LAV derivatives increase in genetic diversity, with evidence of adaptation following passage.

Attenuation of Zika Virus by Passage in Human HeLa Cells

Zika virus (ZIKV) is a mosquito-borne Flavivirus. Previous studies have shown that mosquito-transmitted flaviviruses, including yellow fever, Japanese encephalitis, and West Nile viruses, could be attenuated by serial passaging in human HeLa cells. Therefore, it was hypothesized that wild-type ZIKV would also be attenuated after HeLa cell passaging. A human isolate from the recent ZIKV epidemic was subjected to serial HeLa cell passaging, resulting in attenuated in vitro replication in both Vero and A549 cells. Additionally, infection of AG129 mice with 10 plaque forming units (pfu) of wild-type ZIKV led to viremia and mortality at 12 days, whereas infection with 10³ pfu of HeLa-passage 6 (P6) ZIKV led to lower viremia, significant delay in mortality (median survival: 23 days), and increased cytokine and chemokine responses. Genomic sequencing of HeLa-passaged virus identified two amino acid substitutions as early as HeLa-P3: pre-membrane E87K and nonstructural protein 1 R103K. Furthermore, both substitutions were present in virus harvested from HeLa-P6-infected animal tissue. Together, these data show that, similarly to other mosquito-borne flaviviruses, ZIKV is attenuated following passaging in HeLa cells. This strategy can be used to improve understanding of substitutions that contribute to attenuation of ZIKV and be applied to vaccine development across multiple platforms.