Structural and Nonstructural Genes Contribute to the Genetic Diversity of RNA Viruses

Unraveling flavivirus pathogenesis: from bulk to single-cell RNA-sequencing strategies

The global spread of flaviviruses has triggered major outbreaks worldwide, significantly impacting public health, society, and economies. This has intensified research efforts to understand how flaviviruses interact with their hosts and manipulate the immune system, underscoring the need for advanced research tools. RNA-sequencing (RNA-seq) technologies have revolutionized our understanding of flavivirus infections by offering transcriptome analysis to dissect the intricate dynamics of virus-host interactions. Bulk RNA-seq provides a macroscopic overview of gene expression changes in virus-infected cells, offering insights into infection mechanisms and host responses at the molecular level. Single-cell RNA sequencing (scRNAseq) provides unprecedented resolution by analyzing individual infected cells, revealing remarkable cellular heterogeneity within the host response. A particularly innovative advancement, virus-inclusive single-cell RNA sequencing (viscRNA-seq), addresses the challenges posed by non-polyadenylated flavivirus genomes, unveiling intricate details of virus-host interactions. In this review, we discuss the contributions of bulk RNA-seq, scRNA-seq, and viscRNA-seq to the field, exploring their implications in cell line experiments and studies on patients infected with various flavivirus species. Comprehensive transcriptome analyses from RNA-seq technologies are pivotal in accelerating the development of effective diagnostics and therapeutics, paving the way for innovative treatments and enhancing our preparedness for future outbreaks.

Increased interregional virus exchange and nucleotide diversity outline the expansion of chikungunya virus in Brazil

The emergence and reemergence of mosquito-borne diseases in Brazil such as yellow fever, zika, chikungunya, and dengue have had serious impacts on public health. Concerns have been raised due to the rapid dissemination of the chikungunya virus across the country since its first detection in 2014 in Northeast Brazil. In this work, we carried out on-site training activities in genomic surveillance in partnership with the National Network of Public Health Laboratories that have led to the generation of 422 chikungunya virus genomes from 12 Brazilian states over the past two years (2021-2022), a period that has seen more than 312 thousand chikungunya fever cases reported in the country. These genomes increased the amount of available data and allowed a more comprehensive characterization of the dispersal dynamics of the chikungunya virus East-Central-South-African lineage in Brazil. Tree branching patterns revealed the emergence and expansion of two distinct subclades. Phylogeographic analysis indicated that the northeast region has been the leading hub of virus spread towards other regions. Increased frequency of C > T transitions among the new genomes suggested that host restriction factors from the immune system such as ADAR and AID/APOBEC deaminases might be driving the genetic diversity of the chikungunya virus in Brazil.

Differential Genome Replication of a Unique Single-Amino-Acid Mutation in the Adenovirus-4 Component of the Live Oral Adenovirus Type 4 and Type 7 Vaccine

The FDA-approved Adenovirus Type 4 and Type 7 Vaccine, Live, Oral is highly effective and essential in preventing acute respiratory diseases (ARDs) in U.S. military recruits. Our study revealed the presence of a previously undetected mutation, not found in the wild-type human adenovirus type 4 (HAdV-4) component of the licensed vaccine, which contains an amino acid substitution (P388T) in the pre-terminal protein (pTP). This study demonstrated that replication of the T388 HAdV-4 vaccine mutant virus is favored over the wild type in WI-38 cells, the cell type utilized in vaccine manufacturing. However, results from serial human stool specimens of vaccine recipients support differential genome replication in the gastrointestinal tract (GI), demonstrated by the steady decline of the percentage of mutant T388 vaccine virus. Since vaccine efficacy depends upon GI replication and the subsequent immune response, the mutation can potentially impact vaccine efficacy.

Evaluation of NS4A, NS4B, NS5 and 3'UTR Genetic Determinants of WNV Lineage 1 Virulence in Birds and Mammals

West Nile virus (WNV) is amplified in an enzootic cycle involving birds as amplifying hosts. Because they do not develop high levels of viremia, humans and horses are considered to be dead-end hosts. Mosquitoes, especially from the Culex genus, are vectors responsible for transmission between hosts. Consequently, understanding WNV epidemiology and infection requires comparative and integrated analyses in bird, mammalian, and insect hosts. So far, markers of WNV virulence have mainly been determined in mammalian model organisms (essentially mice), while data in avian models are still missing. WNV Israel 1998 (IS98) is a highly virulent strain that is closely genetically related to the strain introduced into North America in 1999, NY99 (genomic sequence homology > 99%). The latter probably entered the continent at New York City, generating the most impactful WNV outbreak ever documented in wild birds, horses, and humans. In contrast, the WNV Italy 2008 strain (IT08) induced only limited mortality in birds and mammals in Europe during the summer of 2008. To test whether genetic polymorphism between IS98 and IT08 could account for differences in disease spread and burden, we generated chimeric viruses between IS98 and IT08, focusing on the 3' end of the genome (NS4A, NS4B, NS5, and 3'UTR regions) where most of the non-synonymous mutations were detected. In vitro and in vivo comparative analyses of parental and chimeric viruses demonstrated a role for NS4A/NS4B/5'NS5 in the decreased virulence of IT08 in SPF chickens, possibly due to the NS4B-E249D mutation. Additionally, significant differences between the highly virulent strain IS98 and the other three viruses were observed in mice, implying the existence of additional molecular determinants of virulence in mammals, such as the amino acid changes NS5-V258A, NS5-N280K, NS5-A372V, and NS5-R422K. As previously shown, our work also suggests that genetic determinants of WNV virulence can be host-dependent.

Increased interregional virus exchange and nucleotide diversity outline the expansion of the chikungunya virus ECSA lineage in Brazil

The emergence and reemergence of mosquito-borne diseases in Brazil such as Yellow Fever, Zika, Chikungunya, and Dengue have had serious impacts on public health. Concerns have been raised due to the rapid dissemination of the chikungunya virus (CHIKV) across the country since its first detection in 2014 in Northeast Brazil. Faced with this scenario, on-site training activities in genomic surveillance carried out in partnership with the National Network of Public Health Laboratories have led to the generation of 422 CHIKV genomes from 12 Brazilian states over the past two years (2021-2022), a period that has seen more than 312 thousand chikungunya fever cases reported in the country. These new genomes increased the amount of available data and allowed a more comprehensive characterization of the dispersion dynamics of the CHIKV East-Central-South-African (ECSA) lineage in Brazil. Tree branching patterns revealed the emergence and expansion of two distinct subclades. Phylogeographic analysis indicated that the northeast region has been the leading hub of virus spread towards other regions. Increased frequency of C>T transitions among the new genomes suggested that host restriction factors from the immune system such as ADAR and AID/APOBEC deaminases might be driving CHIKV ECSA lineage genetic diversity in Brazil.

Viral intra-host evolutionary dynamics revealed via serial passage of Japanese encephalitis virus in vitro

Analyses of viral inter- and intra-host mutations could better guide the prevention and control of infectious diseases. For a long time, studies on viral evolution have focused on viral inter-host variations. Next-generation sequencing has accelerated the investigations of viral intra-host diversity. However, the theoretical basis and dynamic characteristics of viral intra-host mutations remain unknown. Here, using serial passages of the SA14-14-2 vaccine strain of Japanese encephalitis virus (JEV) as the in vitro model, the distribution characteristics of 1,788 detected intra-host single-nucleotide variations (iSNVs) and their mutated frequencies from 477 deep-sequenced samples were analyzed. Our results revealed that in adaptive (baby hamster kidney (BHK)) cells, JEV is under a nearly neutral selection pressure, and both non-synonymous and synonymous mutations represent an S-shaped growth trend over time. A higher positive selection pressure was observed in the nonadaptive (C6/36) cells, and logarithmic growth in non-synonymous iSNVs and linear growth in synonymous iSNVs were observed over time. Moreover, the mutation rates of the NS4B protein and the untranslated region (UTR) of the JEV are significantly different between BHK and C6/36 cells, suggesting that viral selection pressure is regulated by different cellular environments. In addition, no significant difference was detected in the distribution of mutated frequencies of iSNVs between BHK and C6/36 cells.

Computational Analysis Predicts Correlations among Amino Acids in SARS-CoV-2 Proteomes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a serious global challenge requiring urgent and permanent therapeutic solutions. These solutions can only be engineered if the patterns and rate of mutations of the virus can be elucidated. Predicting mutations and the structure of proteins based on these mutations have become necessary for early drug and vaccine design purposes in anticipation of future viral mutations. The amino acid composition (AAC) of proteomes and individual viral proteins provide avenues for exploitation since AACs have been previously used to predict structure, shape and evolutionary rates. Herein, the frequency of amino acid residues found in 1637 complete proteomes belonging to 11 SARS-CoV-2 variants/lineages were analyzed. Leucine is the most abundant amino acid residue in the SARS-CoV-2 with an average AAC of 9.658% while tryptophan had the least abundance of 1.11%. The AAC and ranking of lysine and glycine varied in the proteome. For some variants, glycine had higher frequency and AAC than lysine and vice versa in other variants. Tryptophan was also observed to be the most intolerant to mutation in the various proteomes for the variants used. A correlogram revealed a very strong correlation of 0.999992 between B.1.525 (Eta) and B.1.526 (Iota) variants. Furthermore, isoleucine and threonine were observed to have a very strong negative correlation of -0.912, while cysteine and isoleucine had a very strong positive correlation of 0.835 at p < 0.001. Shapiro-Wilk normality test revealed that AAC values for all the amino acid residues except methionine showed no evidence of non-normality at p < 0.05. Thus, AACs of SARS-CoV-2 variants can be predicted using probability and z-scores. AACs may be beneficial in classifying viral strains, predicting viral disease types, members of protein families, protein interactions and for diagnostic purposes. They may also be used as a feature along with other crucial factors in machine-learning based algorithms to predict viral mutations. These mutation-predicting algorithms may help in developing effective therapeutics and vaccines for SARS-CoV-2.

The second coming of urban yellow fever in the Americas: looking the past to see the future

Yellow fever (YF) epizootics in South America during the 21st century have an unprecedented recorded magnitude and geographical dispersion. The YF spread progressively involved areas considered previously free of risk reaching the edge of cities with large unvaccinated populations, and urban outbreaks were frequently forecasted. We critically reviewed the initial stages and enhancing contexts of YF urban epidemics since the 17th century in the Americas, and the modeling attempts of YF epidemic risk by of Aedes-Human transmission, to find common factors that increase the probability of these events in the current scenarios. The YF urban outbreaks of the past showed as necessary conditions the multiple introduction by viremic carriers clustered in time and space, coincident with population peaks of Aedes. These conditions are not met in the current outbreaks in the Americas by sylvatic YF cycles, besides the protective impact of vector control campaigns, vaccination coverage, improved surveillance, and case management. Therefore, urban Aedes-Human YF outbreaks in the Americas are still possible but with low probability or very focal transmission, while the conditions reported in the past were avoided, and the surveillance and control measures sustained, including the vaccination of the population at risk.

Molecular Characterization and Phylogenetic Analysis of the 2019 Dengue Outbreak in Wenzhou, China

In 2019, a dengue outbreak occurred with 290 confirmed cases in Wenzhou, a coastal city in southeast China. To identify the origin of the dengue virus (DENV) from this outbreak, viral RNA was extracted from four serum samples and sequenced for whole genome analysis. Then, phylogenetic analysis, gene mutation, secondary structure prediction, selection pressure analysis, and recombination analysis were performed. DENV strains Cam-03 and Cam-11 were isolated from patients traveling from Cambodia, while ZJWZ-18 and ZJWZ-62 strains were isolated from local patients without a record of traveling abroad. The whole genome sequence of all four strains was 10,735 nucleotides long. Phylogenetic tree analysis showed that the four strains belonged to genotype 1 of DENV-1, but the local Wenzhou strains and imported strains clustered in different branches. ZJWZ-18 and ZJWZ-62 were closely related to strain MF033254-Singapore-2016, and Cam-03 and Cam-11 were closely related to strain AB608788-China : Taiwan-1994. A comparison of the coding regions between the local strains and the DENV-1 standard strain (EU848545-Hawaii-1944) showed 82 amino acid mutations between the two strains. A total of 55 amino acid mutations were found between the coding regions of the local and imported strains. The overall secondary structure of the 3' UTR of the local strains had changed: apparent changes in the head and tail position were observed when compared to DENV-1 standard strain. Furthermore, selection pressure analysis and recombination detection using the 4 isolates and 41 reference strains showed two credible positive selection sites and eight credible recombination events, which warrant further studies. This study may enhance the understanding of viral replication, infection, evolution, virulence, and pathogenicity of DENV.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Impact of yellow fever virus envelope protein on wild-type and vaccine epitopes and tissue tropism

The envelope (E) protein of flaviviruses is functionally associated with viral tissue tropism and pathogenicity. For yellow fever virus (YFV), viscerotropic disease primarily involving the liver is pathognomonic for wild-type (WT) infection. In contrast, the live-attenuated vaccine (LAV) strain 17D does not cause viscerotropic disease and reversion to virulence is associated with neurotropic disease. The relationship between structure-function of the E protein for WT strain Asibi and its LAV derivative 17D strain is poorly understood; however, changes to WT and vaccine epitopes have been associated with changes in virulence. Here, a panel of Asibi and 17D infectious clone mutants were generated with single-site mutations at the one membrane residue and each of the eight E protein amino acid substitutions that distinguish the two strains. The mutants were characterized with respect to WT-specific and vaccine-specific monoclonal antibodies (mAbs) binding to virus plus binding of virus to brain, liver, and lung membrane receptor preparations (MRPs) generated from AG129 mice. This approach shows that amino acids in the YFV E protein domains (ED) I and II contain the WT E protein epitope, which overlap with those that mediate YFV binding to mouse liver. Furthermore, amino acids in EDIII associated with the vaccine epitope overlap with those that facilitate YFV binding mouse brain MRPs. Taken together, these data suggest that the YFV E protein is a key determinant in the phenotype of WT and 17D vaccine strains of YFV.

Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus-host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10-3 to 10-6 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3' to 5' exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses.

Japanese encephalitis virus live attenuated vaccine strains display altered immunogenicity, virulence and genetic diversity

Japanese encephalitis virus (JEV) is the etiological agent of Japanese encephalitis (JE). The most commonly used vaccine used to prevent JE is the live-attenuated strain SA14-14-2, which was generated by serial passage of the wild-type (WT) JEV strain SA14. Two other vaccine candidates, SA14-5-3 and SA14-2-8 were derived from SA14. Both were shown to be attenuated but lacked sufficient immunogenicity to be considered effective vaccines. To better contrast the SA14-14-2 vaccine with its less-immunogenic counterparts, genetic diversity, ribavirin sensitivity, mouse virulence and mouse immunogenicity of the three vaccines were investigated. Next generation sequencing demonstrated that SA14-14-2 was significantly more diverse than both SA14-5-3 and SA14-2-8, and was slightly less diverse than WT SA14. Notably, WT SA14 had unpredictable levels of diversity across its genome whereas SA14-14-2 is highly diverse, but genetic diversity is not random, rather the virus only tolerates variability at certain residues. Using Ribavirin sensitivity in vitro, it was found that SA14-14-2 has a lower fidelity replication complex compared to SA14-5-3 and SA14-2-8. Mouse virulence studies showed that SA14-2-8 was the most virulent of the three vaccine strains while SA14-14-2 had the most favorable combination of safety (virulence) and immunogenicity for all vaccines tested. SA14-14-2 contains genetic diversity and sensitivity to the antiviral Ribavirin similar to WT parent SA14, and this genetic diversity likely explains the (1) differences in genomic sequences reported for SA14-14-2 and (2) the encoding of major attenuation determinants by the viral E protein.

Genome Characterization of Yellow Fever Virus Wild-Type Strain Asibi, Parent to Live-Attenuated 17D Vaccine, from Three Different Sources

The yellow fever virus vaccine, 17D, was derived through the serial passage of the wild-type (WT) strain Asibi virus in mouse and chicken tissue. Since its derivation, the mechanism of attenuation of 17D virus has been investigated using three 17D substrains and WT Asibi virus. Although all three substrains of 17D have been sequenced, only one isolate of Asibi has been examined genetically and all interpretation of attenuation is based on this one isolate. Here, we sequenced the genome of Asibi virus from three different laboratories and show that the WT strain is genetically homogenous at the amino acids that distinguish Asibi from 17D vaccine virus.

Flavivirus NS1 and Its Potential in Vaccine Development

The Flavivirus genus contains many important human pathogens, including dengue, Japanese encephalitis (JE), tick-borne encephalitis (TBE), West Nile (WN), yellow fever (YF) and Zika (ZIK) viruses. While there are effective vaccines for a few flavivirus diseases (JE, TBE and YF), the majority do not have vaccines, including WN and ZIK. The flavivirus nonstructural 1 (NS1) protein has an unusual structure-function because it is glycosylated and forms different structures to facilitate different roles intracellularly and extracellularly, including roles in the replication complex, assisting in virus assembly, and complement antagonism. It also plays a role in protective immunity through antibody-mediated cellular cytotoxicity, and anti-NS1 antibodies elicit passive protection in animal models against a virus challenge. Historically, NS1 has been used as a diagnostic marker for the flavivirus infection due to its complement fixing properties and specificity. Its role in disease pathogenesis, and the strong humoral immune response resulting from infection, makes NS1 an excellent target for inclusion in candidate flavivirus vaccines.

Measuring Alphavirus Fidelity Using Non-Infectious Virus Particles

Mutations are incorporated into the genomes of RNA viruses at an optimal frequency and altering this precise frequency has been proposed as a strategy to create live-attenuated vaccines. However, determining the effect of specific mutations that alter fidelity has been difficult because of the rapid selection of the virus population during replication. By deleting residues of the structural polyprotein PE2 cleavage site, E3Δ56-59, in Venezuelan equine encephalitis virus (VEEV) TC-83 vaccine strain, non-infectious virus particles were used to assess the effect of single mutations on mutation frequency without the interference of selection that results from multiple replication cycles. Next-generation sequencing analysis revealed a significantly lower frequency of transversion mutations and overall mutation frequency for the fidelity mutants compared to VEEV TC-83 E3Δ56-59. We demonstrate that deletion of the PE2 cleavage site halts virus infection while making the virus particles available for downstream sequencing. The conservation of the site will allow the evaluation of suspected fidelity mutants across alphaviruses of medical importance.

Limited evolution of the yellow fever virus 17d in a mouse infection model

By infecting mice with the yellow fever virus vaccine strain 17D (YFV-17D; Stamaril®), the dose dependence and evolutionary consequences of neurotropic yellow fever infection was assessed. Highly susceptible AG129 mice were used to allow for a maximal/unlimited expansion of the viral populations. Infected mice uniformly developed neurotropic disease; the virus was isolated from their brains, plaque purified and sequenced. Viral RNA populations were overall rather homogenous [Shannon entropies 0−0.15]. The remaining, yet limited intra-host population diversity (0−11 nucleotide exchanges per genome) appeared to be a consequence of pre-existing clonal heterogeneities (quasispecies) of Stamaril®. In parallel, mice were infected with a molecular clone of YFV-17D which was in vivo launched from a plasmid. Such plasmid-launched YFV-17D had a further reduced and almost clonal evolution. The limited intra-host evolution during unrestricted expansion in a highly susceptible host is relevant for vaccine and drug development against flaviviruses in general. Firstly, a propensity for limited evolution even upon infection with a (very) low inoculum suggests that fractional dosing as implemented in current YF-outbreak control may pose only a limited risk of reversion to pathogenic vaccine-derived virus variants. Secondly, it also largely lowers the chance of antigenic drift and development of resistance to antivirals.

A tug-of-war between severe acute respiratory syndrome coronavirus 2 and host antiviral defence: lessons from other pathogenic viruses

World Health Organization has declared the ongoing outbreak of coronavirus disease 2019 (COVID-19) a Public Health Emergency of International Concern. The virus was named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by the International Committee on Taxonomy of Viruses. Human infection with SARS-CoV-2 leads to a wide range of clinical manifestations ranging from asymptomatic, mild, moderate to severe. The severe cases present with pneumonia, which can progress to acute respiratory distress syndrome. The outbreak provides an opportunity for real-time tracking of an animal coronavirus that has just crossed species barrier to infect humans. The outcome of SARS-CoV-2 infection is largely determined by virus-host interaction. Here, we review the discovery, zoonotic origin, animal hosts, transmissibility and pathogenicity of SARS-CoV-2 in relation to its interplay with host antiviral defense. A comparison with SARS-CoV, Middle East respiratory syndrome coronavirus, community-acquired human coronaviruses and other pathogenic viruses including human immunodeficiency viruses is made. We summarize current understanding of the induction of a proinflammatory cytokine storm by other highly pathogenic human coronaviruses, their adaptation to humans and their usurpation of the cell death programmes. Important questions concerning the interaction between SARS-CoV-2 and host antiviral defence, including asymptomatic and presymptomatic virus shedding, are also discussed.

Whole genome sequencing and genetic variations in several dengue virus type 1 strains from unusual dengue epidemic of 2017 in Vietnam

Background

Dengue hemorrhagic fever is an acute viral infection transmitted by mosquitoes. In the 2017, a dengue epidemic occurred in Hanoi in a short time interval and many cases were serious with associated mortality. This was the largest and unusual dengue fever outbreak in the North of Vietnam over the past 20 years. The objective of the present study was to understand the genetic characteristics of the DENV-1 strain in the 2017 epidemic and its relationship with previous viruses in Vietnam and the rest of the world.

Methods

Complete genomes of 72 DENV-1 from patients in the 2017 epidemic were sequenced using NGS. The full genome sequences were then analyzed to find out the genetic variants in the groups of 72 strains, followed by their comparison with other strains that caused disease in Vietnam previously and several other regions of the world, revealing a genetic relationship between them.

Results

The complete genome sequence of 72 DENV-1 strains comprised 10,697 nucleotides with an open reading frame coding for 3392 amino acids. The genomic analysis revealed different amino acid substitutions in all genes, especially varying at position S75 (Capsid), M125 (PrM), D54 (E), T147, V180 (NS1), G45, Y126, I154 (NS2A), A94 (NS2B), M298 (NS3), K47, V68 (NS4A), I29 (NS4B), and R166, E536, G614, T821 (NS5). The genetic analysis suggested that the viruses were most closely related to the causative virus of the dengue outbreak in Vietnam and Cambodia from 2006 to 2008. These results indicated that DENV-1 from the dengue epidemic 2017 in Northern Vietnam originated from the virus that caused the dengue outbreak during the 2007 to 2008 period in Vietnam.

Conclusion

The present study is the first of its kind to describe complete genome sequence as well as genetic variants and phylogenetic analysis of DENV-1 associated with the unusual dengue epidemic of 2017 in northern Vietnam. These results provide detailed evidence to elucidate the origin, circulation, and genetic evolution of DENV in Vietnam.

Quantitative Analysis of the Substrate Specificity of Human Rhinovirus 3C Protease and Exploration of Its Substrate Recognition Mechanisms

Human rhinovirus 3C protease (HRV 3C-P) is a high-value commercial cysteine protease that could specifically recognize the short peptide sequence of LEVLFQ↓GP. In here, a strategy based on our previous Yeast Endoplasmic Reticulum Sequestration Screening (YESS) approach was developed in Saccharomyces cerevisiae, a model microorganism, to fully characterize the substrate specificity of a typical human virus protease, HRV 3C-P, in a quantitative and fast manner. Our results demonstrated that HRV 3C-P had very high specificity at P1 and P1' positions, only recognizing Gln/Glu at the P1 position and Gly/Ala/Cys/Ser at the P1' position, respectively. Comparably, it exhibited efficient recognition of most residues at the P2' position, except Trp. Further biochemical characterization through site mutagenesis, enzyme structural modeling, and comparison with other 3C proteases indicated that the S1 pocket of HRV 3C-P was constituted by neutral and basic amino acids, in which His160 and Thr141 specifically interacted with Gln or Glu residues at the substrate P1 position. Additionally, the stringent S1' pocket determined its unique property of only accommodating residues without or with short side chains. Based on our characterization, LEVLFQ↓GM was identified as a more favorable substrate than the original LEVLFQ↓GP at high temperature, which might be caused by the conversion of random coils to β-turns in HRV 3C-P along with the temperature increase. Our studies prompted a further understanding of the substrate specificity and recognition mechanism of HRV 3C-P. Besides, the YESS-PSSC combined with the enzyme modeling strategy in this study provides a general strategy for deciphering the substrate specificities of proteases.

Transcriptomic immune profiles of human flavivirus-specific T-cell responses

The Flavivirus genus of viruses includes dengue (DENV), Zika (ZIKV), yellow fever (YFV), Japanese encephalitis (JEV), and West Nile (WNV) viruses. Infections with these species combined are prevalent in tropical and sub-tropical areas, affecting millions of people and ranging from asymptomatic to severe forms of the disease. They therefore pose a serious threat to global public health. Several studies imply a role for T cells in the protection but also pathogenesis against the different flavivirus species. Identifying flavivirus-specific T-cell immune profiles and determining how pre-exposure of one species might affect the immune response against subsequent infections from other species is important to further define the role of T cells in the immune response against infection. Understanding the immune profiles of the flavivirus-specific T-cell response in natural infection is important to understand the T-cell response in the context of vaccination. In this review, we summarize the current knowledge on human immune profiles of flavivirus-specific T-cell reactivity, comparing natural infection with the acute form of the disease and vaccination in different flavivirus infections.

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

Zika virus (ZIKV) strains belong to the East African, West African, and Asian/American phylogenetic lineages. RNA viruses, like ZIKV, exist as populations of genetically-related sequences whose heterogeneity may impact viral fitness, evolution, and virulence. Genetic diversity of representative ZIKVs from each lineage was examined using next generation sequencing (NGS) paired with downstream entropy and single nucleotide variant (SNV) analysis. Comparisons showed that inter-lineage diversity was statistically supported, while intra-lineage diversity. Intra-lineage diversity was significant for East but not West Africa strains. Furthermore, intra-lineage diversity for the Asian/American lineage was not supported for human serum isolates; however, a placenta isolate differed significantly. Relative in the pre-membrane/membrane (prM/M) gene of several ZIKV strains. Additionally, the East African lineage contained a greater number of synonymous SNVs, while a greater number of non-synonymous SNVs were identified for American strains. Further, inter-lineage SNVs were dispersed throughout the genome, whereas intra-lineage non-synonymous SNVs for Asian/American strains clustered within prM/M and NS1 gene. This comprehensive analysis of ZIKV genetic diversity provides a repository of SNV positions across lineages. We posit that increased non-synonymous SNV populations and increased relative genetic diversity of the prM/M and NS1 proteins provides more evidence for their role in ZIKV virulence and fitness.

Japanese encephalitis vaccine-specific envelope protein E138K mutation does not attenuate virulence of West Nile virus

West Nile (WNV) and Japanese encephalitis viruses (JEV) are closely related, mosquito-borne neurotropic flaviviruses. Although there are no licensed human vaccines for WNV, JEV has multiple human vaccines, including the live, attenuated vaccine SA14-14-2. Investigations into determinants of attenuation of JE SA14-14-2 demonstrated that envelope (E) protein mutation E138K was crucial to the attenuation of mouse virulence. As WNV is closely related to JEV, we investigated whether or not the E-E138K mutation would be beneficial to be included in a candidate live attenuated WNV vaccine. Rather than conferring a mouse attenuated phenotype, the WNV E-E138K mutant reverted and retained a wild-type mouse virulence phenotype. Next-generation sequencing analysis demonstrated that, although the consensus sequence of the mutant had the E-E138K mutation, there was increased variation in the E protein, including a single-nucleotide variant (SNV) revertant to the wild-type glutamic acid residue. Modeling of the E protein and analysis of SNVs showed that reversion was likely due to the inability of critical E-protein residues to be compatible electrostatically. Therefore, this mutation may not be reliable for inclusion in candidate live attenuated vaccines in related flaviviruses, such as WNV, and care must be taken in translation of attenuating mutations from one virus to another virus, even if they are closely related.

Japanese encephalitis virus (JEV) and West Nile virus (WNV) are closely related neurotropic viruses—a live attenuated vaccine exists for JEV but not for WNV. A team led by Alan D.T. Barrett at the University of Texas investigated whether a key E-protein mutation (E138K) in the live attenuated JEV vaccine can also attenuate a candidate live WNV vaccine. The mutant WNV vaccine shows essentially identical behavior to the virulent parental strain in vitro but unexpectedly also has unimpaired lethality and neurotropism when mice are challenged intraperitoneally. Sequencing of the mutant WNV vaccine demonstrated increased frequencies of single-nucleotide variants clustered around residue 138—including reversion to the parental strain glutamic acid. E138K mutation is therefore unlikely to be a reliable means to attenuate candidate WNV vaccines.

Attenuation of Live-Attenuated Yellow Fever 17D Vaccine Virus Is Localized to a High-Fidelity Replication Complex

Live-attenuated viral vaccines are highly safe and efficacious but represent complex and often multigenic attenuation mechanisms. Most of these vaccines have been generated empirically by serial passaging of a wild-type (WT) virus in cell culture. One of the safest and most effective live-attenuated vaccines is yellow fever (YF) virus strain 17D, which has been used for over 80 years to control YF disease. The availability of the WT parental strain of 17D, Asibi virus, and large quantities of clinical data showing the effectiveness of the 17D vaccine make this WT parent/vaccine pair an excellent model for investigating RNA virus attenuation. Here, we investigate a mechanism of 17D attenuation and show that the vaccine virus is resistant to the antiviral compound ribavirin. The findings suggest that attenuation is in part due to a low probability of reversion or mutation of the vaccine virus genome to WT, thus maintaining a stable genotype despite external pressures.

The molecular basis of attenuation for live-attenuated vaccines is poorly understood. The yellow fever (YF) 17D vaccine virus was derived from the wild-type, parental strain Asibi virus by serial passage in chicken tissue and has proven to be a very safe and efficacious vaccine. We have previously shown that wild-type Asibi is a typical RNA virus with high genetic diversity, while the 17D vaccine virus has very little genetic diversity. To investigate this further, we treated Asibi and 17D viruses with ribavirin, a GTP analog with strong antiviral activity that increases levels of mutations in the viral genome. As expected, ribavirin treatment introduced mutations into the Asibi virus genome at a very high frequency and decreased viral infectivity while, in contrast, the 17D vaccine virus was resistant to ribavirin, as treatment with the antiviral introduced very few mutations into the genome, and viral infectivity was not lost. The results were confirmed for another YF wild-type parental and vaccine pair, a wild-type French viscerotropic virus and French neurotropic vaccine. Using recombinant Asibi and 17D viruses, ribavirin sensitivity was located to viral nonstructural genes. Thus, two live-attenuated YF vaccine viruses are genetically stable even under intense mutagenic pressure, suggesting that attenuation of live-attenuated YF vaccines is due, at least in part, to fidelity of the replication complex resulting in high genetic stability.

Genotypic and phenotypic characterization of West Nile virus NS5 methyltransferase mutants

Although West Nile virus (WNV) causes annual cases of neurological disease and deaths in humans, a vaccine has not been licensed for human use. Several WNV genes have been targeted for mutagenesis in attempts to generate live attenuated vaccine candidates, including the non-structural protein NS5. Specifically, mutation of WNV NS5-K61A or NS5-E218A in the catalytic tetrad of the methyltransferase decreases enzyme activity of the NS5 protein and correspondingly attenuates the virus in mice. In this report, NS5-K61A, NS5-E218A, and a double mutant encoding both mutations (NS5-K61A/E218A) were compared both in vitro and in vivo. Each single mutant was strongly attenuated in highly susceptible outbred mice, whereas the double mutant unexpectedly was not attenuated. Sequencing analysis demonstrated that the double mutant was capable of reversion at both residues NS5-61 and NS5-218, whereas the genotype of the single mutants did not show evidence of reversion. Overall, either NS5-K61A or NS5-E218A methyltransferase mutations could be potential mutations to include in a candidate live WNV vaccine; however, multiple mutations in the catalytic tetrad of the methyltransferase are not tolerated.

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.