Viral quasispecies evolution

The prevalence of low-level viraemia and its association with virological failure in people living with HIV: a systematic review and meta-analysis

Low-level viraemia (LLV) following antiretroviral therapy (ART) in people living with HIV (PLWH) has not received sufficient attention. To the determine the prevalence of LLV and its association with virological failure (VF), we systematically reviewed evidence-based interventions for PLWH. We searched PubMed, the Cochrane Library, Embase, and Web of Science from inception to 22 May 2024. Cohorts with samples sizes smaller than 1000 in size were excluded. Data from 16 cohort studies, encompassing 13,49,306 PLWH, revealed a pooled prevalence of LLV of 13.81%. Relative risk (RR) and 95% confidence intervals (CI) identified the following risk factors for LLV: viral load (VL) ≥ 105 copies/mL at baseline (1.79, 1.11-2.88), AIDS-defined illness at baseline (1.24, 1.10-1.40), and protease inhibitor-based regimen at ART initiation (1.53, 1.45-1.62) are the risk factors for LLV. Conversely, CD4 count ≥200 cells/μL at baseline (0.90, 0.82-0.98), non-nucleoside reverse transcriptase inhibitor-based regimen (0.81, 0.68-0.96) and the integrase strand transfer inhibitor (INSTI)-based regimen (0.60, 0.42-0.85) were associated with a reduced risk of LLV. Pooling the adjusted hazard ratio (aHR) and the 95% CI, we found that LLV increased the risk of VF with rising VL among 96,711 PLWH (aHR 2.77, 95% CI 2.03-3.76) and increased the risk of all-cause mortality at high VL levels among 14,229 PLWH (aHR 1.66, 95% CI 1.16-2.37). Therefore, the prevalence of LLV in PLWH should not be overlooked. This study aims to guide better management strategies to improve clinical outcomes in patients with LLV.

A method for in-depth analysis of circular DNA virus populations by unambiguously profiling the low abundant virus variants and partial genomic components

Severe epidemic outbreaks of diseases associated with newly emerging strains of single-stranded DNA (ssDNA) viruses have led to serious economic losses of numerous important food crops. While the current mitigation strategies are mostly relying on the deployment of genetic resistance in crop varieties, the constantly evolving virus populations have the potential to rapidly break virus resistance. Therefore, the development of diagnostic tools enabling early detection of virus variants associated with hypervirulence and/or expansion to new host species is urgently needed as an effective mitigation solution. Here, we introduce a novel approach by designing a pipeline that allows accurately identifying and characterizing the full-length sequence variants of viral circular DNA genomes utilizing Nanopore sequencing technology and the bioinformatics tool Genome Detective. We demonstrate that the pipeline is suitable to provide an accurate and in-depth analysis of monopartite Tomato yellow leaf curl Sardinia virus (TYLCSV) and multipartite Banana bunchy top virus (BBTV) ssDNA virus populations resulting in the profiling of high- and low-frequency virus variants with ≥1% relative abundance. The approach also enabled the unambiguous detection and characterization of four TYLCSV partial genomic sequences as well as several partial genomic sequences for each BBTV genomic component not previously reported and accumulating during infection.

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.

Human Riboviruses: A Comprehensive Study

The urgency to understand the complex interactions between viruses, their animal reservoirs, and human populations has been necessitated by the continuous spread of zoonotic viral diseases as evidenced in epidemics and pandemics throughout human history. Riboviruses are involved in some of the most prevalent human diseases, responsible for causing epidemics and pandemics. These viruses have an animal origin and have been known to cross the inter-species barrier time and time again, eventually infecting human beings. Their evolution has been a long road to harbour important adaptations for increasing fitness, mutability and virulence; a result of natural selection and mutation pressure, making these viruses highly infectious and difficult to counter. Accumulating favourable mutations in the course, they imitate the GC content and codon usage patterns of the host for maximising the chances of infection. A myriad of viral and host factors determine the fate of specific viral infections, which may include virus protein and host receptor compatibility, host restriction factors and others. Thus, understanding the biology, transmission and molecular mechanisms of Riboviruses is essential for the development of effective antiviral treatments, vaccine development and strategies to prevent and control viral infections. Keeping these aspects in mind, this review aims to provide a holistic approach towards understanding Riboviruses.

Deciphering mixed infections by plant RNA virus and reconstructing complete genomes simultaneously present within-host

Local co-circulation of multiple phylogenetic lineages is particularly likely for rapidly evolving pathogens in the current context of globalisation. When different phylogenetic lineages co-occur in the same fields, they may be simultaneously present in the same host plant (i.e. mixed infection), with potentially important consequences for disease outcome. This is the case in Burkina Faso for the rice yellow mottle virus (RYMV), which is endemic to Africa and a major constraint on rice production. We aimed to decipher the distinct RYMV isolates that simultaneously infect a single rice plant and to sequence their genomes. To this end, we tested different sequencing strategies, and we finally combined direct cDNA ONT (Oxford Nanopore Technology) sequencing with the bioinformatics tool RVhaplo. This method was validated by the successful reconstruction of two viral genomes that were less than a hundred nucleotides apart (out of a genome of 4450nt length, i.e. 2-3%), and present in artificial mixes at a ratio of up to a 99/1. We then used this method to subsequently analyze mixed infections from field samples, revealing up to three RYMV isolates within one single rice plant sample from Burkina Faso. In most cases, the complete genome sequences were obtained, which is particularly important for a better estimation of viral diversity and the detection of recombination events. The method described thus allows to identify various haplotypes of RYMV simultaneously infecting a single rice plant, obtaining their full-length sequences, as well as a rough estimate of relative frequencies within the sample. It is efficient, cost-effective, as well as portable, so that it could further be implemented where RYMV is endemic. Prospects include unravelling mixed infections with other RNA viruses that threaten crop production worldwide.

Human enteroviruses and the long road to acute flacid paralysis eradication

Enteroviruses (EVs) are a highly diverse group of viruses multiplying primarily in the gastrointestinal tract and/or the upper respiratory tract, initially distributed in two separate genera: Enterovirus and Rhinovirus, respectively. According to the similarities in genome organization and particle structure, rhinovirus species were later reclassified as also belonging to genus Enterovirus. Human EV infections are usually asymptomatic or causing mild clinical manifestations. Nevertheless, some EV infections may derive in severe neural complications, including acute flaccid paralysis (AFP) such as poliomyelitis, whose etiological agent is poliovirus, a member of the Enterovirus C species. The inactivated polio vaccine (IPV) and particularly the oral attenuated polio vaccine (OPV) have contributed to the virtual eradication of the disease. However, sustained global circulation of vaccine-derived poliovirus 2 (cVDPV2), originated from the genetic instability of OPV strain 2 and intertypic recombination between Sabin OPV strains and members of the Enterovirus C species, still causes outbreaks of AFP worldwide. In addition, humanitarian crises, in particular armed conflicts, hamper polio vaccination campaigns and facilitate the occurrence of cases. Additionally, besides poliovirus, other EV may also cause AFP, among them EV A71 or EV D68, and it is highly advisable to implement wastewater surveillance to elucidate the occurrence of not only polioviruses, but also of other EV susceptible to derive in serious neural complications, since the screening of viral RNA in cerebrospinal fluid samples in patients suffering from AFP is not a reliable diagnostic tool.

Tracking cryptic SARS-CoV-2 hospital outbreak through quasispecies analysis

BACKGROUND

Since the beginning of the pandemic, contact tracing has been one of the most relevant issues to understand SARS-CoV-2 transmission dynamics and, in this context, the analysis of quasispecies may turn out to be a useful tool for outbreak investigations. Analysis of the intra-host single nucleotide variants (iSNVs) found in the nsp2, ORF3, and ORF7 genes of SARS-CoV-2 was conducted in order to correctly identify virus transmission chain among patients hospitalized in Brescia Civic Hospital.

METHODS

During the period between August and October 2023, 13 nasopharyngeal specimens, collected from patients admitted to Brescia Civic Hospital, were tested for SARS-CoV-2 positivity and molecularly characterized. Firstly, a phylogenetic analysis was performed to evaluate if they were epidemiologically linked and, then, the Beta-binomial method was used to estimate the transmission bottleneck size (Nb) and quantify the number of viral particles transmitted from one individual (donor) to another (recipient).

RESULTS

According to the molecular characterization of specimens, we identified two transmission clusters in the cardiology unit: the first cluster concerned patients tested positive for the HV.1/EG.5.1.6 lineage, while the second cluster concerned patients tested positive for the FL.10.1 lineage. Moreover, evaluating the bottleneck size, we were able to solve SARS-CoV-2 transmission chain among infected patients.

CONCLUSION

Our method shows that it is possible to conduct a tracing study using a genomic approach based on iSNVs analysis.

Optimized Directed Virus Evolution to Accelerate the Generation of Oncolytic Coxsackievirus B3 Adapted to Resistant Colorectal Cancer Cells

Recently, we demonstrated that the oncolytic Coxsackievirus B3 (CVB3) strain PD-H can be efficiently adapted to resistant colorectal cancer cells through dose-dependent passaging in colorectal cancer cells. However, the method is time-consuming, which limits its clinical applicability. Here, we investigated whether the manufacturing time of the adapted virus can be reduced by replacing the dose-based passaging with volume-based passaging. For this purpose, the murine colorectal carcinoma cell line MC38, resistant to PD-H-induced lysis, was initially infected with PD-H at 0.1 multiplicity of infection (MOI). For subsequent passages, 15-30 µL of a 1:10 dilution of the cell culture supernatant was transferred to fresh MC38 cells early after virus-induced cell lysis became visible. By virus passage 10, complete cell lysis of MC38 cells was achieved. Sequencing of the passage 10 virus (P-10) revealed two nucleotide substitutions in the 5' UTR and six amino acid changes in the viral polyprotein compared to the PD-H founder. P-10, however, consisted of a heterogeneous virus population. Therefore, the detected mutations were introduced into the cDNA of PD-H, from which the recombinant virus PD-MC38 was generated. PD-MC38 exhibited significantly enhanced replication and lytic activity in MC38 cells compared to PD-H, whereas its oncolytic activity in other colorectal cancer cell lines was comparable to or even lower than that of PD-H. These findings demonstrate that volume-based passaging is suitable to generate tumor cell-specific adapted PD-H. Moreover, compared to the dose-dependent passaging, volume-based passaging significantly reduced the time required to generate the adapted virus.

Dynamic expedition of leading mutations in SARS-CoV-2 spike glycoproteins

The continuous evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the recent pandemic, has generated countless new variants with varying fitness. Mutations of the spike glycoprotein play a particularly vital role in shaping its evolutionary trajectory, as they have the capability to alter its infectivity and antigenicity. We present a time-resolved statistical method, Dynamic Expedition of Leading Mutations (deLemus), to analyze the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein. The proposed L -index of the deLemus method is effective in quantifying the mutation strength of each amino acid site and outlining evolutionarily significant sites, allowing the comprehensive characterization of the evolutionary mutation pattern of the spike glycoprotein.

DiMA: sequence diversity dynamics analyser for viruses

Sequence diversity is one of the major challenges in the design of diagnostic, prophylactic, and therapeutic interventions against viruses. DiMA is a novel tool that is big data-ready and designed to facilitate the dissection of sequence diversity dynamics for viruses. DiMA stands out from other diversity analysis tools by offering various unique features. DiMA provides a quantitative overview of sequence (DNA/RNA/protein) diversity by use of Shannon's entropy corrected for size bias, applied via a user-defined k-mer sliding window to an input alignment file, and each k-mer position is dissected to various diversity motifs. The motifs are defined based on the probability of distinct sequences at a given k-mer alignment position, whereby an index is the predominant sequence, while all the others are (total) variants to the index. The total variants are sub-classified into the major (most common) variant, minor variants (occurring more than once and of incidence lower than the major), and the unique (singleton) variants. DiMA allows user-defined, sequence metadata enrichment for analyses of the motifs. The application of DiMA was demonstrated for the alignment data of the relatively conserved Spike protein (2,106,985 sequences) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the relatively highly diverse pol gene (2637) of the human immunodeficiency virus-1 (HIV-1). The tool is publicly available as a web server (https://dima.bezmialem.edu.tr), as a Python library (via PyPi) and as a command line client (via GitHub).

A measles virus collective infectious unit that caused lethal human brain disease includes many locally restricted and few widespread copy-back defective genomes

During virus replication in cultured cells, copy-back defective viral genomes (cbDVGs) can arise. CbDVGs are powerful inducers of innate immune responses in vitro, but their occurrence and impact on natural infections of human hosts remain poorly defined. We asked whether cbDVGs were generated in the brain of a patient who succumbed to subacute sclerosing panencephalitis (SSPE) about 20 years after acute measles virus (MeV) infection. Previous analyses of 13 brain specimens of this patient indicated that a collective infectious unit (CIU) drove lethal MeV spread. In this study, we identified 276 replication-competent cbDVG species, each present in over 100 copies in the brain. Six species were detected in multiple forebrain locations, implying that they travelled long-distance with the CIU. The cbDVG to full-length genomes ratio was often close to 1 (0.6-1.74). Most cbDVGs were 324-2,000 bases in length, corresponding to 2%-12% of the full-length genome; all are predicted to have complementary terminal sequences. If improperly encapsidated, these sequences have the potential to form double-stranded structures that can induce innate immune responses. To assess this, we examined the transcriptome of all brain specimens. Several interferon and inflammatory response genes were upregulated, but upregulation levels did not correlate with cbDVG levels in the specimens. Thus, the CIU that drove MeV pathogenesis in this brain includes, in addition to two complementary full-length genome populations, many locally restricted and few widespread cbDVG species. The widespread cbDVG species may have been positively selected but how they impacted pathogenesis remains to be determined.IMPORTANCECopy-back defective viral genomes (cbDVGs) can drive virus-host interactions. They can suppress virus replication directly, by competing with full-length genomes, or indirectly by stimulating antiviral immunity. In vitro, cbDVG can slow down infections and promote persistence, but there is limited documentation of their presence in human hosts or of their impact on disease. We had the unique opportunity to analyze the brain of a patient who succumbed to subacute sclerosing panencephalitis, a rare but lethal consequence of measles. We detected more than 270 distinct cbDVG species; most were restricted to one specimen, but several reached all lobes of the forebrain, suggesting positive selection. Our analyses provide the missing knowledge of the diversity of cbDVG in a natural infection of a human host. They also reveal that a collective infectious unit that caused lethal human brain disease includes few widespread cbDVG, in addition to two ubiquitous complementary full-length genome populations.

The impact of immune cell interactions on virus quasi-species formation

The process of viral infection spreading in tissues was influenced by various factors, including virus replication within host cells, transportation, and the immune response. Reaction-diffusion systems provided a suitable framework for examining this process. In this work, we studied a nonlocal reaction-diffusion system of equations that modeled the distribution of viruses based on their genotypes and their interaction with the immune response. It was shown that the infection may persist at a certain level alongside a chronic immune response, exhibiting spatially uniform or oscillatory behavior. Finally, the immune cells may become entirely depleted, leading to a high viral load persisting in the tissue. Numerical simulations were employed to elucidate the nonlinear dynamics and pattern formation inherent in the nonlocal model.

Error thresholds in the presence of epistatic interactions

Models for viral populations with high replication error rates (such as RNA viruses) rely on the quasispecies concept, in which mutational pressure beyond the so-called "error threshold" leads to a loss of essential genetic information and population collapse, an effect known as the "error catastrophe." We explain how crossing this threshold, as a result of increasing mutation rates, can be understood as a second-order phase transition, even in the presence of lethal mutations. In particular, we show that, in fitness landscapes with a single peak, this collapse is equivalent to a ferroparamagnetic transition, where the back-mutation rate plays the role of the external magnetic field. We then generalize this framework to rugged fitness landscapes, like the ones that arise from epistatic interactions, and provide numerical evidence that there is a transition from a high average fitness regime to a low average fitness one, similarly to single-peaked landscapes. The onset of the transition is heralded by a sudden change in the susceptibility to variations in the mutation rate. We use insight from replica symmetry breaking mechanisms in spin glasses, in particular by considering the fluctuations of the genotype similarity distribution as the order parameter.

Discovery of novel and known viruses associated with toxigenic and non-toxigenic bloom forming diatoms from the Northern Adriatic Sea

Algal blooms impact trophic interactions, community structure and element fluxes. Despite playing an important role in the demise of phytoplankton blooms, only few viruses infecting diatoms have been cultured. Pseudo-nitzschia is a widespread diatom genus that commonly blooms in coastal waters and contains toxin-producing species. This study describes the characterization of a novel virus infecting the toxigenic species Pseudo-nitzschia galaxiae isolated from the northern Adriatic Sea. The ssRNA virus PnGalRNAV has 29.5 nm ± 1.2 nm icosahedral virions and a genome size of 8.8 kb. It belongs to the picorna-like Marnaviridae family and shows high specificity for P. galaxiae infecting two genetically and morphologically distinct strains. We found two genetically distinct types of this virus and screening of the global virome database revealed matching sequences from the Mediterranean region and China, suggesting its global distribution. Another virus of similar shape and size infecting Pseudo-nitzschia calliantha was found, but its genome could not be determined. In addition, we have obtained and characterized a new virus that infects Chaetoceros tenuissimus. The replicase protein of this virus is very similar to the previously described ChTenDNAV type-II virus, but it has a unique genome and infection pattern. Our study is an important contribution to the collective diatom virus culture collection and will allow further investigation into how these viruses control diatom bloom termination, carbon export and toxin release in the case of Pseudo-nitzschia.

Chronic infections can generate SARS-CoV-2-like bursts of viral evolution without epistasis

Multiple SARS-CoV-2 variants have arisen during the first years of the pandemic, often bearing many new mutations. Several explanations have been offered for the surprisingly sudden emergence of multiple mutations that enhance viral fitness, including cryptic transmission, spillover from animal reservoirs, epistasis between mutations, and chronic infections. Here, we simulated pathogen evolution combining within-host replication and between-host transmission. We found that, under certain conditions, chronic infections can lead to SARS-CoV-2-like bursts of mutations even without epistasis. Chronic infections can also increase the global evolutionary rate of a pathogen even in the absence of clear mutational bursts. Overall, our study supports chronic infections as a plausible origin for highly mutated SARS-CoV-2 variants. More generally, we also describe how chronic infections can influence pathogen evolution under different scenarios.

Enhanced production of recombinant coxsackievirus A16 using a serum-free HEK293A suspension culture system for bivalent enterovirus vaccine development

Coxsackievirus A16 (CVA16) is one of the primary pathogens that causes hand, foot, and mouth disease (HFMD) in young children. In previous studies, CVA16 vaccine development has encountered several challenges, such as inefficient replication of the CVA16 virus in present culture systems, the induction of only mild neutralizing antibody titers, and neutralizing antibodies induced by certain vaccine candidates that are unable to protect against CVA16 viral challenge. In this study, we constructed a DNA-launched CVA16 infectious clone (CVA16ic) based on the genomic sequence of the CVA16 N5079 strain to minimize interference from viral quasispecies. The biochemical properties of this CVA16ic strain were similar to those of its parental strain. Serum-free HEK293A suspension cells, which produced higher virus titers than Vero cells, were demonstrated to improve CVA16 production yields. In addition, our study showed that inactivated EV-A71 antigens could enhance the immunogenicity of inactivated CVA16 mature/full particles (F-particles), suggesting that a bivalent CVA16 and EV-A71 vaccine may be an effective strategy for CVA16 vaccine development. These findings are expected to provide novel strategies and accelerate the development of bivalent HFMD vaccines.

Nanopore sequencing provides snapshots of the genetic variation within salmonid alphavirus-3 (SAV3) during an ongoing infection in Atlantic salmon (Salmo salar) and brown trout (Salmo trutta)

Frequent RNA virus mutations raise concerns about evolving virulent variants. The purpose of this study was to investigate genetic variation in salmonid alphavirus-3 (SAV3) over the course of an experimental infection in Atlantic salmon and brown trout. Atlantic salmon and brown trout parr were infected using a cohabitation challenge, and heart samples were collected for analysis of the SAV3 genome at 2-, 4- and 8-weeks post-challenge. PCR was used to amplify eight overlapping amplicons covering 98.8% of the SAV3 genome. The amplicons were subsequently sequenced using the Nanopore platform. Nanopore sequencing identified a multitude of single nucleotide variants (SNVs) and deletions. The variation was widespread across the SAV3 genome in samples from both species. Mostly, specific SNVs were observed in single fish at some sampling time points, but two relatively frequent (i.e., major) SNVs were observed in two out of four fish within the same experimental group. Two other, less frequent (i.e., minor) SNVs only showed an increase in frequency in brown trout. Nanopore reads were de novo clustered using a 99% sequence identity threshold. For each amplicon, a number of variant clusters were observed that were defined by relatively large deletions. Nonmetric multidimensional scaling analysis integrating the cluster data for eight amplicons indicated that late in infection, SAV3 genomes isolated from brown trout had greater variation than those from Atlantic salmon. The sequencing methods and bioinformatics pipeline presented in this study provide an approach to investigate the composition of genetic diversity during viral infections.

Phylodynamic and Epistatic Analysis of Coxsackievirus A24 and Its Variant

Coxsackievirus A24 (CV-A24) is a human enterovirus that causes acute flaccid paralysis. However, a Coxsackievirus A24 variant (CV-A24v) is the most common cause of eye infections. The causes of these variable pathogenicity and tissue tropism remain unclear. To elucidate the phylodynamics of CV-A24 and CV-A24v, we analyzed a dataset of 66 strains using Bayesian phylodynamic approach, along with detailed sequence variation and epistatic analyses. Six CV-A24 strains available in GenBank and 60 CV-A24v strains, including 11 Taiwanese strains, were included in this study. The results revealed striking differences between CV-A24 and CV-A24v exhibiting long terminal branches in the phylogenetic tree, respectively. CV-A24v presented distinct ladder-like clustering, indicating immune escape mechanisms. Notably, 10 genetic recombination events in the 3D regions were identified. Furthermore, 11 missense mutation signatures were detected to differentiate CV-A24 and CV-A24v; among these mutations, the F810Y substitution may significantly affect the secondary structure of the GH loop of VP1 and subsequently affect the epitopes of the capsid proteins. In conclusion, this study provides critical insights into the evolutionary dynamics and epidemiological characteristics of CV-A24 and CV-A24v, and highlights the differences in viral evolution and tissue tropism.

Adapting the rhizome concept to an extended definition of viral quasispecies and the implications for molecular evolution

The rhizome concept proposed by Gilles Deleuze and Félix Guattari offers a novel perspective on the organization and interdependence of complex constellations of heterogeneous entities, their mapping and their ruptures. The emphasis of the present study is placed on the dynamics of contacts and communication among such entities that arise from experimentation, without any favored hierarchy or origin. When applied to biological evolution, the rhizome concept integrates all types of heterogeneity resulting from "symbiotic" relationships among living beings (or their genomic material), horizontal genetic transfer, recombination and mutation, and breaks away from the approach that gives rise to the phylogenetic tree of life. It has already been applied to describe the dynamics and evolution of RNA viruses. Thus, here we introduce a novel framework for the interpretation the viral quasispecies concept, which explains the evolution of RNA virus populations as the result of dynamic interconnections and multifaceted interdependence between highly heterogeneous viral sequences and its inherently heterogeneous host cells. The rhizome network perspective underlines even further the medical implications of the broad mutant spectra of viruses that are in constant flow, given the multiple pathways they have available for fitness loss and gain.

Selective and non-selective evolutionary signatures found in the simplest replicative biological entities

Mitoviruses, which are considered evolutionary relics of extinct alpha-proteobacteria RNA phages, represent one of the simplest self-replicating biological systems. This study aims to quantitatively describe genomes and identify potential genomic signatures that support the protein phylogenetic-based classification criterion. Genomic variables, such as mononucleotide and dinucleotide composition, codon usage bias, and minimal free energy derived from optimized predicted RNA secondary structure, were analyzed. From the values obtained, the main evolutionary pressures were discussed, indicating that natural selection plays a significant role in shaping mitovirus genomes. However, neutral evolution also makes a significant contribution. This study reveals a significant discovery of structural divergence in Kvaramitovirus. The energy minimization approach employed to study 2D folding in this study reveals a distinct spatial organization of their genomes, providing evidence for the hypothesis of a single evolutionary event of circularization in the most recent common ancestor of the lineage. This hypothesis was discussed in light of recent discoveries by other researchers that partially support the existence of mitoviruses with circular genomes. Finally, this study represents a significant advancement in the understanding of mitoviruses, as it quantitatively describes the nucleotide sequence at the family and genus taxonomic levels. Additionally, we provide hypotheses that can be experimentally validated to inspire new research and address the gaps in knowledge of this fascinating, basally divergent RNA virus lineage.

Sequence analysis of isolated strains of herpes zoster virus among patients with shingles

Background and Objectives

Herpes zoster, or shingles, is caused by the varicella-zoster virus (VZV), which initially presents as chickenpox in children. VZV is a global health concern, especially in winter and spring, affecting 10-20% of adults over 50 and posing a 30% risk for the general population. This study used PCR to detect VZV, confirming results with duplicated DNA samples and identifying 234 bp fragments by targeting the gpB gene.

Materials and Methods

This study examined 50 herpes zoster cases from October 2020 to April 2021, involving 30 males and 20 females aged 10 to 90, diagnosed by dermatologists. Data were collected via a questionnaire. PCR detected VZV by amplifying the gpB and MCP genes from skin lesion samples. Six positive 234-bp PCR products were sequenced at Macrogen Inc. in Seoul, South Korea.

Results

Six DNA samples with 234 bp amplicons were sequenced, showing 99-100% similarity to human alpha herpesvirus sequences in the gpB gene. NCBI BLAST matched these sequences to a reference (GenBank acc. MT370830.1), assigning accession numbers LC642111, LC642112, and LC642113. Eight nucleic acid substitutions caused amino acid changes in the gpB protein: isoleucine to threonine, serine to isoleucine, and threonine to Proline. These variants were deposited in NCBI GenBank as gpB3 samples.

Conclusion

The study found high sequence similarity to known VZV sequences, identifying six nucleic acid variations and eight SNPs. Notable amino acid changes in the gpB protein were deposited in NCBI GenBank as the gpB3 sample.

Evolution and diversity of the hepatitis B virus genome: Clinical implications

Hepatitis B virus (HBV) infection remains a significant global health burden. The genetic variation of HBV is complex. HBV can be divided into nine genotypes, which show significant differences in geographical distribution, clinical manifestations, transmission routes and treatment response. In recent years, substantial progress has been made through various research methods in understanding the development, pathogenesis, and antiviral treatment response of clinical disease associated with HBV genetic variants. This progress provides important theoretical support for a deeper understanding of the natural history of HBV infection, virus detection, drug treatment, vaccine development, mother-to-child transmission, and surveillance management. This review summarizes the mechanisms of HBV diversity, discusses methods used to detect viral diversity in current studies, and the impact of viral genome variation during infection on the development of clinical disease.

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.

All the single cells: if you like it then you should put some virus on it

Across a rich 70-year history, single-cell virology has revealed the impact of host and pathogen heterogeneity during virus infections. Recent technological innovations have enabled higher-resolution analyses of cellular and viral heterogeneity. Furthermore, single-cell analysis has revealed extreme phenotypes and provided additional insights into host-pathogen dynamics. Using a single-cell approach to explore fundamental virology questions, contemporary researchers have contributed to a revival of interest in single-cell virology with increased insights and enthusiasm.

The Characteristic of HBV Quasispecies Is Related to Occult HBV Infection of Infants Born to Highly Viremic Mothers

Although a combination of immunoprophylaxis and antiviral therapy can effectively prevent mother-to-child transmission (MTCT) of hepatitis B virus (HBV), a considerable number of infants born to highly viremic mothers still develop occult HBV infection (OBI). To uncover the virological factor and risk predictor for OBI in infants, we found that the diversity and complexity of maternal HBV quasispecies in the case group were lower than those in the control group. Mutations with significant differences between the two groups were most enriched in the NTCPbd and PreC regions. Genetic distance at the amino-acid level of the PreC region, especially the combination of three amino-acid mutations in the PreC region, could strongly predict the risk of OBI in infants. HBV quasispecies in OBI infants were highly complex, and the non-synonymous substitutions were mainly found in the RT and HBsAg regions. The sK47E (rtQ55R) and sP49L mutations in OBI infants might contribute to OBI through inhibiting the production of HBV DNA and HBsAg, respectively. This study found the potential virological factors and risk predictors for OBI in infants born to highly viremic mothers, which might be helpful for controlling OBI in infants.

In vitro and in vivo antiviral effects of CLEVir-X against porcine reproductive and respiratory syndrome virus

The aim of this study was to investigate the in vitro and in vivo antiviral effects of CLEVir-X, against porcine reproductive and respiratory syndrome virus (PRRSV). CLEVir-X is a nucleoside analogue and a dialdehyde form of xanthosine. CLEVir-X demonstrated antiviral action during the in vitro portion of this experiment with its inosine monophosphate dehydrogenase (IMPDH) inhibition against PRRSV. The anti-PRRSV effect of CLEVir-X was recovered through supplementation with guanosine. This suggests that PRRSV replication may be regulated through IMPDH and its guanosine biosynthetic pathway. CLEVir-X treatment in cultures resulted in mutation frequency increase of up to 7.8-fold within the viral genomes (e.g. ORF6) compared to their parallel, untreated cultures. The incorporation of CLEVir-X into the viral genome causes lethal mutagenesis and subsequent decrease in specific infectivity. During the in vivo antiviral experiment, 21-day-old pigs began oral administration of 5 mL of phosphate buffered saline containing CLEVir-X (with purity of 68 % and dosage of 40 mg/kg body weight). This treatment was provided twice daily at 9:00AM and 5:00PM for 14 days. Pigs were simultaneously intranasally inoculated with PRRSV at the beginning of CLEVir-X treatment (21 days of age). Several beneficial effects from the oral administration of CLEVir-X were observed including reduction of body temperature, alleviation of respiratory clinical signs, decreased PRRSV load in both blood and lung tissues, and mitigation of lung interstitial pneumonia lesions. The results of the present study demonstrated that CLEVir-X has mutagenic and nonmutagenic modes of antiviral action against PRRSV based on both in vitro and in vivo antiviral experiments.

Recombinant Viruses from the Picornaviridae Family Occurring in Racing Pigeons

Viruses from Picornaviridae family are known pathogens of poultry, although the information on their occurrence and pathogenicity in pigeons is scarce. In this research, efforts are made to broaden the knowledge on Megrivirus B and Pigeon picornavirus B prevalence, phylogenetic relationship with other avian picornaviruses and their possible connection with enteric disease in racing pigeons. As a result of Oxford Nanopore Sequencing, five Megrivirus and two pigeon picornavirus B-like genome sequences were recovered, among which three recombinant strains were detected. The recombinant fragments represented an average of 10.9% and 25.5% of the genome length of the Pigeon picornavirus B and Megrivirus B reference strains, respectively. The phylogenetic analysis revealed that pigeons are carriers of species-specific picornaviruses. TaqMan qPCR assays revealed 7.8% and 19.0% prevalence of Megrivirus B and 32.2% and 39.7% prevalence of Pigeon picornavirus B in the group of pigeons exhibiting signs of enteropathy and in the group of asymptomatic pigeons, respectively. In turn, digital droplet PCR showed a considerably higher number of genome copies of both viruses in sick than in asymptomatic pigeons. The results of quantitative analysis leave the role of picornaviruses in enteropathies of pigeons unclear.

RNA three-dimensional structure drives the sequence organization of potato spindle tuber viroid quasispecies

RNA viruses and viroids exist and evolve as quasispecies due to error-prone replication. Quasispecies consist of a few dominant master sequences alongside numerous variants that contribute to genetic diversity. Upon environmental changes, certain variants within quasispecies have the potential to become the dominant sequences, leading to the emergence of novel infectious strains. However, the emergence of new infectious variants remains unpredictable. Using mutant pools prepared by saturation mutagenesis of selected stem and loop regions, our study of potato spindle tuber viroid (PSTVd) demonstrates that mutants forming local three-dimensional (3D) structures similar to the wild type (WT) are more likely to accumulate in PSTVd quasispecies. The selection mechanisms underlying this biased accumulation are likely associated with cell-to-cell movement and long-distance trafficking. Moreover, certain trafficking-defective PSTVd mutants can be spread by functional sister genomes in the quasispecies. Our study reveals that the RNA 3D structure of stems and loops constrains the evolution of viroid quasispecies. Mutants with a structure similar to WT have a higher likelihood of being maintained within the quasispecies and can potentially give rise to novel infectious variants. These findings emphasize the potential of targeting RNA 3D structure as a more robust approach to defend against viroid infections.

Cell-cell communication and initial population composition shape the structure of potato spindle tuber viroid quasispecies

RNA viruses and viroids replicate with high mutation rates, forming quasispecies, population of variants centered around dominant sequences. The mechanisms governing quasispecies remain unclear. Plasmodesmata regulate viroid movement and were hypothesized to impact viroid quasispecies. Here, we sequenced the progeny of potato spindle tuber viroid intermediate (PSTVd-I) strain from mature guard cells lacking plasmodesmal connections and from in vitro-cultivated mesophyll cell protoplasts from systemic leaves of early-infected tomato (Solanum lycopersicum) plants. Remarkably, more variants accumulated in guard cells compared to whole leaves. Similarly, after extended cell culture, we observed more variants in cultivated mesophyll protoplasts. Coinfection and single-cell sequencing experiments demonstrated that the same plant cell can be infected multiple times by the same or different PSTVd sequences. To study the impact of initial population composition on PSTVd-I quasispecies, we conducted coinfections with PSTVd-I and variants. Two inoculum ratios (10:1 or 1:10) established quasispecies with or without PSTVd-I as the master sequence. In the absence of the master sequence, the percentage of novel variants initially increased. Moreover, a 1:1 PSTVd-I/variant RNA ratio resulted in PSTVd-I dominating (>50%), while the variants reached 20%. After PSTVd-I-only infection, the variants reached around 10%, while after variant-only infection, the variants were significantly more than 10%. These results emphasize the role of cell-to-cell communication and initial population composition in shaping PSTVd quasispecies.

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.

The dynamic variation position and predominant quasispecies of hepatitis B virus: Novel predictors of early hepatocarcinoma

To find the predictors of early HCC based on the dynamic changes of HBV quasispecies, this study utilizing the second-generation sequencing (NGS) and high-order multiplex droplet digital PCR (ddPCR) technology to examine the HBV quasispecies in serum of total 247 subjects recruited from high-incidence area of HCC. In the discovery stage, 15 non-synonymous Single Nucleotide Polymorphisms (SNPs) with higher variant proportion in HCC case group were founded (all P<0.05). Furthermore, the variant proportions in some of these SNPs were observed changing regularly within 5 years before the onset of HCC, and 5 of them located in HBX, 2 in HBS and 2 in HBC. The HBV predominant quasispecies and their consensus sequences were identified by genetic evolution analysis, in which the high HBS and HBC quasispecies heterogeneity were found associated with the forming of multifarious quasispecies clones, and the HBX gene had the highest proportion of predominant quasispecies (46.7 % in HBX vs 12.7 % and 13.8 % in HBS and HBC respectively) with the key variations (G1512A, A1630G, T1753C/G/A, A1762T and G1764A) determined. In the validation stage, we confirmed that the combined double mutations of G1512A+A1630G, A1762T+G1764A, and the combined triple mutations of T1753C/G/A + A1762T+G1764A, all expressed higher in early HCC cases when comparing with control group (all P<0.05). We also demonstrated the advantages of ddPCR using in multi-variations detection in large-sample for early HCC surveillance and screening. So we think that the dynamic of key HBV variation positions and their different combinations determined by quasispecies anlysis in this study can act as the novel predictors of early hepatocarcinoma and suitable to popularize and apply in HCC screening.

Quasispecies productivity

The quasispecies theory is a helpful concept in the explanation of RNA virus evolution and behaviour, with a relevant impact on methods used to fight viral diseases. It has undergone some adaptations to integrate new empirical data, especially the non-deterministic nature of mutagenesis, and the variety of behavioural motifs in cooperation, competition, communication, innovation, integration, and exaptation. Also, the consortial structure of quasispecies with complementary roles of memory genomes of minority populations better fits the empirical data than did the original concept of a master sequence and its mutant spectra. The high productivity of quasispecies variants generates unique sequences that never existed before and will never exist again. In the present essay, we underline that such sequences represent really new ontological entities, not just error copies of previous ones. Their primary unique property, the incredible variant production, is suggested here as quasispecies productivity, which replaces the error-replication narrative to better fit into a new relationship between mankind and living nature in the twenty-first century.

Evolutionary deletions within the SARS-CoV-2 genome as signature trends for virus fitness and adaptation

Coronaviruses are large RNA viruses that can infect and spread among humans and animals. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for coronavirus disease 2019, has evolved since its first detection in December 2019. Deletions are a common occurrence in SARS-CoV-2 evolution, particularly in specific genomic sites, and may be associated with the emergence of highly competent lineages. While deletions typically have a negative impact on viral fitness, some persist and become fixed in viral populations, indicating that they may confer advantageous benefits for the virus's adaptive evolution. This work presents a literature review and data analysis on structural losses in the SARS-CoV-2 genome and the potential relevance of specific signatures for enhanced viral fitness and spread.

Sequence analysis of the Spike, RNA-dependent RNA polymerase, and protease genes reveals a distinct evolutionary pattern of SARS-CoV-2 variants circulating in Yogyakarta and Central Java provinces, Indonesia

During the Covid-19 pandemic, the resurgence of SARS-CoV-2 was due to the development of novel variants of concern (VOC). Thus, genomic surveillance is essential to monitor continuing evolution of SARS-CoV-2 and to track the emergence of novel variants. In this study, we performed phylogenetic, mutation, and selection pressure analyses of the Spike, nsp12, nsp3, and nsp5 genes of SARS-CoV-2 isolates circulating in Yogyakarta and Central Java provinces, Indonesia from May 2021 to February 2022. Various bioinformatics tools were employed to investigate the evolutionary dynamics of distinct SARS-CoV-2 isolates. During the study period, 213 and 139 isolates of Omicron and Delta variants were identified, respectively. Particularly in the Spike gene, mutations were significantly more abundant in Omicron than in Delta variants. Consistently, in all of four genes studied, the substitution rates of Omicron were higher than that of Delta variants, especially in the Spike and nsp12 genes. In addition, selective pressure analysis revealed several sites that were positively selected in particular genes, implying that these sites were functionally essential for virus evolution. In conclusion, our study demonstrated a distinct evolutionary pattern of SARS-CoV-2 variants circulating in Yogyakarta and Central Java provinces, Indonesia.

The GA-Hecate Peptide inhibits the ZIKV Replicative Cycle in Different Steps and can Inhibit the Flavivirus NS2B-NS3 Protease after Cell Infection

BACKGROUND

Peptide drugs are advantageous because they are subject to rational design and exhibit highly diverse structures and broad biological activities. The NS2B-NS3 protein is a particularly promising flavivirus therapeutic target, with extensive research on the development of inhibitors as therapeutic candidates, and was used as a model in this work to determine the mechanism by which GA-Hecate inhibits ZIKV replication.

OBJECTIVE

The present study aimed to evaluate the potential of GA-Hecate, a new antiviral developed by our group, against the Brazilian Zika virus and to evaluate the mechanism of action of this compound on the flavivirus NS2B-NS3 protein.

METHODS

Solid-phase peptide Synthesis, High-Performance Liquid Chromatography, and Mass Spectrometry were used to obtain, purify, and characterize the synthesized compound. Real-time and enzymatic assays were used to determine the antiviral potential of GA-Hecate against ZIKV.

RESULTS

The RT-qPCR results showed that GA-Hecate decreased the number of ZIKV RNA copies in the virucidal, pre-treatment, and post-entry assays, with 5- to 6-fold fewer RNA copies at the higher nontoxic concentration in Vero cells (HNTC: 10 μM) than in the control cells. Enzymatic and kinetic assays indicated that GA-Hecate acts as a competitive ZIKV NS2B-NS3 protease inhibitor with an IC50 of 32 nM and has activity against the yellow fever virus protease.

CONCLUSION

The results highlight the antiviral potential of the GA-Hecate bioconjugate and open the door for the development of new antivirals.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

SARS-CoV-2 genomics and impact on clinical care for COVID-19

The emergence and worldwide spread of SARS-CoV-2 during the COVID-19 pandemic necessitated the adaptation and rapid deployment of viral WGS and analysis techniques that had been previously applied on a more limited basis to other viral pathogens, such as HIV and influenza viruses. The need for WGS was driven in part by the low mutation rate of SARS-CoV-2, which necessitated measuring variation along the entire genome sequence to effectively differentiate lineages and characterize viral evolution. Several WGS approaches designed to maximize throughput and accuracy were quickly adopted by surveillance labs around the world. These broad-based SARS-CoV-2 genomic sequencing efforts revealed ongoing evolution of the virus, highlighted by the successive emergence of new viral variants throughout the course of the pandemic. These genomic insights were instrumental in characterizing the effects of viral mutations on transmissibility, immune escape and viral tropism, which in turn helped guide public health policy, the use of monoclonal antibody therapeutics and vaccine development strategies. As the use of direct-acting antivirals for the treatment of COVID-19 became more widespread, the potential for emergence of antiviral resistance has driven ongoing efforts to delineate resistance mutations and to monitor global sequence databases for their emergence. Given the critical role of viral genomics in the international effort to combat the COVID-19 pandemic, coordinated efforts should be made to expand global genomic surveillance capacity and infrastructure towards the anticipation and prevention of future pandemics.

Characterisation of human astrovirus in a diarrhoea outbreak using nanopore and Sanger sequencing protocols

Human astroviruses (HAstV) are etiologic agents of acute gastroenteritis that most often afflict young children and elderly adults. Most studies of HAstV have focused on epidemiology. In this study, we collected 10 stool samples from a diarrhea outbreak from a diarrhea sentinel surveillance hospital in Beijing. Samples were evaluated immediately using parallel multiplex RT-qPCR and nanopore sequencing, and were then amplified by designed primers and Sanger sequencing to obtain whole genome sequences. Six isolates were categorized as HAstV-5 and subjected to whole genome analysis to characterize their genetic variation and evolution. Full genome analysis revealed low genetic variation (99.38-100% identity) among isolates. Phylogenetic analysis showed that all isolates were closely related to domestic strains Yu/1-CHN and 2013/Fuzhou/85. The recombination breakpoint of the six isolates was located at 2741 bp in the overlap region of ORF1a and ORF1b, similar to those of Yu/1-CHN and 2013/Fuzhou/85. Overall, our study highlights the combined use of RT-qPCR and sequencing as an important tool in rapid diagnosis and acquisition of whole genome sequences of HAstV.

From viral democratic genomes to viral wild bunch of quasispecies

The tremendous majority of RNA genomes from pathogenic viruses analyzed and deposited in databases are consensus or "democratic" genomes. They represent the genomes most frequently found in the clinical samples of patients but do not account for the huge genetic diversity of coexisting genomes, which is better described as quasispecies. A viral quasispecies is defined as the dynamic distribution of nonidentical but closely related mutants, variants, recombinant, or reassortant viral genomes. Viral quasispecies have collective behavior and dynamics and are the subject of internal interactions that comprise interference, complementation, or cooperation. In the setting of SARS-CoV-2 infection, intrahost SARS-CoV-2 genetic diversity was recently notably reported for immunocompromised, chronically infected patients, for patients treated with monoclonal antibodies targeting the viral spike protein, and for different body compartments of a single patient. A question that deserves attention is whether such diversity is generated postinfection from a clonal genome in response to selection pressure or is already present at the time of infection as a quasispecies. In the present review, we summarize the data supporting that hosts are infected by a "wild bunch" of viruses rather than by multiple virions sharing the same genome. Each virion in the "wild bunch" may have different virulence and tissue tropisms. As the number of viruses replicated during host infections is huge, a viral quasispecies at any time of infection is wide and is also influenced by host-specific selection pressure after infection, which accounts for the difficulty in deciphering and predicting the appearance of more fit variants and the evolution of epidemics of novel RNA viruses.

Open questions in the social lives of viruses

Social interactions among viruses occur whenever multiple viral genomes infect the same cells, hosts, or populations of hosts. Viral social interactions range from cooperation to conflict, occur throughout the viral world, and affect every stage of the viral lifecycle. The ubiquity of these social interactions means that they can determine the population dynamics, evolutionary trajectory, and clinical progression of viral infections. At the same time, social interactions in viruses raise new questions for evolutionary theory, providing opportunities to test and extend existing frameworks within social evolution. Many opportunities exist at this interface: Insights into the evolution of viral social interactions have immediate implications for our understanding of the fundamental biology and clinical manifestation of viral diseases. However, these opportunities are currently limited because evolutionary biologists only rarely study social evolution in viruses. Here, we bridge this gap by (1) summarizing the ways in which viruses can interact socially, including consequences for social evolution and evolvability; (2) outlining some open questions raised by viruses that could challenge concepts within social evolution theory; and (3) providing some illustrative examples, data sources, and conceptual questions, for studying the natural history of social viruses.

Reovirus genomic diversity confers plasticity for protease utility during adaptation to intracellular uncoating

IMPORTANCE

Reoviruses infect many mammals and are widely studied as a model system for enteric viruses. However, most of our reovirus knowledge comes from laboratory strains maintained on immortalized L929 cells. Herein, we asked whether naturally circulating reoviruses possess the same genetic and phenotypic characteristics as laboratory strains. Naturally circulating reoviruses obtained from sewage were extremely diverse genetically. Moreover, sewage reoviruses exhibited poor fitness on L929 cells and relied heavily on gut proteases for viral uncoating and productive infection compared to laboratory strains. We then examined how naturally circulating reoviruses might adapt to cell culture conditions. Within three passages, virus isolates from the parental sewage population were selected, displaying improved fitness and intracellular uncoating in L929 cells. Remarkably, selected progeny clones were present at 0.01% of the parental population. Altogether, using reovirus as a model, our study demonstrates how the high genetic diversity of naturally circulating viruses results in rapid adaptation to new environments.

Diversity and Complexity of Internally Deleted Viral Genomes in Influenza A Virus Subpopulations with Enhanced Interferon-Inducing Phenotypes

Influenza A virus (IAV) populations harbor large subpopulations of defective-interfering particles characterized by internally deleted viral genomes. These internally deleted genomes have demonstrated the ability to suppress infectivity and boost innate immunity, rendering them promising for therapeutic and immunogenic applications. In this study, we aimed to investigate the diversity and complexity of the internally deleted IAV genomes within a panel of plaque-purified avian influenza viruses selected for their enhanced interferon-inducing phenotypes. Our findings unveiled that the abundance and diversity of internally deleted viral genomes were contingent upon the viral subculture and plaque purification processes. We observed a heightened occurrence of internally deleted genomes with distinct junctions in viral clones exhibiting enhanced interferon-inducing phenotypes, accompanied by additional truncation in the nonstructural 1 protein linker region (NS1Δ76-86). Computational analyses suggest the internally deleted IAV genomes can encode a broad range of carboxy-terminally truncated and intrinsically disordered proteins with variable lengths and amino acid composition. Further research is imperative to unravel the underlying mechanisms driving the increased diversity of internal deletions within the genomes of viral clones exhibiting enhanced interferon-inducing capacities and to explore their potential for modulating cellular processes and immunity.

Whole Genome Sequencing of Infectious Bursal Disease Viruses Isolated from a Californian Outbreak Unravels the Underlying Virulence Markers and Highlights Positive Selection Incidence

Outbreaks of the immunosuppressive infectious bursal disease (IBD) are frequently reported worldwide, despite the vaccination regimes. A 2009 Californian IBD outbreak caused by rA and rB isolates was described as very virulent (vv) IBD virus (IBDV); however, molecular factors beyond this virulence were not fully uncovered. Therefore, segments of both isolates were amplified, successfully cloned, whole genome sequenced by Next Generation Sequencing, genotyped, and the leading virulence factors were entirely investigated in terms of phylogenetic and amino acid analysis and protein modeling for positive selection orientation and interaction analysis. rA and rB isolates displayed the highest amino acid identity (97.84-100%) with Genotype 3 strains. Interestingly, rA and rB contained all virulence hallmarks of hypervariable (HVR), including 222A, 242I, 249Q, 256I, 284A, 286T, 294I, 299S, and 318G, as well as the serine-rich heptapeptide sequence. Moreover, we pinpointed the A3B2 genotype of rA and rB, predominant in non-reassortants, and we highlighted the absence of recombination events. Furthermore, gene-wise phylogenetic analysis showed the entire genes of rA and rB clustered with the vvIBDVs and emphasized their share in IBDV virulence. VP5 showed a virulence marker, MLSL (amino acid sequence). VP2 encountered three significant novel mutations apart from the HVR, including G163E in rA and Y173C and V178A in rB, all residing within interacting motifs. VP4 contained 168Y, 173N, 203S, and 239D characteristic for the vv phenotype. A235V mutation was detected at the dsRNA binding domain of VP3. In VP1, the TDN triplet and the mutation (V4I) were detected, characteristic of hypervirulence occurring at the N-terminus responsible for protein priming. Although selection analysis revealed seven sites, codon 222 was the only statistically significant selection site. The VP2 modeling of rA and rB highlighted great structure fitness, with 96.14% Ramachandran favored positioning including the 222A, i.e., not influencing the structure stability. The 222A was found to be non-interface surface residue, associated with no interaction with the attachment-mediated ligand motif. Our findings provide pivotal insights into the evolution and underlying virulence factors and will assist in the development of control strategies via sequence-based continuous monitoring for the early detection of novel vv strains.

Establishment and Application of CRISPR-Cas12a-Based Recombinase Polymerase Amplification and a Lateral Flow Dipstick and Fluorescence for the Detection and Distinction of Deformed Wing Virus Types A and B

Deformed wing virus (DWV) is one of the important pathogens of the honey bee (Apis mellifera), which consists of three master variants: types A, B, and C. Among them, DWV types A (DWV-A) and B (DWV-B) are the most prevalent variants in honey bee colonies and have been linked to colony decline. DWV-A and DWV-B have different virulence, but it is difficult to distinguish them via traditional methods. In this study, we established a visual detection assay for DWV-A and DWV-B using recombinase polymerase amplification (RPA) and a lateral flow dipstick (LFD) coupled with the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein (Cas) 12a fluorescence system (RPA-CRISPR-Cas12a-LFD). The limit of detection of this system was ~6.5 × 100 and 6.2 × 101 copies/μL for DWV-A and DWV-B, respectively. The assays were specific and non-cross-reactive against other bee viruses, and the results could be visualized within 1 h. The assays were validated by extracting cDNA from 36 clinical samples of bees that were suspected to be infected with DWV. The findings were consistent with those of traditional reverse transcription-quantitative polymerase chain reaction, and the RPA-CRISPR-Cas12a assay showed the specific, sensitive, simple, and appropriate detection of DWV-A and DWV-B. This method can facilitate the visual and qualitative detection of DWV-A and DWV-B as well as the monitoring of different subtypes, thereby providing potentially better control and preventing current and future DWV outbreaks.

Genetic diversification patterns in swine influenza A virus (H1N2) in vaccinated and nonvaccinated animals

Influenza A viruses (IAVs) are characterized by having a segmented genome, low proofreading polymerases, and a wide host range. Consequently, IAVs are constantly evolving in nature causing a threat to animal and human health. In 2009 a new human pandemic IAV strain arose in Mexico because of a reassortment between two strains previously circulating in pigs; Eurasian "avian-like" (EA) swine H1N1 and "human-like" H1N2, highlighting the importance of swine as adaptation host of avian to human IAVs. Nowadays, although of limited use, a trivalent vaccine, which include in its formulation H1N1, H3N2, and, H1N2 swine IAV (SIAV) subtypes, is one of the most applied strategies to reduce SIAV circulation in farms. Protection provided by vaccines is not complete, allowing virus circulation, potentially favoring viral evolution. The evolutionary dynamics of SIAV quasispecies were studied in samples collected at different times from 8 vaccinated and 8 nonvaccinated pigs, challenged with H1N2 SIAV. In total, 32 SIAV genomes were sequenced by next-generation sequencing, and subsequent variant-calling genomic analysis was carried out. Herein, a total of 364 de novo single nucleotide variants (SNV) were found along all genetic segments in both experimental groups. The nonsynonymous substitutions proportion found was greater in vaccinated animals suggesting that H1N2 SIAV was under positive selection in this scenario. The impact of each substitution with an allele frequency greater than 5% was hypothesized according to previous literature, particularly in the surface glycoproteins hemagglutinin and neuraminidase. The H1N2 SIAV quasispecies evolution capacity was evidenced, observing different evolutionary trends in vaccinated and nonvaccinated animals.

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.

Influenza A Virus Multicycle Replication Yields Comparable Viral Population Emergence in Human Respiratory and Ocular Cell Types

While primarily considered a respiratory pathogen, influenza A virus (IAV) is nonetheless capable of spreading to, and replicating in, numerous extrapulmonary tissues in humans. However, within-host assessments of genetic diversity during multicycle replication have been largely limited to respiratory tract tissues and specimens. As selective pressures can vary greatly between anatomical sites, there is a need to examine how measures of viral diversity may vary between influenza viruses exhibiting different tropisms in humans, as well as following influenza virus infection of cells derived from different organ systems. Here, we employed human primary tissue constructs emulative of the human airway or corneal surface, and we infected both with a panel of human- and avian-origin IAV, inclusive of H1 and H3 subtype human viruses and highly pathogenic H5 and H7 subtype viruses, which are associated with both respiratory disease and conjunctivitis following human infection. While both cell types supported productive replication of all viruses, airway-derived tissue constructs elicited greater induction of genes associated with antiviral responses than did corneal-derived constructs. We used next-generation sequencing to examine viral mutations and population diversity, utilizing several metrics. With few exceptions, generally comparable measures of viral diversity and mutational frequency were detected following homologous virus infection of both respiratory-origin and ocular-origin tissue constructs. Expansion of within-host assessments of genetic diversity to include IAV with atypical clinical presentations in humans or in extrapulmonary cell types can provide greater insight into understanding those features most prone to modulation in the context of viral tropism. IMPORTANCE Influenza A virus (IAV) can infect tissues both within and beyond the respiratory tract, leading to extrapulmonary complications, such as conjunctivitis or gastrointestinal disease. Selective pressures governing virus replication and induction of host responses can vary based on the anatomical site of infection, yet studies examining within-host assessments of genetic diversity are typically only conducted in cells derived from the respiratory tract. We examined the contribution of influenza virus tropism on these properties two different ways: by using IAV associated with different tropisms in humans, and by infecting human cell types from two different organ systems susceptible to IAV infection. Despite the diversity of cell types and viruses employed, we observed generally similar measures of viral diversity postinfection across all conditions tested; these findings nonetheless contribute to a greater understanding of the role tissue type contributes to the dynamics of virus evolution within a human host.

ACE2 receptor polymorphism in humans and animals increases the risk of the emergence of SARS-CoV-2 variants during repeated intra- and inter-species host-switching of the virus

Like other coronaviruses, SARS-CoV-2 has ability to spread through human-to-human transmission and to circulate from humans to animals and from animals to humans. A high frequency of SARS-CoV-2 mutations has been observed in the viruses isolated from both humans and animals, suggesting a genetic fitness under positive selection in both ecological niches. The most documented positive selection force driving SARS-CoV-2 mutations is the host-specific immune response. However, after electrostatic interactions with lipid rafts, the first contact between the virus and host proteins is the viral spike-cellular receptor binding. Therefore, it is likely that the first level of selection pressure impacting viral fitness relates to the virus's affinity for its receptor, the angiotensin I converting enzyme 2 (ACE2). Although sufficiently conserved in a huge number of species to support binding of the viral spike with enough affinity to initiate fusion, ACE2 is highly polymorphic both among species and within a species. Here, we provide evidence suggesting that when the viral spike-ACE2 receptor interaction is not optimal, due to host-switching, mutations can be selected to improve the affinity of the spike for the ACE2 expressed by the new host. Notably, SARS-CoV-2 is mutation-prone in the spike receptor binding domain (RBD), allowing a better fit for ACE2 orthologs in animals. It is possibly that this may also be true for rare human alleles of ACE2 when the virus is spreading to billions of people. In this study, we present evidence that human subjects expressing the rare E₃₂₉G allele of ACE2 with higher allele frequencies in European populations exhibit a improved affinity for the SARS-CoV-2 spike N₅₀₁Y variant of the virus. This may suggest that this viral N₅₀₁Y variant emerged in the human population after SARS-CoV-2 had infected a human carrying the rare E₃₂₉G allele of ACE2. In addition, this viral evolution could impact viral replication as well as the ability of the adaptive humoral response to control infection with RBD-specific neutralizing antibodies. In a shifting landscape, this ACE2-driven genetic drift of SARS-CoV-2 which we have named the 'boomerang effect', could complicate the challenge of preventing COVID with a SARS-CoV-2 spike-derived vaccine.