Mutation rates among RNA viruses

Identifying the unknown: Application of molecular epidemiology tools to identify clustering and HIV transmission routes in Poland

BACKGROUND

Understanding the dynamics of HIV-1 transmission is essential for developing effective screening and intervention strategies. Viral genetic sequences provide valuable information that can be used to infer the history and patterns of viral transmission.

PURPOSE

Our study explores the structure and dynamics of HIV transmissions in Poland from 1999 to 2022 to elucidate key patterns related with national epidemics.

METHODS

To understand the temporal dynamics of transmission routes we examined HIV pol sequence data from 5705 Polish PWH. The HIV-TRAnsmission Cluster Engine (HIV-TRACE) was utilized to identify potential links between different risk groups and putative links to individuals with unreported transmission risk.

RESULTS

Our analyses generated 503 clusters, containing 3942 individuals, and identified 13,917 putative links. Approximately 69.1 % of the sequences formed clusters. In the dataset 32.2 % of individuals were reported MSM transmission route, 7.9 % by heterosexual, and 5.6 % by PWID transmissions. The transmission route was unknown for 54.2 % of patients. Putative transmissions from MSM to all other groups revealed that 45.1 % of links lead to people with unregistered transmission mode. For heterosexual patients, 40.2 % of connections were directed to patients lacking information on infection routes and 30.5 % to MSM individuals. Our analysis unveiled that 45.1 % of cases with unreported transmission routes may be identified as MSM, while 3.5 % might be potential non-disclosed MSM.

CONCLUSIONS

Genetic linkages can provide valuable insights into the transmission dynamics among individuals, even in cases where transmission risk information is missing or unreported. The observed association between MSM and unreported cases highlights the potential of molecular epidemiology to complete missing patient data.

SARS-CoV-2 biological clones are genetically heterogeneous and include clade-discordant residues

Defective genomes are part of SARS-CoV-2 quasispecies. High-resolution, ultra-deep sequencing of bulk RNA from viral populations does not distinguish RNA mutations, insertions, and deletions in viable genomes from those in defective genomes. To quantify SARS-CoV-2 infectious variant progeny, virus from four individual plaques (biological clones) of a preparation of isolate USA-WA1/2020, formed on Vero E6 cell monolayers, was subjected to further biological cloning to yield 9 second-generation and 15 third-generation sub-clones. Consensus genomic sequences of the biological clones and sub-clones included an average of 2.8 variations per viable genome, relative to the consensus sequence of the parental USA-WA1/2020 virus. This value is 6.5-fold lower than the estimates for biological clones of other RNA viruses such as bacteriophage Qβ, foot-and-mouth disease virus, or hepatitis C virus in cell culture. The mutant spectrum complexity of the nsp12 (polymerase)- and spike (S)-coding region was unique in the progeny of each of 10 third-generation sub-clones; they shared 2.4% of the total of 164 different mutations and deletions scored in the 3,719 genomic residues that were screened. The presence of minority out-of-frame deletions revealed the ease of defective genome production from an individual infectious genome. Several low-frequency point mutations and deletions were clade-discordant in that they were not typical of USA-WA1/2020 but served to define the consensus sequences of future SARS-CoV-2 clades. Implications for SARS-CoV-2 adaptability and COVID-19 control of the viable genome heterogeneity and the generation of complex mutant spectra from individual genomes are discussed.IMPORTANCESequencing of biological clones is a means to identify mutations, insertions, and deletions located in viable genomes. This distinction is particularly important for viral populations, such as those of SARS-CoV-2, that contain large proportions of defective genomes. By sequencing biological clones and sub-clones, we quantified the heterogeneity of the viable complement of USA-WA1/2020 to be lower than exhibited by other RNA viruses. This difference may be due to a reduced mutation rate or to limited tolerance of the large coronavirus genome to incorporate mutations and deletions and remain functional or a combination of both influences. The presence of clade-discordant residues in the progeny of individual biological sub-clones suggests limitations in the occupation of sequence space by SARS-CoV-2. However, the complex and unique mutant spectra that are rapidly generated from individual genomes suggest an aptness to confront selective constraints.

Receptor Binding for the Entry Mechanisms of SARS-CoV-2: Insights from the Original Strain and Emerging Variants

Since its emergence in late 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continuously evolved, giving rise to multiple variants that have significantly altered the trajectory of the COVID-19 pandemic. These variants have resulted in multiple waves of the pandemic, exhibiting characteristic mutations in the spike (S) protein that may have affected receptor interaction, tissue tropism, and cell entry mechanisms. While the virus was shown to primarily utilize the angiotensin-converting enzyme 2 (ACE2) receptor and host proteases such as transmembrane serine protease 2 (TMPRSS2) for entry into host cells, alterations in the S protein have resulted in changes to receptor binding affinity and use of alternative receptors, potentially expanding the virus's ability to infect different cell types or tissues, contributing to shifts in clinical presentation. These changes have been linked to variations in disease severity, the emergence of new clinical manifestations, and altered transmission dynamics. In this paper, we overview the evolving receptor utilization strategies of SARS-CoV-2, focusing on how mutations in the S protein may have influenced viral entry mechanisms and clinical outcomes across the ongoing pandemic waves.

Whole genome sequencing and molecular detection of potato virus X in Bangladesh

Potato Virus X (PVX) is a significant viral pathogen affecting potato (Solanum tuberosum) crops globally, yet its molecular characterization in Bangladesh remains limited. This study presents the first whole genome sequence (WGS) and molecular analysis of PVX isolated from potato plants in Gazipur, Bangladesh. Initial virus detection was performed using DAS-ELISA on symptomatic potato leaves, followed by RT-PCR targeting the coat protein (CP) gene, which confirmed PVX presence in 'Patnai' and 'Challisha' potato varieties through a 562 bp amplicon. The WGS of the Patnai-PVX isolate was determined to be 6,435 nucleotides long and deposited in GenBank (accession: PQ527059). Genome analysis identified five major open reading frames encoding the RNA-dependent RNA polymerase (RdRp), triple gene block proteins (TGBp1, TGBp2, TGBp3), and CP. Basic Local Alignment Search Tool X (BLASTX) analysis revealed high sequence similarity with PVX isolates from neighboring regions, suggesting evolutionary conservation. Mutation analysis identified 265 SNPs, predominantly synonymous mutations, indicating maintained protein-coding integrity despite genetic variation. Fewer non-synonymous mutations were detected, potentially affecting viral protein functions and pathogenicity. Phylogenetic analysis based on the entire genome sequence placed the Bangladeshi isolate (PQ527059.1) in a well-supported clade (bootstrap value 99%) with isolates from Peru (MT752634.1, MT752612.1, MT752621.1), highlighting potential international transmission routes while also exhibiting unique genetic markers indicative of regional specificity. This comprehensive molecular characterization enhances understanding of PVX genetic diversity and evolution in Bangladesh, providing valuable insights for developing effective virus management strategies in potato cultivation.

Multifractal analysis and support vector machine for the classification of coronaviruses and SARS-CoV-2 variants

This study presents a novel approach for the classification of coronavirus species and variants of SARS-CoV-2 using Chaos Game Representation (CGR) and 2D Multifractal Detrended Fluctuation Analysis (2D MF-DFA). By extracting fractal parameters from CGR images, we constructed a state space that effectively distinguishes different species and variants. Our method achieved [Formula: see text] accuracy in species classification, with a notable [Formula: see text] accuracy for SARS-CoV-2 variants despite their genetic similarities. Using a Support Vector Machine (SVM) as a classifier further enhanced the performance. This approach, which requires fewer steps than most existing methods, offers an efficient and effective tool for viral classification, with implications for bioinformatics, public health, and vaccine development.

The coronavirus nsp14 exoribonuclease interface with the cofactor nsp10 is essential for efficient virus replication and enzymatic activity

Coronaviruses (CoVs) encode non-structural proteins (nsp's) 1-16, which assemble to form replication-transcription complexes that function in viral RNA synthesis. All CoVs encode a proofreading 3'-5' exoribonuclease in non-structural protein 14 (nsp14-ExoN) that mediates proofreading and high-fidelity replication and is critical for other roles in replication and pathogenesis. The in vitro enzymatic activity of nsp14-ExoN is enhanced in the presence of the cofactor nsp10. We introduced alanine substitutions in nsp14 of murine hepatitis virus (MHV) at the nsp14-nsp10 interface and recovered mutant viruses with a range of impairments in replication and in vitro biochemical exonuclease activity. Two of these substitutions, nsp14 K7A and D8A, had impairments intermediate between wild type-MHV nsp14 and the known ExoN(-) D89A/E91A nsp14 catalytic inactivation mutant. All introduced nsp14-nsp10 interface alanine substitutions impaired in vitro exonuclease activity. Passage of the K7A and D8A mutant viruses selected second-site non-synonymous mutations in nsp14 associated with improved mutant virus replication and exonuclease activity. These results confirm the essential role of the nsp14-nsp10 interaction for efficient enzymatic activity and virus replication, identify proximal and long-distance determinants of nsp14-nsp10 interaction, and support targeting the nsp14-nsp10 interface for viral inhibition and attenuation.IMPORTANCECoronavirus replication requires assembly of a replication transcription complex composed of nsp's, including polymerase, helicase, exonuclease, capping enzymes, and non-enzymatic cofactors. The coronavirus nsp14 exoribonuclease mediates several functions in the viral life cycle including genomic and subgenomic RNA synthesis, RNA recombination, RNA proofreading and high-fidelity replication, and native resistance to many nucleoside analogs. The nsp-14 exonuclease activity in vitro requires the non-enzymatic cofactor nsp10, but the determinants and importance of the nsp14-nsp10 interactions during viral replication have not been defined. Here we show that for the coronavirus murine hepatitis virus, nsp14 residues at the nsp14-nsp10 interface are essential for efficient viral replication and in vitro exonuclease activity. These results shed new light on the requirements for protein interactions within the coronavirus replication transcription complex, and they may reveal novel non-active-site targets for virus inhibition and attenuation.

A general and biomedical perspective of viral quasispecies

Viral quasispecies refers to the complex and dynamic mutant distributions (also termed mutant spectra, clouds, or swarms) that arise as a result of high error rates during RNA genome replication. The mutant spectrum of individual RNA virus populations is modified by continuous generation of variant genomes, competition and interactions among them, environmental influences, bottleneck events, and bloc transmission of viral particles. Quasispecies dynamics provides a new perspective on how viruses adapt, evolve, and cause disease, and sheds light on strategies to combat them. Molecular flexibility, together with ample opportunity of mutant cloud traffic in our global world, are key ingredients of viral disease emergences, as exemplified by the recent COVID-19 pandemic. In the present article, we present a brief overview of the molecular basis of mutant swarm formation and dynamics, and how the latter relates to viral disease and epidemic spread. We outline future challenges derived of the highly diverse cellular world in which viruses are necessarily installed.

New N6-Substituted Adenine Derivatives with High Antiviral Activity against RNA-Containing Viruses

In this work, two new compounds, N6-(4,5-dimethoxyphenyl)adenine and N6-(3,5-di-trifluoromethylphenyl)adenine, with a broad range of antiviral activity against RNA viruses were identified. We showed that these compounds exhibit pronounced antiviral activity against human poliovirus types 1, 2, and 3, belonging to enterovirus C species. Both compounds also demonstrated pronounced antiviral activity against Coxsackie viruses B3, B5, and B6, belonging to enterovirus B species. In addition, the compounds demonstrated antiviral activity against Newcastle disease virus, which belongs to the paramyxovirus genus. The compounds discovered in this work can subsequently serve as prototypes for the development of new antiviral drugs against epidemiologically significant human RNA viruses.

The Trajectory of KoRV-A Evolution Indicates Initial Integration into the Koala Germline Genome Near Coffs Harbour

Background

Koala Retrovirus-A is a gamma-retrovirus that is spreading across wild koala populations through horizontal and vertical transmission, contributing significantly to genomic diversity across and even within koala populations. Previous studies have estimated that KoRV-A initially integrated into the koala genome less than 50,000 years ago, but the precise origins and the patterns of spread after its endogenization remain unclear.

Results

In this study, we analyzed germline insertions of KoRV-A using whole-genome sequencing data from 405 wild koalas, representing nearly the species' entire geographic range. Our findings reveal an evolutionary trajectory for KoRV-A, suggesting that the initial endogenization might occur near Coffs Harbour on the Mid-north coast of NSW around the middle of the koala's range. As KoRV-A spread, certain subtypes emerged and became prevalent, two of which recombined with an ancient endogenous retrovirus, PhER, resulting in distinct recombination variants in northern and southern koala populations. Additionally, we identified a geographic barrier north of Sydney, which may have slowed the southward spread of KoRV-A into Sydney and beyond.

Conclusions

Our study proposes a comprehensive evolutionary pathway for KoRV-A, beginning with its initial endogenization near Coffs Harbour and highlighting barriers and diversification events that have shaped its distribution and impact on koala populations.

Effect of 8-Oxo-1,N6-Ethenoadenine Derivatives on the Activity of RNA Polymerases from SARS-CoV-2 and Escherichia coli

Bacterial and viral RNA polymerases are promising targets for the development of new transcription inhibitors. One of the potential blockers of RNA synthesis is 7,8-dihydro-8-oxo-1,N6-ethenoadenine (oxo-εA), a synthetic compound that combines two adenine modifications: 8-oxoadenine and 1,N6-ethenoadenine. In this study, we synthesized oxo-εA triphosphate (oxo-εATP) and showed that it could be incorporated by the RNA-dependent RNA polymerase of SARS-CoV-2 into synthesized RNA opposite template residues A and G in the presence of Mn2+ ions. Escherichia coli RNA polymerase incorporated oxo-εATP opposite A residues in the template DNA strand. The presence of oxo-εA instead of adenine in the template DNA strand completely stopped transcription at the modified nucleotide. At the same time, oxo-εATP did not suppress RNA synthesis by both RNA polymerases in the presence of unmodified nucleotides. Therefore, the oxo-εA modification significantly disrupts nucleotide base pairing during RNA synthesis by RNA polymerases of different classes, and the corresponding nucleotide derivatives cannot be used as potential antiviral or antibacterial transcription inhibitors.

Impact of Different Foot and Mouth Disease Vaccine Schemes in Cross-Neutralization Against Heterologous Serotype O Strains in Cattle

The high antigenic variability of the foot-and-mouth disease virus (FMDV) represents a challenge for developing prophylactic strategies, stressing the need for research into vaccines offering broad protection against a range of virus strains. Here, the heterotypic cross-reaction using different vaccine schemes against serotype O strains was studied, evaluating the impact of revaccination, antigen dose, and incorporation of additional FMDV serotypes. Naïve cattle were immunized with seven distinct FMDV vaccines, receiving three doses of the same formulation at 0, 28, and 56 days post-primary vaccination (dpv). Serum samples were collected up to 70 dpv and tested by a virus-neutralizing test against serotype O strains from a South American lineage and two strains representative of two Asian lineages. Our results showed that vaccines containing the ME-SA topotype O1/Campos strain developed cross-neutralizing responses against the two Asian viruses after the first vaccination. In contrast, significant heterotypic neutralizing antibody titers against the homologous topotype strain were only found after the second vaccination, indicating that the phylogenic relationship may differ from the antigenic profiles for these two viruses. The amount of the O1/Campos strain and the revaccination were essential factors for neutralization against the homologous- and heterologous-type O FMDV viruses. The strain composition of the vaccine was only relevant for cross-neutralization against one of the Asian strains, suggesting potential intra-serotypic divergences for this pattern.

Highly specific SARS-CoV-2 main protease (Mpro) mutations against the clinical antiviral ensitrelvir selected in a safe, VSV-based system

In the SARS-CoV-2 pandemic, the so far two most effective approved antivirals are the protease inhibitors nirmatrelvir, in combination with ritonavir (Paxlovid) and ensitrelvir (Xocova). However, antivirals and indeed all antimicrobial drugs are sooner or later challenged by resistance mutations. Studying such mutations is essential for treatment decisions and pandemic preparedness. At the same time, generating resistant viruses to assess mutants is controversial, especially with pathogens of pandemic potential like SARS-CoV-2. To circumvent gain-of-function research with non-attenuated SARS-CoV-2, a previously developed safe system based on a chimeric vesicular stomatitis virus dependent on the SARS-CoV-2 main protease (VSV-Mpro) was used to select mutations against ensitrelvir. Ensitrelvir is clinically especially relevant due to its single-substance formulation, avoiding drug-drug interactions by the co-formulated CYP3A4 inhibitor ritonavir in Paxlovid. By treating VSV-Mpro with ensitrelvir, highly-specific resistant mutants against this inhibitor were selected, while being still fully or largely susceptible to nirmatrelvir. We then confirmed several ensitrelvir-specific mutants in gold standard enzymatic assays and SARS-CoV-2 replicons. These findings indicate that the two inhibitors can have distinct viral resistance profiles, which could determine treatment decisions.

A Sensitivity and Consistency Comparison Between Next-Generation Sequencing and Sanger Sequencing in HIV-1 Pretreatment Drug Resistance Testing

Next-generation sequencing (NGS) for HIV drug resistance (DR) testing has an increasing number of applications for the detection of low-abundance drug-resistant variants (LA-DRVs) in regard to its features as a quasi-species. However, there is less information on its detection performance in DR detection with NGS. To determine the feasibility of using NGS technology in LA-DRV detection for HIV-1 pretreatment drug resistance, 80 HIV-infected individuals who had never undergone antiretroviral therapy were subjected to both NGS and Sanger sequencing (SS) in HIV-1 drug resistance testing. The results reported in this study show that NGS exhibits higher sensitivity for drug resistance identification than SS at a 5% detection threshold. NGS showed a better consistency compared with that of SS for both protease inhibitors (PIs) and integrase inhibitors (INSTIs), with a figure amounting to more than 90%, but worse consistency in nucleotide reverse transcriptase inhibitors (NRTIs), with a consistency ranging from only 61.25% to 87.50%. The consistency of non-nucleotide reverse transcriptase inhibitors (NNRTIs) between NGS and SS was around 85%. NGS showed the highest sensitivity of 87.0% at a 5% threshold. The application of NGS technology in HIV-1 genotype resistance detection in different populations infected with HIV requires further documentation and validation.

Distinct negative-sense RNA viruses induce a common set of transcripts encoding proteins forming an extensive network

The large group of negative-strand RNA viruses (NSVs) comprises many important pathogens. To identify conserved patterns in host responses, we systematically compared changes in the cellular RNA levels after infection of human hepatoma cells with nine different NSVs of different virulence degrees. RNA sequencing experiments indicated that the amount of viral RNA in host cells correlates with the number of differentially expressed host cell transcripts. Time-resolved differential gene expression analysis revealed a common set of 178 RNAs that are regulated by all NSVs analyzed. A newly developed open access web application allows downloads and visualizations of all gene expression comparisons for individual viruses over time or between several viruses. Most of the genes included in the core set of commonly differentially expressed genes (DEGs) encode proteins that serve as membrane receptors, signaling proteins and regulators of transcription. They mainly function in signal transduction and control immunity, metabolism, and cell survival. One hundred sixty-five of the DEGs encode host proteins from which 47 have already been linked to the regulation of viral infections in previous studies and 89 proteins form a complex interaction network that may function as a core hub to control NSV infections.IMPORTANCEThe infection of cells with negative-strand RNA viruses leads to the differential expression of many host cell RNAs. The differential spectrum of virus-regulated RNAs reflects a large variety of events including anti-viral responses, cell remodeling, and cell damage. Here, these virus-specific differences and similarities in the regulated RNAs were measured in a highly standardized model. A newly developed app allows interested scientists a wide range of comparisons and visualizations.

The coronavirus nsp14 exoribonuclease interface with the cofactor nsp10 is essential for efficient virus replication and enzymatic activity

Coronaviruses (CoVs) encode nonstructural proteins (nsps) 1-16, which assemble to form replication-transcription complexes that function in viral RNA synthesis. All CoVs encode a proofreading 3'-5' exoribonuclease (ExoN) in nsp14 (nsp14-ExoN) that mediates proofreading and high-fidelity replication and is critical for other roles in replication and pathogenesis. The in vitro enzymatic activity of nsp14 ExoN is enhanced in the presence of the cofactor nsp10. We introduced alanine substitutions in nsp14 of murine hepatitis virus (MHV) at the nsp14-10 interface and recovered mutant viruses with a range of impairments in replication and in vitro biochemical exonuclease activity. Two of these substitutions, nsp14 K7A and D8A, had impairments intermediate between WT-MHV nsp14 and the known ExoN(-) D89A/E91A nsp14 catalytic inactivation mutant. All introduced nsp14-10 interface alanine substitutions impaired in vitro exonuclease activity. Passage of the K7A and D8A mutant viruses selected second-site non-synonymous mutations in nsp14 associated with improved mutant virus replication and exonuclease activity. These results confirm the essential role of the nsp14-nsp10 interaction for efficient enzymatic activity and virus replication, identify proximal and long-distance determinants of nsp14-nsp10 interaction, and support targeting the nsp14-10 interface for viral inhibition and attenuation.

IMPORTANCE

Coronavirus replication requires assembly of a replication transcription complex composed of nonstructural proteins (nsp), including polymerase, helicase, exonuclease, capping enzymes, and non-enzymatic cofactors. The coronavirus nsp14 exoribonuclease mediates several functions in the viral life cycle including genomic and subgenomic RNA synthesis, RNA recombination, RNA proofreading and high-fidelity replication, and native resistance to many nucleoside analogs. The nsp-14 exonuclease activity in vitro requires the non-enzymatic co-factor nsp10, but the determinants and importance the nsp14-10 interactions during viral replication have not been defined. Here we show that for the coronavirus murine hepatitis virus, nsp14 residues at the nsp14-10 interface are essential for efficient viral replication and in vitro exonuclease activity. These results shed new light on the requirements for protein interactions within the coronavirus replication transcription complex, and they may reveal novel non active-site targets for virus inhibition and attenuation.

Development of genogroup-specific ELISAs based on the VP1 protein to detect antibodies to GIV and GVI feline norovirus

Feline norovirus (FNoV) is a potential pathogen of feline gastroenteritis and has two genogroups (GIV and GVI). Few epidemiological studies have been conducted on FNoV. We designed two enzyme-linked immunosorbent assays (ELISAs) to identify genogroup-specific FNoV antibodies for serological surveillance. Analysis of sera from cats experimentally infected with FNoV GIV or GVI and from specific-pathogen-free (SPF) cats confirmed that the two recombinant proteins used in the assay react in a genogroup-specific manner. Of the 183 samples tested, 6.6% were positive for GIV and 26.2% were positive for GVI. Antibodies to both FNoV genogroups were detected in sera collected in 2005, seven years before FNoV was first reported.

A comprehensive study of SARS-CoV-2 main protease (Mpro) inhibitor-resistant mutants selected in a VSV-based system

Nirmatrelvir was the first protease inhibitor specifically developed against the SARS-CoV-2 main protease (3CLpro/Mpro) and licensed for clinical use. As SARS-CoV-2 continues to spread, variants resistant to nirmatrelvir and other currently available treatments are likely to arise. This study aimed to identify and characterize mutations that confer resistance to nirmatrelvir. To safely generate Mpro resistance mutations, we passaged a previously developed, chimeric vesicular stomatitis virus (VSV-Mpro) with increasing, yet suboptimal concentrations of nirmatrelvir. Using Wuhan-1 and Omicron Mpro variants, we selected a large set of mutants. Some mutations are frequently present in GISAID, suggesting their relevance in SARS-CoV-2. The resistance phenotype of a subset of mutations was characterized against clinically available protease inhibitors (nirmatrelvir and ensitrelvir) with cell-based, biochemical and SARS-CoV-2 replicon assays. Moreover, we showed the putative molecular mechanism of resistance based on in silico molecular modelling. These findings have implications on the development of future generation Mpro inhibitors, will help to understand SARS-CoV-2 protease inhibitor resistance mechanisms and show the relevance of specific mutations, thereby informing treatment decisions.

Genetic diversity accelerates canine distemper virus adaptation to ferrets

RNA viruses adapt rapidly to new host environments by generating highly diverse genome sets, so-called "quasispecies." Minor genetic variants promote their rapid adaptation, allowing for the emergence of drug-resistance or immune-escape mutants. Understanding these adaptation processes is highly relevant to assessing the risk of cross-species transmission and the safety and efficacy of vaccines and antivirals. We hypothesized that genetic memory within a viral genome population facilitates rapid adaptation. To test this, we investigated the adaptation of the Morbillivirus canine distemper virus to ferrets and compared an attenuated, Vero cell-adapted virus isolate with its recombinant derivative over consecutive ferret passages. Although both viruses adapted to the new host, the reduced initial genetic diversity of the recombinant virus resulted in delayed disease onset. The non-recombinant virus gradually increased the frequencies of beneficial mutations already present at very low frequencies in the input virus. In contrast, the recombinant virus first evolved de novo mutations to compensate for the initial fitness impairments. Importantly, while both viruses evolved different sets of mutations, most mutations found in the adapted non-recombinant virus were identical to those found in a previous ferret adaptation experiment with the same isolate, indicating that mutations present at low frequency in the original virus stock serve as genetic memory. An arginine residue at position 519 in the carboxy terminus of the nucleoprotein shared by all adapted viruses was found to contribute to pathogenesis in ferrets. Our work illustrates the importance of genetic diversity for adaptation to new environments and identifies regions with functional relevance.IMPORTANCEWhen viruses encounter a new host, they can rapidly adapt to this host and cause disease. How these adaptation processes occur remains understudied. Morbilliviruses have high clinical and veterinary relevance and are attractive model systems to study these adaptation processes. The canine distemper virus is of particular interest, as it exhibits a broader host range than other morbilliviruses and frequently crosses species barriers. Here, we compared the adaptation of an attenuated virus and its recombinant derivative to that of ferrets. Pre-existing mutations present at low frequency allowed faster adaptation of the non-recombinant virus compared to the recombinant virus. We identified a common point mutation in the nucleoprotein that affected the pathogenesis of both viruses. Our study shows that genetic memory facilitates environmental adaptation and that erasing this genetic memory by genetic engineering results in delayed and different adaptation to new environments, providing an important safety aspect for the generation of live-attenuated vaccines.

Genetic variations and gene expression profiles of Rice Black-streaked dwarf virus (RBSDV) in different host plants and insect vectors: insights from RNA-Seq analysis

Rice black-streaked dwarf virus (RBSDV) is an etiological agent of a destructive disease infecting some economically important crops from the Gramineae family in Asia. While RBSDV causes high yield losses, genetic characteristics of replicative viral populations have not been investigated within different host plants and insect vectors. Herein, eleven publicly available RNA-Seq datasets from Chinese RBSDV-infected rice, maize, and viruliferous planthopper (Laodelphax striatellus) were obtained from the NCBI database. The patterns of SNP and RNA expression profiles of expected RBSDV populations were analyzed by CLC Workbench 20 and Geneious Prime software. These analyses discovered 2,646 mutations with codon changes in RBSDV whole transcriptome and forty-seven co-mutated hotspots with high variant frequency within the crucial regions of S5-1, S5-2, S6, S7-1, S7-2, S9, and S10 open reading frames (ORFs) which are responsible for some virulence and host range functions. Moreover, three joint mutations are located on the three-dimensional protein of P9-1. The infected RBSDV-susceptible rice cultivar KTWYJ3 and indigenous planthopper datasets showed more co-mutated hotspot numbers than others. Our analyses showed the expression patterns of viral genomic fragments varied depending on the host type. Unlike planthopper, S5-1, S2, S6, and S9-1 ORFs, respectively had the greatest read numbers in host plants; and S5-2, S9-2, and S7-2 were expressed in the lowest level. These findings underscore virus/host complexes are effective in the genetic variations and gene expression profiles of plant viruses. Our analysis revealed no evidence of recombination events. Interestingly, the negative selection was observed at 12 RBSDV ORFs, except for position 1015 in the P1 protein, where a positive selection was detected. The research highlights the potential of SRA datasets for analysis of the virus cycle and enhances our understanding of RBSDV's genetic diversity and host specificity.

Recreating the biological steps of viral infection on a cell-free bioelectronic platform to profile viral variants of concern

Viral mutations frequently outpace technologies used to detect harmful variants. Given the continual emergence of SARS-CoV-2 variants, platforms that can identify the presence of a virus and its propensity for infection are needed. Our electronic biomembrane sensing platform recreates distinct SARS-CoV-2 host cell entry pathways and reports the progression of entry as electrical signals. We focus on two necessary entry processes mediated by the viral Spike protein: virus binding and membrane fusion, which can be distinguished electrically. We find that closely related variants of concern exhibit distinct fusion signatures that correlate with trends in cell-based infectivity assays, allowing us to report quantitative differences in their fusion characteristics and hence their infectivity potentials. We use SARS-CoV-2 as our prototype, but we anticipate that this platform can extend to other enveloped viruses and cell lines to quantifiably assess virus entry.

Characterization of a Molecular Clone of Deformed Wing Virus B

Honey bees (Apis mellifera) play a crucial role in agriculture through their pollination activities. However, they have faced significant health challenges over the past decades that can limit colony performance and even lead to collapse. A primary culprit is the parasitic mite Varroa destructor, known for transmitting harmful bee viruses. Among these viruses is deformed wing virus (DWV), which impacts bee pupae during their development, resulting in either pupal demise or in the emergence of crippled adult bees. In this study, we focused on DWV master variant B. DWV-B prevalence has risen sharply in recent decades and appears to be outcompeting variant A of DWV. We generated a molecular clone of a typical DWV-B strain to compare it with our established DWV-A clone, examining RNA replication, protein expression, and virulence. Initially, we analyzed the genome using RACE-PCR and RT-PCR techniques. Subsequently, we conducted full-genome RT-PCR and inserted the complete viral cDNA into a bacterial plasmid backbone. Phylogenetic comparisons with available full-length sequences were performed, followed by functional analyses using a live bee pupae model. Upon the transfection of in vitro-transcribed RNA, bee pupae exhibited symptoms of DWV infection, with detectable viral protein expression and stable RNA replication observed in subsequent virus passages. The DWV-B clone displayed a lower virulence compared to the DWV-A clone after the transfection of synthetic RNA, as evidenced by a reduced pupal mortality rate of only 20% compared to 80% in the case of DWV-A and a lack of malformations in 50% of the emerging bees. Comparable results were observed in experiments with low infection doses of the passaged virus clones. In these tests, 90% of bees infected with DWV-B showed no clinical symptoms, while 100% of pupae infected with DWV-A died. However, at high infection doses, both DWV-A and DWV-B caused mortality rates exceeding 90%. Taken together, we have generated an authentic virus clone of DWV-B and characterized it in animal experiments.

Accurate Recapitulation of Chikungunya Virus Complete Coding Sequence Phylogeny Using Variable Genome Regions for Genomic Surveillance

Chikungunya virus (CHIKV) is transmitted by mosquito bites and causes chikungunya fever (CHIKF). CHIKV has a single-stranded RNA genome and belongs to a single serotype with three genotypes. The Asian lineage has recently emerged in the Western Hemisphere, likely due to travel-associated introduction. Genetic variation accumulates in the CHIKV genome as the virus replicates, creating new lineages. Whole genome sequencing is ideal for studying virus evolution and spread but is expensive and complex. This study investigated whether specific, highly variable regions of the CHIKV genome could recapitulate the phylogeny obtained with a complete coding sequence (CDS). Our results revealed that concatenated highly variable regions accurately reconstructed CHIKV phylogeny, exhibiting statistically indistinguishable branch lengths and tree confidence compared to CDS. In addition, these regions adequately inferred the evolutionary relationships among CHIKV isolates from the American outbreak with similar results to the CDS. This finding suggests that highly variable regions can effectively capture the evolutionary relationships among CHIKV isolates, offering a simpler approach for future studies. This approach could be particularly valuable for large-scale surveillance efforts.

The Development of New Primer Sets for the Amplification and Sequencing of the Envelope Gene of All Dengue Virus Serotypes

Dengue virus (DENV) poses a significant threat to global health, infecting approximately 390 million people annually. This virus comprises four serotypes capable of causing severe disease. Genetic analyses are crucial for understanding the epidemiology, evolution, and spread of DENV. Although previous studies have focused on the envelope protein-coding (E) gene, only a few primers can efficiently detect and amplify the viral genes from multiple endemic countries simultaneously. In this study, we designed degenerate primer pairs for each DENV serotype to amplify and sequence the entire E gene, using globally representative sequences for each serotype. These primers were validated using DENV isolates from various Asian countries and demonstrated broad-spectrum detection capabilities and high-quality sequences. The primers provide effective tools for genetic analysis in the regions affected by dengue, aiding strain identification and epidemiological studies during outbreaks.

Dynamics of emergence and genetic diversity of dengue virus in Reunion Island from 2012 to 2022

BACKGROUND

Dengue is a major public health concern in Reunion Island, marked by recurrent epidemics, including successive outbreaks of dengue virus serotypes 1 and 2 (DENV1 and DENV2) with over 70,000 cases confirmed since 2017.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we used Oxford Nanopore NGS technology for sequencing virologically-confirmed samples and clinical isolates collected between 2012 and 2022 to investigate the molecular epidemiology and evolution of DENV in Reunion Island. Here, we generated and analyzed a total of 499 DENV1, 360 DENV2, and 18 DENV3 sequences. By phylogenetic analysis, we show that different genotypes and variants of DENV have circulated in the past decade that likely originated from Seychelles, Mayotte and Southeast Asia and highly affected areas in Asia and Africa.

CONCLUSIONS/SIGNIFICANCE

DENV sequences from Reunion Island exhibit a high genetic diversity which suggests regular introductions of new viral lineages from various Indian Ocean islands. The insights from our phylogenetic analysis may inform local health authorities about the endemicity of DENV variants circulating in Reunion Island and may improve dengue management and surveillance. This work emphasizes the importance of strong local coordination and collaboration to inform public health stakeholders in Reunion Island, neighboring areas, and mainland France.

Antigenic cooperation in viral populations: Transformation of functions of intra-host viral variants

In this paper, we study intra-host viral adaptation by antigenic cooperation - a mechanism of immune escape that serves as an alternative to the standard mechanism of escape by continuous genomic diversification and allows to explain a number of experimental observations associated with the establishment of chronic infections by highly mutable viruses. Within this mechanism, the topology of a cross-immunoreactivity network forces intra-host viral variants to specialize for complementary roles and adapt to the host's immune response as a quasi-social ecosystem. Here we study dynamical changes in immune adaptation caused by evolutionary and epidemiological events. First, we show that the emergence of a viral variant with altered antigenic features may result in a rapid re-arrangement of the viral ecosystem and a change in the roles played by existing viral variants. In particular, it may push the population under immune escape by genomic diversification towards the stable state of adaptation by antigenic cooperation. Next, we study the effect of a viral transmission between two chronically infected hosts, which results in the merging of two intra-host viral populations in the state of stable immune-adapted equilibrium. In this case, we also describe how the newly formed viral population adapts to the host's environment by changing the functions of its members. The results are obtained analytically for minimal cross-immunoreactivity networks and numerically for larger populations.

Whole Genome Sequencing of DENV-2 isolated from Aedes aegypti mosquitoes in Esmeraldas, Ecuador. Genomic epidemiology of genotype III Southern Asian-American in the country

Ecuador is a tropical country reporting Dengue virus (DENV) outbreaks with areas of hyperendemic viral transmission. Entomo-virological surveillance and monitoring effort conducted in the Northwestern border province of Esmeraldas in April 2022, five pools of female Aedes aegypti mosquitoes from a rural community tested positive for DENV serotype 2 by RT-qPCR. One pool was sequenced by Illumina MiSeq, and it corresponded to genotype III Southern Asian-American. Comparison with other genomes revealed genetic similarity to a human DENV genome sequenced in 2021, also from Esmeraldas. Potential introduction events to the country could have originated from Colombia, considering the vicinity of the collection sites to the neighboring country and high human movement. The inclusion of genomic information complements entomo-virological surveillance, providing valuable insights into genetic variants. This contribution enhances our understanding of Dengue virus (DENV) epidemiology in rural areas and guides evidence-based decisions for surveillance and interventions.

AltaiR: a C toolkit for alignment-free and temporal analysis of multi-FASTA data

BACKGROUND

Most viral genome sequences generated during the latest pandemic have presented new challenges for computational analysis. Analyzing millions of viral genomes in multi-FASTA format is computationally demanding, especially when using alignment-based methods. Most existing methods are not designed to handle such large datasets, often requiring the analysis to be divided into smaller parts to obtain results using available computational resources.

FINDINGS

We introduce AltaiR, a toolkit for analyzing multiple sequences in multi-FASTA format using exclusively alignment-free methodologies. AltaiR enables the identification of singularity and similarity patterns within sequences and computes static and temporal dynamics without restrictions on the number or size of input sequences. It automatically filters low-quality, biased, or deviant data. We demonstrate AltaiR's capabilities by analyzing more than 1.5 million full severe acute respiratory virus coronavirus 2 sequences, revealing interesting observations regarding viral genome characteristics over time, such as shifts in nucleotide composition, decreases in average Kolmogorov sequence complexity, and the evolution of the smallest sequences not found in the human host.

CONCLUSIONS

AltaiR can identify temporal characteristics and trends in large numbers of sequences, making it ideal for scenarios involving endemic or epidemic outbreaks with vast amounts of available sequence data. Implemented in C with multithreading and methodological optimizations, AltaiR is computationally efficient, flexible, and dependency-free. It accepts any sequence in FASTA format, including amino acid sequences. The complete toolkit is freely available at https://github.com/cobilab/altair.

Genetic Diversity of Grapevine Virus A in Three Australian Vineyards Using Amplicon High Throughput Sequencing (Amplicon-HTS)

Shiraz disease (SD) is one of the most destructive viral diseases of grapevines in Australia and is known to cause significant economic loss to local growers. Grapevine virus A (GVA) was reported to be the key pathogen associated with this disease. This study aimed to better understand the diversity of GVA variants both within and between individual SD and grapevine leafroll disease (LRD) affected grapevines located at vineyards in South Australia. Amplicon high throughput sequencing (Amplicon-HTS) combined with median-joining networks (MJNs) was used to analyze the variability in specific gene regions of GVA variants. Several GVAII variant groups contain samples from both vineyards studied, suggesting that these GVAII variants were from a common origin. Variant groups analyzed by MJNs using the overall data set denote that there may be a possible relationship between variant groups of GVA and the geographical location of the grapevines.

Evaluation of recombination detection methods for viral sequencing

Recombination is a key evolutionary driver in shaping novel viral populations and lineages. When unaccounted for, recombination can impact evolutionary estimations or complicate their interpretation. Therefore, identifying signals for recombination in sequencing data is a key prerequisite to further analyses. A repertoire of recombination detection methods (RDMs) have been developed over the past two decades; however, the prevalence of pandemic-scale viral sequencing data poses a computational challenge for existing methods. Here, we assessed eight RDMs: PhiPack (Profile), 3SEQ, GENECONV, recombination detection program (RDP) (OpenRDP), MaxChi (OpenRDP), Chimaera (OpenRDP), UCHIME (VSEARCH), and gmos; to determine if any are suitable for the analysis of bulk sequencing data. To test the performance and scalability of these methods, we analysed simulated viral sequencing data across a range of sequence diversities, recombination frequencies, and sample sizes. Furthermore, we provide a practical example for the analysis and validation of empirical data. We find that RDMs need to be scalable, use an analytical approach and resolution that is suitable for the intended research application, and are accurate for the properties of a given dataset (e.g. sequence diversity and estimated recombination frequency). Analysis of simulated and empirical data revealed that the assessed methods exhibited considerable trade-offs between these criteria. Overall, we provide general guidelines for the validation of recombination detection results, the benefits and shortcomings of each assessed method, and future considerations for recombination detection methods for the assessment of large-scale viral sequencing data.

Genomic stability of self-inactivating rabies

Transsynaptic viral vectors provide means to gain genetic access to neurons based on synaptic connectivity and are essential tools for the dissection of neural circuit function. Among them, the retrograde monosynaptic ΔG-Rabies has been widely used in neuroscience research. A recently developed engineered version of the ΔG-Rabies, the non-toxic self-inactivating (SiR) virus, allows the long term genetic manipulation of neural circuits. However, the high mutational rate of the rabies virus poses a risk that mutations targeting the key genetic regulatory element in the SiR genome could emerge and revert it to a canonical ΔG-Rabies. Such revertant mutations have recently been identified in a SiR batch. To address the origin, incidence and relevance of these mutations, we investigated the genomic stability of SiR in vitro and in vivo. We found that "revertant" mutations are rare and accumulate only when SiR is extensively amplified in vitro, particularly in suboptimal production cell lines that have insufficient levels of TEV protease activity. Moreover, we confirmed that SiR-CRE, unlike canonical ΔG-Rab-CRE or revertant-SiR-CRE, is non-toxic and that revertant mutations do not emerge in vivo during long-term experiments.

A Simulation Framework for Modeling the Within-Patient Evolutionary Dynamics of SARS-CoV-2

The global impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for rarely acting positive selection are best performed via comparison of empirical data with simulated data wherein commonly acting evolutionary factors, including mutation and recombination, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. Although there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intrahost evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them with existing empirical data. Of these, 592 models (∼5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intrahost SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed toward strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.

Evolution of the coronavirus spike protein in the full-length genome and defective viral genome under diverse selection pressures

How coronaviruses evolve by altering the structures of their full-length genome and defective viral genome (DVG) under dynamic selection pressures has not been studied. In this study, we aimed to experimentally identify the dynamic evolutionary patterns of the S protein sequence in the full-length genome and DVG under diverse selection pressures, including persistence, innate immunity and antiviral drugs. The evolutionary features of the S protein sequence in the full-length genome and in the DVG under diverse selection pressures are as follows: (i) the number of nucleotide (nt) mutations does not necessarily increase with the number of selection pressures; (ii) certain types of selection pressure(s) can lead to specific nt mutations; (iii) the mutated nt sequence can be reverted to the wild-type nt sequence under the certain type of selection pressure(s); (iv) the DVG can also undergo mutations and evolve independently of the full-length genome; and (v) DVG species are regulated during evolution under diverse selection pressures. The various evolutionary patterns of the S protein sequence in the full-length genome and DVG identified in this study may contribute to coronaviral fitness under diverse selection pressures.

Design principles and applications of synthetic self-replicating RNAs

With the advent of ever more sophisticated methods for the in vitro synthesis and the in vivo delivery of RNAs, synthetic mRNAs have gained substantial interest both for medical applications, as well as for biotechnology. However, in most biological systems exogeneous mRNAs possess only a limited half-life, especially in fast dividing cells. In contrast, viral RNAs can extend their lifetime by actively replicating inside their host. As such they may serve as scaffolds for the design of synthetic self-replicating RNAs (srRNA), which can be used to increase both the half-life and intracellular concentration of coding RNAs. Synthetic srRNAs may be used to enhance recombinant protein expression or induce the reprogramming of differentiated cells into pluripotent stem cells but also to create cell-free systems for research based on experimental evolution. In this article, we discuss the applications and design principles of srRNAs used for cellular reprogramming, mRNA-based vaccines and tools for synthetic biology. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA in Disease and Development > RNA in Development RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

A comprehensive study of SARS-CoV-2 main protease (Mpro) inhibitor-resistant mutants selected in a VSV-based system

Nirmatrelvir was the first protease inhibitor (PI) specifically developed against the SARS-CoV-2 main protease (3CLpro/Mpro) and licensed for clinical use. As SARS-CoV-2 continues to spread, variants resistant to nirmatrelvir and other currently available treatments are likely to arise. This study aimed to identify and characterize mutations that confer resistance to nirmatrelvir. To safely generate Mpro resistance mutations, we passaged a previously developed, chimeric vesicular stomatitis virus (VSV-Mpro) with increasing, yet suboptimal concentrations of nirmatrelvir. Using Wuhan-1 and Omicron Mpro variants, we selected a large set of mutants. Some mutations are frequently present in GISAID, suggesting their relevance in SARS-CoV-2. The resistance phenotype of a subset of mutations was characterized against clinically available PIs (nirmatrelvir and ensitrelvir) with cell-based and biochemical assays. Moreover, we showed the putative molecular mechanism of resistance based on in silico molecular modelling. These findings have implications on the development of future generation Mpro inhibitors, will help to understand SARS-CoV-2 protease-inhibitor-resistance mechanisms and show the relevance of specific mutations in the clinic, thereby informing treatment decisions.

Genetic variability of rice stripe virus after its pandemic in Jiangsu

BACKGROUND

Rice stripe virus (RSV) caused a serious disease pandemic in rice in East China between 2001 and 2010. The continuous integrated managements reduced virus epidemic year by year until it was non-epidemic. As an RNA virus, its genetic variability after undergoing a long-term non-epidemic period was meaningful to study. While in 2019, the sudden occurrence of RSV in Jiangsu provided an opportunity for the study.

METHODS AND RESULTS

The complete genome of JY2019, an RSV isolate from Jiangyan, was determined. A genotype profile of 22 isolates from China, Japan and Korea indicated that the isolates from Yunnan formed the subtype II, and other isolates clustered the subtype I. RNA 1-3 of JY2019 isolate well-clustered in the subtype I clade, and RNA 4 was also in subtype I, but it had a slight separation from other intra-group isolates. After phylogenetic analyses, it was considered NSvc4 gene contributed to the tendency, because it exhibited an obvious trend towards the subtype II (Yunnan) group. High sequence identity (100%) of NSvc4 between JY2019 and barnyardgrass isolate from different regions demonstrated genetic variation of NSvc4 was consistent in RSV natural populations in Jiangsu in the non-epidemic period. In the phylogenetic tree of all 74 NSvc4 genes, JY2019 belonged to a minor subtype Ib, suggesting the subtype Ib isolates might have existed in natural populations before the non-epidemic period, but not a dominant population.

CONCLUSIONS

Our results suggested that NSvc4 gene was susceptible to selection pressure, and the subtype Ib might be more adaptable for the interaction between RSV and hosts in the non-epidemic ecological conditions.

Deep spatial profiling of Venezuelan equine encephalitis virus reveals increased genetic diversity amidst neuroinflammation and cell death during brain infection

Venezuelan equine encephalitis virus (VEEV) causes a febrile illness that can progress to neurological disease with the possibility of death in human cases. The evaluation and optimization of therapeutics that target brain infections demands knowledge of the host's response to VEEV, the dynamics of infection, and the potential for within-host evolution of the virus. We hypothesized that selective pressures during infection of the brain may differ temporally and spatially and so we investigated the dynamics of the host response, viral transcript levels, and genetic variation of VEEV TC-83 in eight areas of the brain in mice over 7 days post-infection (dpi). Viral replication increased throughout the brain until 5-6 dpi and decreased thereafter with neurons as the main site of viral replication. Low levels of genetic diversity were noted on 1 dpi and were followed by an expansion in the genetic diversity of VEEV and nonsynonymous (Ns) mutations that peaked by 5 dpi. The pro-inflammatory response and the influx of immune cells mirrored the levels of virus and correlated with substantial damage to neurons by 5 dpi and increased activation of microglial cells and astrocytes. The prevalence and dynamics of Ns mutations suggest that the VEEV is under selection within the brain and that progressive neuroinflammation may play a role in acting as a selective pressure. IMPORTANCE Treatment of encephalitis in humans caused by Venezuelan equine encephalitis virus (VEEV) from natural or aerosol exposure is not available, and hence, there is a great interest to address this gap. In contrast to natural infections, therapeutic treatment of infections from aerosol exposure will require fast-acting drugs that rapidly penetrate the blood-brain barrier, engage sites of infection in the brain and mitigate the emergence of drug resistance. Therefore, it is important to understand not only VEEV pathogenesis, but the trafficking of the viral population within the brain, the potential for within-host evolution of the virus, and how VEEV might evolve resistance.

Maintenance of a host-specific minority mutation in the West Nile virus NS3

West Nile virus (WNV), the most prevalent arthropod-borne virus (arbovirus) in the United States, is maintained in a cycle between Culex spp. mosquitoes and birds. Arboviruses exist within hosts and vectors as a diverse set of closely related genotypes. In theory, this genetic diversity can facilitate adaptation to distinct environments during host cycling, yet host-specific fitness of minority genotypes has not been assessed. Utilizing WNV deep-sequencing data, we previously identified a naturally occurring, mosquito-biased substitution, NS3 P319L. Using both cell culture and experimental infection in natural hosts, we demonstrated that this substitution confers attenuation in vertebrate hosts and increased transmissibility by mosquitoes. Biochemical assays demonstrated temperature-sensitive ATPase activity consistent with host-specific phenotypes. Together these data confirm the maintenance of host-specific minority variants in arbovirus mutant swarms, suggest a unique role for NS3 in viral fitness, and demonstrate that intrahost sequence data can inform mechanisms of host-specific adaptation.

High-resolution mapping reveals the mechanism and contribution of genome insertions and deletions to RNA virus evolution

RNA viruses rapidly adapt to selective conditions due to the high intrinsic mutation rates of their RNA-dependent RNA polymerases (RdRps). Insertions and deletions (indels) in viral genomes are major contributors to both deleterious mutational load and evolutionary novelty, but remain understudied. To characterize the mechanistic details of their formation and evolutionary dynamics during infection, we developed a hybrid experimental-bioinformatic approach. This approach, called MultiMatch, extracts insertions and deletions from ultradeep sequencing experiments, including those occurring at extremely low frequencies, allowing us to map their genomic distribution and quantify the rates at which they occur. Mapping indel mutations in adapting poliovirus and dengue virus populations, we determine the rates of indel generation and identify mechanistic and functional constraints shaping indel diversity. Using poliovirus RdRp variants of distinct fidelity and genome recombination rates, we demonstrate tradeoffs between fidelity and Indel generation. Additionally, we show that maintaining translation frame and viral RNA structures constrain the Indel landscape and that, due to these significant fitness effects, Indels exert a significant deleterious load on adapting viral populations. Conversely, we uncover positively selected Indels that modulate RNA structure, generate protein variants, and produce defective interfering genomes in viral populations. Together, our analyses establish the kinetic and mechanistic tradeoffs between misincorporation, recombination, and Indel rates and reveal functional principles defining the central role of Indels in virus evolution, emergence, and the regulation of viral infection.

Molecular Population Genetics of Aspen Mosaic-Associated Virus in Finland and Sweden

Aspen mosaic-associated virus (AsMaV) is a newly identified Emaravirus, in the family Fimoviridae, Bunyavirales, associated with mosaic symptoms in aspen trees (Populus tremula). Aspen trees are widely distributed in Europe and understanding the population structure of AsMaV may aid in the development of better management strategies. The virus genome consists of five negative-sense single-stranded RNA (-ssRNA) molecules. To investigate the genetic diversity and population parameters of AsMaV, different regions of the genome were amplified and analyzed and full-length sequence of the divergent isolates were cloned and sequenced. The results show that RNA3 or nucleoprotein is a good representative for studying genetic diversity in AsMaV. Developed RT-PCR-RFLP was able to identify areas with a higher number of haplotypes and could be applied for screening the large number of samples. In general, AsMaV has a conserved genome and based on the phylogenetic studies, geographical structuring was observed in AsMaV isolates from Sweden and Finland, which could be attributed to founder effects. The genome of AsMaV is under purifying selection but not distributed uniformly on genomic RNAs. Distant AsMaV isolates displayed amino acid sequence variations compared to other isolates, and bioinformatic analysis predicted potential post-translational modification sites in some viral proteins.

A simulation framework for modeling the within-patient evolutionary dynamics of SARS-CoV-2

The global impact of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for positive selection are best performed via comparison of empirical data to simulated data wherein evolutionary factors, including mutation and recombination rates, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. While there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intra-host evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them to existing empirical data. Of these, 592 models (~5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intra-host SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed towards strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.

Flaviviruses and the Traveler: Around the World and to Your Stage. A Review of West Nile, Yellow Fever, Dengue, and Zika Viruses for the Practicing Pathologist

Flaviviruses are a genus of single-stranded RNA viruses that impose an important and growing burden to human health. There are over 3 billion individuals living in areas where flaviviruses are endemic. Flaviviruses and their arthropod vectors (which include mosquitoes and ticks) take advantage of global travel to expand their distribution and cause severe disease in humans, and they can be grouped according to their vector and pathogenicity. The mosquito-borne flaviviruses cause a spectrum of diseases from encephalitis to hepatitis and vascular shock syndrome, congenital abnormalities, and fetal death. Neurotropic infections such as Zika virus and West Nile virus cross the blood-brain barrier and infect neurons and other cells, leading to meningoencephalitis. In the hemorrhagic fever clade, there are yellow fever virus, the prototypical hemorrhagic fever virus that infects hepatocytes, and dengue virus, which infects cells of the reticuloendothelial system and can lead to a dramatic plasma cell leakage and shock syndrome. Zika virus also causes congenital infections and fetal death and is the first and only example of a teratogenic arbovirus in humans. Diagnostic testing for flaviviruses broadly includes the detection of viral RNA in serum (particularly within the first 10 days of symptoms), viral isolation by cell culture (rarely performed due to complexity and biosafety concerns), and histopathologic evaluation with immunohistochemistry and molecular testing on formalin-fixed paraffin-embedded tissue blocks. This review focuses on 4 mosquito-borne flaviviruses-West Nile, yellow fever, dengue, and Zika virus-and discusses the mechanisms of transmission, the role of travel in geographic distribution and epidemic emergence, and the clinical and histopathologic features of each. Finally, prevention strategies such as vector control and vaccination are discussed.

Novel approaches for the rapid development of rationally designed arbovirus vaccines

Vector-borne diseases, including those transmitted by mosquitoes, account for more than 17% of infectious diseases worldwide. This number is expected to rise with an increased spread of vector mosquitoes and viruses due to climate change and man-made alterations to ecosystems. Among the most common, medically relevant mosquito-borne infections are those caused by arthropod-borne viruses (arboviruses), especially members of the genera Flavivirus and Alphavirus. Arbovirus infections can cause severe disease in humans, livestock and wildlife. Severe consequences from infections include congenital malformations as well as arthritogenic, haemorrhagic or neuroinvasive disease. Inactivated or live-attenuated vaccines (LAVs) are available for a small number of arboviruses; however there are no licensed vaccines for the majority of these infections. Here we discuss recent developments in pan-arbovirus LAV approaches, from site-directed attenuation strategies targeting conserved determinants of virulence to universal strategies that utilize genome-wide re-coding of viral genomes. In addition to these approaches, we discuss novel strategies targeting mosquito saliva proteins that play an important role in virus transmission and pathogenesis in vertebrate hosts. For rapid pre-clinical evaluations of novel arbovirus vaccine candidates, representative in vitro and in vivo experimental systems are required to assess the desired specific immune responses. Here we discuss promising models to study attenuation of neuroinvasion, neurovirulence and virus transmission, as well as antibody induction and potential for cross-reactivity. Investigating broadly applicable vaccination strategies to target the direct interface of the vertebrate host, the mosquito vector and the viral pathogen is a prime example of a One Health strategy to tackle human and animal diseases.

Non-equilibrium virus particle dynamics: Microsecond MD simulations of the complete Flock House virus capsid under different conditions

Flock House virus (FHV) is an animal virus and considered a model system for non-enveloped viruses. It has a small, icosahedral capsid (T=3) and a bipartite positive-sense RNA genome. We present an extensive study of the FHV capsid dynamics from all-atom molecular dynamics simulations of the complete capsid. The simulations explore different biologically relevant conditions (neutral/low pH, with/without RNA in the capsid) using the CHARMM force field. The results show that low pH destabilizes the capsid, causing radial expansion, and RNA stabilizes the capsid. The finding of low pH destabilization is biologically relevant because the capsid is exposed to low pH in the endosome, where conformational changes occur leading to genome release. We also observe structural changes at the fivefold and twofold symmetry axes that likely relate to the externalization of membrane active γ peptides through the fivefold vertex and extrusion of RNA at the twofold axis. Simulations using the Amber force field at neutral pH are also performed and display similar characteristics to the CHARMM simulations.

Iron Saturation Drives Lactoferrin Effects on Oxidative Stress and Neurotoxicity Induced by HIV-1 Tat

The Trans-Activator of Transcription (Tat) of Human Immunodeficiency Virus (HIV-1) is involved in virus replication and infection and can promote oxidative stress in human astroglial cells. In response, host cells activate transcription of antioxidant genes, including a subunit of System Xc- cystine/glutamate antiporter which, in turn, can trigger glutamate-mediated excitotoxicity. Here, we present data on the efficacy of bovine Lactoferrin (bLf), both in its native (Nat-bLf) and iron-saturated (Holo-bLf) forms, in counteracting oxidative stress in U373 human astroglial cells constitutively expressing the viral protein (U373-Tat). Our results show that, dependent on iron saturation, both Nat-bLf and Holo-bLf can boost host antioxidant response by up-regulating System Xc- and the cell iron exporter Ferroportin via the Nuclear factor erythroid 2-related factor (Nrf2) pathway, thus reducing Reactive Oxygen Species (ROS)-mediated lipid peroxidation and DNA damage in astrocytes. In U373-Tat cells, both forms of bLf restore the physiological internalization of Transferrin (Tf) Receptor 1, the molecular gate for Tf-bound iron uptake. The involvement of astrocytic antioxidant response in Tat-mediated neurotoxicity was evaluated in co-cultures of U373-Tat with human neuronal SH-SY5Y cells. The results show that the Holo-bLf exacerbates Tat-induced excitotoxicity on SH-SY5Y, which is directly dependent on System-Xc- upregulation, thus highlighting the mechanistic role of iron in the biological activities of the glycoprotein.

Viruses of a key coral symbiont exhibit temperature-driven productivity across a reefscape

Viruses can affect coral health by infecting their symbiotic dinoflagellate partners (Symbiodiniaceae). Yet, viral dynamics in coral colonies exposed to environmental stress have not been studied at the reef scale, particularly within individual viral lineages. We sequenced the viral major capsid protein (mcp) gene of positive-sense single-stranded RNA viruses known to infect symbiotic dinoflagellates ('dinoRNAVs') to analyze their dynamics in the reef-building coral, Porites lobata. We repeatedly sampled 54 colonies harboring Cladocopium C15 dinoflagellates, across three environmentally distinct reef zones (fringing reef, back reef, and forereef) around the island of Moorea, French Polynesia over a 3-year period and spanning a reef-wide thermal stress event. By the end of the sampling period, 28% (5/18) of corals in the fringing reef experienced partial mortality versus 78% (14/18) of corals in the forereef. Over 90% (50/54) of colonies had detectable dinoRNAV infections. Reef zone influenced the composition and richness of viral mcp amino acid types ('aminotypes'), with the fringing reef containing the highest aminotype richness. The reef-wide thermal stress event significantly increased aminotype dispersion, and this pattern was strongest in the colonies that experienced partial mortality. These findings demonstrate that dinoRNAV infections respond to environmental fluctuations experienced in situ on reefs. Further, viral productivity will likely increase as ocean temperatures continue to rise, potentially impacting the foundational symbiosis underpinning coral reef ecosystems.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

Bioinformatics Analysis of Mutations Sheds Light on the Evolution of Dengue NS1 Protein With Implications in the Identification of Potential Functional and Druggable Sites

Non-structural protein (NS1) is a 350 amino acid long conserved protein in the dengue virus. Conservation of NS1 is expected due to its importance in dengue pathogenesis. The protein is known to exist in dimeric and hexameric states. The dimeric state is involved in its interaction with host proteins and viral replication, and the hexameric state is involved in viral invasion. In this work, we performed extensive structure and sequence analysis of NS1 protein, and uncovered the role of NS1 quaternary states in its evolution. A three-dimensional modeling of unresolved loop regions in NS1 structure is performed. "Conserved" and "Variable" regions within NS1 protein were identified from sequences obtained from patient samples and the role of compensatory mutations in selecting destabilizing mutations were identified. Molecular dynamics (MD) simulations were performed to extensively study the effect of a few mutations on NS1 structure stability and compensatory mutations. Virtual saturation mutagenesis, predicting the effect of every individual amino acid substitution on NS1 stability sequentially, revealed virtual-conserved and variable sites. The increase in number of observed and virtual-conserved regions across NS1 quaternary states suggest the role of higher order structure formation in its evolutionary conservation. Our sequence and structure analysis could enable in identifying possible protein-protein interfaces and druggable sites. Virtual screening of nearly 10,000 small molecules, including FDA-approved drugs, permitted us to recognize six drug-like molecules targeting the dimeric sites. These molecules could be promising due to their stable interactions with NS1 throughout the simulation.

Brome mosaic virus detected in Kansas wheat co-infected with other common wheat viruses

Wheat breeders are developing new virus-resistant varieties; however, it is assumed that only a few viruses or well-known viruses are present in the field. New sequencing technology is allowing for better determination of natural field virus populations. For three years, 2019-2021, Kansas wheat field surveys were conducted to determine the constituents of natural field virus populations using nanopore sequencing. During analysis, brome mosaic virus (BMV) was identified for the first time in Kansas but was in association with other wheat viruses. Brome mosaic virus was identified from 29 out of 47 different Kansas counties sampled and 44% of the total samples. BMV was found co-infected with wheat streak mosaic virus (WSMV) and Triticum mosaic virus (TriMV) in 27.8% of the samples, with WSMV only (13.9%) and co-infected with WSMV + TriMV + High Plains wheat mosaic emaravirus (HPWMoV) (13.9%). RNA genomes of Kansas BMV isolates had 99.4 to 100% nucleotide and amino acid sequence identity, respectively, to each other. RNA2a possessed relatively high divergence (π = 0.01) compared to RNA1a and RNA3a (π = 0.004). Coding regions of all BMV RNAs were considered negative for purifying selection pressure as nonsynonymous and synonymous nucleotide ratio was less than one (dNs/dS >1). The identification of BMV in Kansas virus populations adds another layer of complexity to plant breeding. This work provides information to improve tools to aid in monitoring, detecting, and determining the variation within BMV.

Fidelity of Ribonucleotide Incorporation by the SARS-CoV-2 Replication Complex

The SARS-CoV-2 coronavirus has caused a global pandemic. Despite the initial success of vaccines at preventing infection, genomic variation has led to the proliferation of variants capable of higher infectivity. Mutations in the SARS-CoV-2 genome are the consequence of replication errors, highlighting the importance of understanding the determinants of SARS-CoV-2 replication fidelity. The RNA-dependent RNA polymerase (RdRp) is the central catalytic subunit for SARS-CoV-2 RNA replication and genome transcription. Here, we report the fidelity of ribonucleotide incorporation by SARS-CoV-2 RdRp (nsp12), along with its co-factors nsp7/nsp8, using steady-state kinetic analysis. Our analysis suggests that in the absence of the proofreading subunit (nsp14), the nsp12/7/8 complex has a surprisingly low base substitution fidelity (10-1-10-3). This is orders of magnitude lower than the fidelity reported for other coronaviruses (10-6-10-7), highlighting the importance of proofreading for faithful SARS-CoV-2 replication. We performed a mutational analysis of all reported SARS-CoV-2 genomes and identified mutations in both nsp12 and nsp14 that appear likely to lower viral replication fidelity through mechanisms that include impairing the nsp14 exonuclease activity or its association with the RdRp. Our observations provide novel insight into the mechanistic basis of replication fidelity in SARS-CoV-2 and the potential effect of nsp12 and nsp14 mutations on replication fidelity, informing the development of future antiviral agents and SARS-CoV-2 vaccines.