Zika Virus Genome Research Articles

Genomic variability in Zika virus in GBS cases in Colombia.

Major clusters of Guillain-Barré Syndrome (GBS) emerged during the Zika virus (ZIKV) outbreaks in the South Pacific and the Americas from 2014 to 2016. The factors contributing to GBS susceptibility in ZIKV infection remain unclear, although considerations of viral variation, patient susceptibility, environmental influences, and other potential factors have been hypothesized. Studying the role of viral genetic factors has been challenging due to the low viral load and rapid viral clearance from the blood after the onset of Zika symptoms. The prolonged excretion of ZIKV in urine by the time of GBS onset, when the virus is no longer present in the blood, provides an opportunity to unravel whether specific ZIKV mutations are related to the development of GBS in certain individuals. This study aimed to investigate the association between specific ZIKV genotypes and the development of GBS, taking advantage of a unique collection of ZIKV-positive urine samples obtained from GBS cases and controls during the 2016 ZIKV outbreak in Colombia. Utilizing Oxford-Nanopore technology, we conducted complete genome sequencing of ZIKV in biological samples from 15 patients with GBS associated with ZIKV and 17 with ZIKV infection without neurological complications. ZIKV genotypes in Colombia exhibited distribution across three clades (average bootstrap of 90.9±14.9%), with two clades dominating the landscape. A comparative analysis of ZIKV genomes from GBS and non-neurological complications, alongside 1368 previously reported genomes, revealed no significant distinctions between the two groups. Both genotypes were similarly distributed among observed clades in Colombia. Furthermore, no variations were identified in the amino acid composition of the viral genome between the two groups. Our findings suggest that GBS in ZIKV infection is perhaps associated with patient susceptibility and/or other para- or post-infectious immune-mediated mechanisms rather than with specific ZIKV genome variations.

Analysis of Preferred Codon Usage in the ZIKA Virus Genome and Their Implications for Genome Evolution and Vaccine Design

Analysis of Preferred Codon Usage in the ZIKA Virus Genome and Their Implications for Genome Evolution and Vaccine Design

An accumulated mutation gained in mosquito cells enhances Zika virus virulence and fitness in mice.

Zika virus (ZIKV) is an important arbovirus with a global impact. Experimental evolution by serial passaging of ZIKV in susceptible cells has led to the identification of a panel of critical amino acid substitutions with specific functions. Herein, we identified a mosquito cell-derived substitution, H401Y, in the ZIKV E protein via experimental evolution. The H401Y substitution significantly enhanced viral virulence and fitness in mammal cells and mice. Notably, the H401Y substitution has been detected in recent mosquito and human isolates from regions spanning Asia to the Americas. Our work elucidates unrecognized virulence determinant in the ZIKV genome that warrants urgent attention. Moreover, the findings underscore the critical need for extensive molecular surveillance and rigorous clinical observation to establish the potential impact in natural circulation. These endeavors are crucial for unraveling the potential of mutation to act as a catalyst for future epidemics, thereby preempting the public health challenges it may pose.

Structure of an RNA G-quadruplex from the West Nile virus genome

Potential G-quadruplex sites have been identified in the genomes of DNA and RNA viruses and proposed as regulatory elements. The genus Orthoflavivirus contains arthropod-transmitted, positive-sense, single-stranded RNA viruses that cause significant human disease globally. Computational studies have identified multiple potential G-quadruplex sites that are conserved across members of this genus. Subsequent biophysical studies established that some G-quadruplexes predicted in Zika and tickborne encephalitis virus genomes can form and known quadruplex binders reduced viral yields from cells infected with these viruses. The susceptibility of RNA to degradation and the variability of loop regions have made structure determination challenging. Despite these difficulties, we report a high-resolution structure of the NS5-B quadruplex from the West Nile virus genome. Analysis reveals two stacked tetrads that are further stabilized by a stacked triad and transient noncanonical base pairing. This structure expands the landscape of solved RNA quadruplex structures and demonstrates the diversity and complexity of biological quadruplexes. We anticipate that the availability of this structure will assist in solving further viral RNA quadruplexes and provides a model for a conserved antiviral target in Orthoflavivirus genomes.

PfAgo-Based Zika Virus Detection.

As a mosquito-borne flavivirus, Zika virus (ZIKV) has been identified as a global health threat. The virus has been linked to severe congenital disabilities, including microcephaly and other congenital malformations, resulting in fatal intrauterine death. Therefore, developing sensitive and specific methods for the early detection and accurate diagnosis of the ZIKV is essential for controlling its spread and mitigating its impact on public health. Herein, we set up a novel nucleic acid detection system based on Pyrococcus furiosus Argonaute (PfAgo)-mediated nucleic acid detection, targeting the non-structural protein 5 (NS5) region of the ZIKV genome (abbreviated ZIKV-PAND). Without preamplification with the polymerase chain reaction (PCR), the minimum detection concentration (MDC) of ZIKV-PAND was about 10 nM. When introducing an amplification step, the MDC can be dramatically decreased to the aM level (8.3 aM), which is comparable to qRT-PCR assay (1.6 aM). In addition, the diagnostic findings from the analysis of simulated clinical samples or Zika virus samples using ZIKV-PAND show a complete agreement of 100% with qRT-PCR assays. This correlation can aid in the implementation of molecular testing for clinical diagnoses and the investigation of ZIKV infection on an epidemiological scale.

Evidence of Zika Virus Reinfection by Genome Diversity and Antibody Response Analysis, Brazil.

We generated 238 Zika virus (ZIKV) genomes from 135 persons in Brazil who had samples collected over 1 year to evaluate virus persistence. Phylogenetic inference clustered the genomes together with previously reported ZIKV strains from northern Brazil, showing that ZIKV has been remained relatively stable over time. Temporal phylogenetic analysis revealed limited within-host diversity among most ZIKV-persistent infected associated samples. However, we detected unusual virus temporal diversity from >5 persons, uncovering the existence of divergent genomes within the same patient. All those patients showed an increase in neutralizing antibody levels, followed by a decline at the convalescent phase of ZIKV infection. Of interest, in 3 of those patients, titers of neutralizing antibodies increased again after 6 months of ZIKV infection, concomitantly with real-time reverse transcription PCR re-positivity, supporting ZIKV reinfection events. Altogether, our findings provide evidence for the existence of ZIKV reinfection events.

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

ObjectivesWe 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. MethodsWe 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. ResultsThe 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. ConclusionFive 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.

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.

Repurposing of Doramectin as a New Anti-Zika Virus Agent.

Zika virus (ZIKV), belonging to the Flavivirus family and mainly transmitted by mosquitoes, causes a variety of adverse outcomes, including Guillain-Barré syndrome, microcephaly, and meningoencephalitis. However, there are no approved vaccines or drugs available for ZIKV. The discovery and research on drugs for ZIKV are still essential. In this study, we identified doramectin, an approved veterinary antiparasitic drug, as a novel anti-ZIKV agent (EC50 value from 0.85 μM to 3.00 μM) with low cytotoxicity (CC50 > 50 μM) in multiple cellular models. The expression of ZIKV proteins also decreased significantly under the treatment of doramectin. Further study showed that doramectin directly interacted with the key enzyme for ZIKV genome replication, RNA-dependent RNA polymerase (RdRp), with a stronger affinity (Kd = 16.9 μM), which may be related to the effect on ZIKV replication. These results suggested that doramectin might serve as a promising drug candidate for anti-ZIKV.

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.

Correction for Li et al., "Different Degrees of 5'-to-3' DAR Interactions Modulate Zika Virus Genome Cyclization and Host-Specific Replication".

Correction for Li et al., "Different Degrees of 5'-to-3' DAR Interactions Modulate Zika Virus Genome Cyclization and Host-Specific Replication".

Zika virus RNA structure controls its unique neurotropism by bipartite binding to Musashi-1

Human RNA binding protein Musashi-1 (MSI1) plays a critical role in neural progenitor cells (NPCs) by binding to various host RNA transcripts. The canonical MSI1 binding site (MBS), A/GU(1-3)AG single-strand motif, is present in many RNA virus genomes, but only Zika virus (ZIKV) genome has been demonstrated to bind MSI1. Herein, we identified the AUAG motif and the AGAA tetraloop in the Xrn1-resistant RNA 2 (xrRNA2) as the canonical and non-canonical MBS, respectively, and both are crucial for ZIKV neurotropism. More importantly, the unique AGNN-type tetraloop is evolutionally conserved, and distinguishes ZIKV from other known viruses with putative MBSs. Integrated structural analysis showed that MSI1 binds to the AUAG motif and AGAA tetraloop of ZIKV in a bipartite fashion. Thus, our results not only identified an unusual viral RNA structure responsible for MSI recognition, but also revealed a role for the highly structured xrRNA in controlling viral neurotropism.

Comparative mutational analysis of the Zika virus genome from different geographical locations and its effect on the efficacy of Zika virus-specific neutralizing antibodies.

The Zika virus (ZIKV), which originated in Africa, has become a significant global health threat. It is an RNA virus that continues to mutate and accumulate multiple mutations in its genome. These genetic changes can impact the virus's ability to infect, cause disease, spread, evade the immune system, and drug resistance. In this study genome-wide analysis of 175 ZIKV isolates deposited at the National Center for Biotechnology Information (NCBI), was carried out. The comprehensive mutational analysis of these isolates was carried out by DNASTAR and Clustal W software, which revealed 257 different substitutions at the proteome level in different proteins when compared to the reference sequence (KX369547.1). The substitutions were capsid (17/257), preM (17/257), envelope (44/257), NS1 (34/257), NS2A (30/257), NS2B (11/257), NS3 (37/257), NS4A (6/257), 2K (1/257), NS4B (15/257), and NS5 (56/257). Based on the coexisting mutational analysis, the MN025403.1 isolate from Guinea was identified as having 111 substitutions in proteins and 6 deletions. The effect of coexisting/reoccurring mutations on the structural stability of each protein was also determined by I-mutant and MUpro online servers. Furthermore, molecular docking and simulation results showed that the coexisting mutations (I317V and E393D) in Domain III (DIII) of the envelope protein enhanced the bonding network with ZIKV-specific neutralizing antibodies. This study, therefore, highlighted the rapid accumulation of different substitutions in various ZIKV proteins circulating in different geographical regions of the world. Surveillance of such mutations in the respective proteins will be helpful in the development of effective ZIKV vaccines and neutralizing antibody engineering.

Therapeutics for flaviviral infections

Flaviviruses are vector-borne pathogens capable of causing devastating human diseases. The re-emergence of Zika in 2016 notoriously led to a widescale epidemic in the Americas. New daunting evidence suggests that a single mutation in Zika virus genome may increase transmission and pathogenesis, further highlighting the need to be prepared for flavivirus outbreaks. Dengue, in particular infects about 400 million people each year, leading to reoccurring local outbreaks. Public health efforts to mitigate flavivirus transmission is largely dependent on vector control strategies, as only a limited number of flavivirus vaccines have been developed thus far. There are currently no commercially available antivirals for flaviviruses, leaving supportive care as the primary treatment option. In this review, we will briefly paint a broad picture of the flavivirus landscape in terms of therapeutics, with particular focus on viral targets, promising novel compounds entering the drug discovery pipeline, as well as model systems for evaluating drug efficacy.

Zika Virus Infection and Development of Drug Therapeutics

Zika virus (ZIKV) is an emerging flavivirus that is associated with neurological complications, such as neuroinflammatory Guillain Barré Syndrome in adults and microcephaly in newborns, and remains a potentially significant and international public health concern. The World Health Organization is urging the development of novel antiviral therapeutic strategies against ZIKV, as there are no clinically approved vaccines or drugs against this virus. Given the public health crisis that is related to ZIKV cases in the last decade, efficient strategies should be identified rapidly to combat or treat ZIKV infection. Several promising strategies have been reported through drug repurposing studies, de novo design, and the high-throughput screening of compound libraries in only a few years. This review summarizes the genome and structure of ZIKV, viral life cycle, transmission cycle, clinical manifestations, cellular and animal models, and antiviral drug developments, with the goal of increasing our understanding of ZIKV and ultimately defeating it.

CpG content in the Zika virus genome affects infection phenotypes in the adult brain and fetal lymph nodes.

Increasing the number of CpG dinucleotides in RNA viral genomes, while preserving the original amino acid composition, leads to impaired infection which does not cause disease. Beneficially, impaired infection evokes antiviral host immune responses providing a cutting-edge vaccine approach. For example, we previously showed that CpG-enriched Zika virus variants cause attenuated infection phenotypes and protect against lethal challenge in mice. While CpG recoding is an emerging and promising vaccine approach, little is known about infection phenotypes caused by recoded viruses in vivo, particularly in non-rodent species. Here, we used well-established mouse and porcine models to study infection phenotypes of the CpG-enriched neurotropic and congenital virus—Zika virus, directly in the target tissues—the brain and placenta. Specifically, we used the uttermost challenge and directly injected mice intracerebrally to compare infection phenotypes caused by wild-type and two CpG-recoded Zika variants and model the scenario where vaccine strains breach the blood-brain barrier. Also, we directly injected porcine fetuses to compare in utero infection phenotypes and model the scenario where recoded vaccine strains breach the placental barrier. While overall infection kinetics were comparable between wild-type and recoded virus variants, we found convergent phenotypical differences characterized by reduced pathology in the mouse brain and reduced replication of CpG-enriched variants in fetal lymph nodes. Next, using next-generation sequencing for the whole virus genome, we compared the stability of de novo introduced CpG dinucleotides during prolonged virus infection in the brain and placenta. Most de novo introduced CpG dinucleotides were preserved in sequences of recoded Zika viruses showing the stability of vaccine variants. Altogether, our study emphasized further directions to fine-tune the CpG recoding vaccine approach for better safety and can inform future immunization strategies.

Characterization of m6 A modifications in the contemporary Zika virus genome and host cellular transcripts.

Zika virus (ZIKV) suddenly evolved from a neglected arthropod-borne flavivirus into a pandemic pathogen during 2015-2016. A panel of amino acid mutations has been shown to be responsible for the enhanced neurovirulence and transmissibility of ZIKV. Recent studies have demonstrated that ZIKV genomic RNA is modified by host N6-methyladenosine (m6 A) machinery during viral replication in host cells, and the m6 A profiles vary among different isolates and different host cells. In the present study, using a contemporary Asian ZIKV strain isolated in 2019 (SZ1901) as a model, we profiled m6 A modifications on both the viral genome RNA and cellular transcripts from the ZIKV-infected human hepatocarcinoma cell line Huh7. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) identified a unique m6 A map in the genome of ZIKV strain SZ1901 that is different from all previous isolates. Meanwhile, ZIKV infection induced m6 A upregulation in the CDS regions but downregulation in the 3' untranslated region of host RNA transcripts. The m6 A peak intensity in the majority of host genes was downregulated, including ISG-related genes. Overall, our study describes unique viral and host m6 A profiles in contemporary ZIKV-infected Huh7 cells, highlighting the complexity and importance of m6 A modification during viral infection.

Intrinsic features of Zika Virus non-structural proteins NS2A and NS4A in the regulation of viral replication.

Zika virus (ZIKV) is a mosquito-borne flavivirus and can cause neurodevelopmental disorders in fetus. As a neurotropic virus, ZIKV persistently infects neural tissues during pregnancy but the viral pathogenesis remains largely unknown. ZIKV has a positive-sense and single-stranded RNA genome, which encodes 7 non-structural (NS) proteins, participating in viral replication and dysregulation of host immunity. Like those in many other viruses, NS proteins are considered to be products evolutionarily beneficiary to viruses and some are virulence factors. However, we found that some NS proteins encoded by ZIKV genome appeared to function against the viral replication. In this report we showed that exogenously expressed ZIKV NS2A and NS4A inhibited ZIKV infection by inhibiting viral RNA replication in microglial cells and astrocytes. To understand how viral NS proteins suppressed viral replication, we analyzed the transcriptome of the microglial cells and astrocytes and found that expression of NS4A induced the upregulation of ISGs, including MX1/2, OAS1/2/3, IFITM1, IFIT1, IFI6, IFI27, ISG15 or BST2 through activating the ISGF3 signaling pathway. Upregulation of these ISGs seemed to be related to the inhibition of ZIKV replication, since the anti-ZIKV function of NS4A was partially attenuated when the cells were treated with Abrocitinib, an inhibitor of the ISGF3 signaling pathway, or were knocked down with STAT2. Aborting the protein expression of NS4A, but not its nucleic acid, eliminated the antiviral activity of NS4A effectively. Dynamic expression of viral NS proteins was examined in ZIKV-infected microglial cells and astrocytes, which showed comparatively NS4A occurred later than other NS proteins during the infection. We hypothesize that NS4A may possess intrinsic features to serve as a unique type of pathogen associated molecular pattern (PAMP), detectable by the cells to induce an innate immune response, or function with other mechanisms, to restrict the viral replication to a certain level as a negative feedback, which may help ZIKV maintain its persistent infection in fetal neural tissues.

Modulation of Zika virus replication via glycosphingolipids

The enveloped positive-sense RNA viruses including Zika virus (ZIKV) need host lipids to successfully replicate. The nature of the lipids and the replication step(s) where lipids are utilized often vary amongst viruses. In this study, we demonstrate that ZIKV particle envelope is significantly enriched in distinct sphingolipid species. To determine the role of sphingolipids in ZIKV replication, we leveraged a panel of sphingolipid-deficient cell lines. Notably, knockout of glucosylceramide and lactosylceramide synthase encoding genes (GCSKO; B4G5KO) resulted in a marked decrease in ZIKV titers. GCSKO or pharmacological inhibition of GCS also led to a significant decrease in ZIKV genome replication. Further analysis indicated that GCSKO reduced intracellular virus titers but had minimal impact on ZIKV binding. Restoration of B4G5 expression in B4G5KO cells or supplementing PDMP-treated cells with glucosylceramide led to a significant rescue of ZIKV replication. Altogether, our findings suggest that ZIKV needs glycosphingolipids to facilitate virus replication.

A Retrospective Overview of Zika Virus Evolution in the Midwest of Brazil.

ABSTRACTSince the introduction of the Zika virus (ZIKV) into Brazil in 2015, its transmission dynamics have been intensively studied in many parts of the country, although much is still unknown about its circulation in the midwestern states. Here, using nanopore technology, we obtained 23 novel partial and near-complete ZIKV genomes from the state of Goiás, located in the Midwest of Brazil. Genomic, phylogenetic, and epidemiological approaches were used to retrospectively explore the spatiotemporal evolution of the ZIKV-Asian genotype in this region. As a likely consequence of a gradual accumulation of herd immunity, epidemiological data revealed a decline in the number of reported cases over 2018 to 2021. Phylogenetic reconstructions revealed that multiple independent introductions of the Asian lineage have occurred in Goiás over time and revealed a complex transmission dynamic between epidemic seasons. Together, our results highlight the utility of genomic, epidemiological, and evolutionary methods to understand mosquito-borne epidemics.IMPORTANCE Despite the considerable morbidity and mortality of arboviral infections in Brazil, such as Zika, chikungunya, dengue fever, and yellow fever, our understanding of these outbreaks is hampered by the limited availability of genomic data to track and control the epidemic. In this study, we provide a retrospective reconstruction of the Zika virus transmission dynamics in the state of Goiás by analyzing genomic data from areas in Midwest Brazil not covered by other previous studies. Our study provides an understanding of how ZIKV initiates transmission in this region and reveals a complex transmission dynamic between epidemic seasons. Together, our results highlight the utility of genomic, epidemiological, and evolutionary methods to understand mosquito-borne epidemics, revealing how this toolkit can be used to help policymakers prioritize areas to be targeted, especially in the context of finite public health resources.