Rates of evolutionary change in viruses: patterns and determinants

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.

Analysis of the genetic diversity in RNA-directed RNA polymerase sequences: implications for an automated RNA virus classification system

RNA viruses are characterized by a broad host range and high levels of genetic diversity. Despite a recent expansion in the known virosphere following metagenomic sequencing, our knowledge of the species rank genetic diversity of RNA viruses, and how often they are misassigned and misclassified, is limited. We performed a clustering analysis of 7801 RNA-directed RNA polymerase (RdRp) sequences representing 1897 established RNA virus species. From this, we identified substantial genetic divergence within some virus species and inconsistency in RNA virus assignment between the GenBank database and The International Committee on Taxonomy of Viruses (ICTV). In particular, 27.57% virus species comprised multiple virus operational taxonomic units (vOTUs), including Alphainfluenzavirus influenzae, Mammarenavirus lassaense, Apple stem pitting virus, and Rotavirus A, with each having over 100 vOTUs. In addition, the distribution of average amino acid identity between vOTUs within single assigned species showed a relatively low threshold: <90% and sometimes <50%. However, when only exemplar sequences from virus species were analyzed, 1889 of the ICTV-designated RNA virus species (99.58%) were clustered into a single vOTU. Clustering of the RdRp sequences from different virus species also revealed that 17 vOTUs contained two distinct virus species. These potential misassignments were confirmed by phylogenetic analysis. A further analysis of average nucleotide identity (ANI) values ranging from 70% to 97.5% revealed that at an ANI of 82.5%, 1559 (82.18%) of the 1897 virus species could be correctly clustered into one single vOTU. However, at ANI values >82.5%, an increasing number of species were clustered into two or more vOTUs. In sum, we have identified some inconsistency and misassignment of the RNA virus species based on the analysis of RdRp sequences alone, which has important implications for the development of an automated RNA virus classification system.

Computational Analysis of the Accumulation of Mutations in Therapeutically Important RNA Viral Proteins During Pandemics with Special Emphasis on SARS-CoV-2

Single stranded RNA viruses are primary causative agents for pandemics, causing extensive morbidity and mortality worldwide. A pivotal question in pandemic preparedness and therapeutic intervention is what are the specific mutations which are more likely to emerge during such global health crises? This study aims to identify markers for mutations with the highest probability of emergence in these pandemics, focusing on the SARS-CoV-2 spike protein, an essential and therapeutically significant viral protein, starting from sequence information from the onset of the pandemic until July 2022. Quite consistently, we observed that emerged mutations tended to demonstrate a high genetic score, which reflects high similarity of the type of codon required for translation between an amino acid and to the mutated one. Further, this pattern is also observed in therapeutically significant proteins of other ssRNA pandemic viruses, including influenza (HA, NA), spike proteins of Ebola, envelope of Dengue and Chikungunya. We propose that the genetic score serves as an initial indicator, preceding the actual impact of the mutation on viral fitness. Finally, we developed a comprehensive computational pipeline to further explore and predict the subsequent effects of mutations on viral fitness. We believe that our pipeline can narrow down and predict future mutations in therapeutically important viral proteins during a pandemic.

Host-Virus Cophylogenetic Trajectories: Investigating Molecular Relationships between Coronaviruses and Bat Hosts

Bats, with their virus tolerance, social behaviors, and mobility, are reservoirs for emerging viruses, including coronaviruses (CoVs) known for genetic flexibility. Studying the cophylogenetic link between bats and CoVs provides vital insights into transmission dynamics and host adaptation. Prior research has yielded valuable insights into phenomena such as host switching, cospeciation, and other dynamics concerning the interaction between CoVs and bats. Nonetheless, a distinct gap exists in the current literature concerning a comparative cophylogenetic analysis focused on elucidating the contributions of sequence fragments to the co-evolution between hosts and viruses. In this study, we analyzed the cophylogenetic patterns of 69 host-virus connections. Among the 69 host-virus links examined, 47 showed significant cophylogeny based on ParaFit and PACo analyses, affirming strong associations. Focusing on two proteins, ORF1ab and spike, we conducted a comparative analysis of host and CoV phylogenies. For ORF1ab, the specific window ranged in multiple sequence alignment (positions 520-680, 770-870, 2930-3070, and 4910-5080) exhibited the lowest Robinson-Foulds (RF) distance (i.e., 84.62%), emphasizing its higher contribution in the cophylogenetic association. Similarly, within the spike region, distinct window ranges (positions 0-140, 60-180, 100-410, 360-550, and 630-730) displayed the lowest RF distance at 88.46%. Our analysis identified six recombination regions within ORF1ab (positions 360-1390, 550-1610, 680-1680, 700-1710, 2060-3090, and 2130-3250), and four within the spike protein (positions 10-510, 50-560, 170-710, and 230-730). The convergence of minimal RF distance regions with combination regions robustly affirms the pivotal role of recombination in viral adaptation to host selection pressures. Furthermore, horizontal gene transfer reveals prominent instances of partial gene transfer events, occurring not only among variants within the same host species but also crossing host species boundaries. This suggests a more intricate pattern of genetic exchange. By employing a multifaceted approach, our comprehensive strategy offers a nuanced understanding of the intricate interactions that govern the co-evolutionary dynamics between bat hosts and CoVs. This deeper insight enhances our comprehension of viral evolution and adaptation mechanisms, shedding light on the broader dynamics that propel viral diversity.

Interrogating Genomes and Geography to Unravel Multiyear Vesicular Stomatitis Epizootics

We conducted an integrative analysis to elucidate the spatial epidemiological patterns of the Vesicular Stomatitis New Jersey virus (VSNJV) during the 2014-15 epizootic cycle in the United States (US). Using georeferenced VSNJV genomics data, confirmed vesicular stomatitis (VS) disease cases from surveillance, and a suite of environmental factors, our study assessed environmental and phylogenetic similarity to compare VS cases reported in 2014 and 2015. Despite uncertainties from incomplete virus sampling and cross-scale spatial processes, patterns suggested multiple independent re-invasion events concurrent with potential viral overwintering between sequential seasons. Our findings pointed to a geographically defined southern virus pool at the US-Mexico interface as the source of VSNJV invasions and overwintering sites. Phylodynamic analysis demonstrated an increase in virus diversity before a rise in case numbers and a pronounced reduction in virus diversity during the winter season, indicative of a genetic bottleneck and a significant narrowing of virus variation between the summer outbreak seasons. Environment-vector interactions underscored the central role of meta-population dynamics in driving disease spread. These insights emphasize the necessity for location- and time-specific management practices, including rapid response, movement restrictions, vector control, and other targeted interventions.

COVID-19 vaccination: 2023 Taiwan Association of Gerontology and Geriatrics (TAGG) consensus statements

The COVID-19 pandemic remains challenging due to the rapid evolution of the severe acute respiratory syndrome coronavirus 2. This article discusses recent findings on high-risk groups for COVID-19 mortality and morbidity, along with consensus statements from the 2023 Taiwan Association of Gerontology and Geriatrics (TAGG) meeting. It examines evidence on viral mutation mechanisms, emerging variants, and their implications for vaccination strategies. The article underscores advanced age, immunocompromised status, chronic medical conditions, occupational exposure, and socioeconomic disparities as significant risk factors for severe COVID-19 outcomes. TAGG's consensus emphasizes robust vaccination promotion, prioritizing elderly, and immunocompromised groups, individualized multi-dose regimens for immunocompromised patients, and simplified clinical guidelines. Discussions on global and regional recommendations for regular, variant-adapted boosters highlight the non-seasonal nature of COVID-19. Key agreements include escalating domestic preparedness, implementing vigorous risk-based vaccination, and adapting global guidelines to local contexts. Given ongoing viral evolution, proactive adjustment of vaccination policies is essential. Scientific consensus, tailored recommendations, and rapid knowledge dissemination are vital for optimizing COVID-19 protection among vulnerable groups in Taiwan. This article seeks to inform clinical practice and public health policy by summarizing expert-driven vaccination perspectives.

Allosteric Regulation and Inhibition of Coronavirus 3CLpro Revealed by HDX-MS

Coronavirus (CoV) infections have caused contagious and fatal respiratory diseases in humans worldwide. CoV 3-chymotrypsin-like proteases (3CLpro or Mpro) play an important role in viral maturation, and maintenance of their dimeric conformation is crucial for viral activity. Therefore, allosterically regulated dimerization of 3CLpro can be employed as a drug development target. Here, we investigated the allosteric regulatory mechanism of 3CLpro dimerization by using hydrogen/deuterium exchange coupled with mass spectrometry (HDX-MS) technology. We found that the FLAG tag directly coupled to the N-finger of 3CLpro significantly increased HDX kinetics at the dimer interface, and 3CLpro transformed from a dimer to a monomer. The 3CLpro mutants of SARS-CoV-2, which are monomeric, also exhibited increased deuterium exchange. Binding of the allosteric inhibitor Gastrodenol to most betacoronavirus 3CLpros led to increased allosteric deuterium exchange, resulting in the monomeric conformation of the CoV 3CLpro upon binding. Molecular dynamics (MD) simulation analysis further indicated the molecular mechanism of action of Gastrodenol on CoV 3CLpro: binding of Gastrodenol to SARS-CoV-2 3CLpro destroyed the hydrogen bond in the dimer interface. These results suggest that Gastrodenol may be a potential broad-spectrum anti-betacoronavirus drug.

Molecular epidemiology of recurrent zoonotic transmission of mpox virus in West Africa

Nigeria and Cameroon reported their first mpox cases in over three decades in 2017 and 2018 respectively. The outbreak in Nigeria is recognised as an ongoing human epidemic. However, owing to sparse surveillance and genomic data, it is not known whether the increase in cases in Cameroon is driven by zoonotic or sustained human transmission. Notably, the frequency of zoonotic transmission remains unknown in both Cameroon and Nigeria. To address these uncertainties, we investigated the zoonotic transmission dynamics of the mpox virus (MPXV) in Cameroon and Nigeria, with a particular focus on the border regions. We show that in these regions mpox cases are still driven by zoonotic transmission of a newly identified Clade IIb.1. We identify two distinct zoonotic lineages that circulate across the Nigeria-Cameroon border, with evidence of recent and historic cross border dissemination. Our findings support that the complex cross-border forest ecosystems likely hosts shared animal populations that drive cross-border viral spread, which is likely where extant Clade IIb originated. We identify that the closest zoonotic outgroup to the human epidemic circulated in southern Nigeria in October 2013. We also show that the zoonotic precursor lineage circulated in an animal population in southern Nigeria for more than 45 years. This supports findings that southern Nigeria was the origin of the human epidemic. Our study highlights the ongoing MPXV zoonotic transmission in Cameroon and Nigeria, underscoring the continuous risk of MPXV (re)emergence.

Mutations in the SARS-CoV-2 spike receptor binding domain and their delicate balance between ACE2 affinity and antibody evasion

Intensive selection pressure constrains the evolutionary trajectory of SARS-CoV-2 genomes and results in various novel variants with distinct mutation profiles. Point mutations, particularly those within the receptor binding domain (RBD) of SARS-CoV-2 spike (S) protein, lead to the functional alteration in both receptor engagement and monoclonal antibody (mAb) recognition. Here, we review the data of the RBD point mutations possessed by major SARS-CoV-2 variants and discuss their individual effects on ACE2 affinity and immune evasion. Many single amino acid substitutions within RBD epitopes crucial for the antibody evasion capacity may conversely weaken ACE2 binding affinity. However, this weakened effect could be largely compensated by specific epistatic mutations, such as N501Y, thus maintaining the overall ACE2 affinity for the spike protein of all major variants. The predominant direction of SARS-CoV-2 evolution lies neither in promoting ACE2 affinity nor evading mAb neutralization but in maintaining a delicate balance between these two dimensions. Together, this review interprets how RBD mutations efficiently resist antibody neutralization and meanwhile how the affinity between ACE2 and spike protein is maintained, emphasizing the significance of comprehensive assessment of spike mutations.

Molecular evolutionary characteristics of severe acute respiratory syndrome coronavirus 2 and the relatedness of epidemiological and socio-environmental factors

After the first outbreak, SARS-CoV-2 infection continues to occur due to the emergence of new variants. There is limited information available on the comparative evaluation of evolutionary characteristics of SARS-CoV-2 among different countries over time, and its relatedness to epidemiological and socio-environmental factors within those countries. We assessed comparative Bayesian evolutionary characteristics for SARS-CoV-2 in eight countries from 2020 to 2022 using BEAST version 2.6.7. Additionally, the relatedness between virus evolution factors and both epidemiological and socio-environmental factors was analyzed using Pearson's correlation coefficient. The estimated substitution rates in the gene encoding S protein of SARS-CoV-2 exhibited a continuous increase from 2020 to 2022 and were divided into two distinct groups in 2022 (p value < 0.05). Effective population size (Ne) generally showed decreased patterns by time. Notably, the change rates of the substitution rates were negatively correlated with the cumulative vaccination rates in 2021. A strict and rapid vaccination policy in the United Arab Emirates dramatically reduced the evolution of the virus, compared to other countries. Also, the average yearly temperature in countries were negatively correlated with the substitution rates. The changes of six epitopes in SARS-CoV-2 were related to various socio-environmental factors. We figured out comparative virus evolutionary traits and the association of epidemiological and socio-environmental factors especially cumulative vaccination rates and average temperature.

The human microbiota is a beneficial reservoir for SARS-CoV-2 mutations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations are rapidly emerging. In particular, beneficial mutations in the spike (S) protein, which can either make a person more infectious or enable immunological escape, are providing a significant obstacle to the prevention and treatment of pandemics. However, how the virus acquires a high number of beneficial mutations in a short time remains a mystery. We demonstrate here that variations of concern may be mutated due in part to the influence of the human microbiome. We searched the National Center for Biotechnology Information database for homologous fragments (HFs) after finding a mutation and the six neighboring amino acids in a viral mutation fragment. Among the approximate 8,000 HFs obtained, 61 mutations in S and other outer membrane proteins were found in bacteria, accounting for 62% of all mutation sources, which is 12-fold higher than the natural variable proportion. A significant proportion of these bacterial species-roughly 70%-come from the human microbiota, are mainly found in the lung or gut, and share a composition pattern with COVID-19 patients. Importantly, SARS-CoV-2 RNA-dependent RNA polymerase replicates corresponding bacterial mRNAs harboring mutations, producing chimeric RNAs. SARS-CoV-2 may collectively pick up mutations from the human microbiota that change the original virus's binding sites or antigenic determinants. Our study clarifies the evolving mutational mechanisms of SARS-CoV-2.

IMPORTANCE

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations are rapidly emerging, in particular advantageous mutations in the spike (S) protein, which either increase transmissibility or lead to immune escape and are posing a major challenge to pandemic prevention and treatment. However, how the virus acquires a high number of advantageous mutations in a short time remains a mystery. Here, we provide evidence that the human microbiota is a reservoir of advantageous mutations and aids mutational evolution and host adaptation of SARS-CoV-2. Our findings demonstrate a conceptual breakthrough on the mutational evolution mechanisms of SARS-CoV-2 for human adaptation. SARS-CoV-2 may grab advantageous mutations from the widely existing microorganisms in the host, which is undoubtedly an "efficient" manner. Our study might open a new perspective to understand the evolution of virus mutation, which has enormous implications for comprehending the trajectory of the COVID-19 pandemic.

Accumulation Dynamics of Defective Genomes during Experimental Evolution of Two Betacoronaviruses

Virus-encoded replicases often generate aberrant RNA genomes, known as defective viral genomes (DVGs). When co-infected with a helper virus providing necessary proteins, DVGs can multiply and spread. While DVGs depend on the helper virus for propagation, they can in some cases disrupt infectious virus replication, impact immune responses, and affect viral persistence or evolution. Understanding the dynamics of DVGs alongside standard viral genomes during infection remains unclear. To address this, we conducted a long-term experimental evolution of two betacoronaviruses, the human coronavirus OC43 (HCoV-OC43) and the murine hepatitis virus (MHV), in cell culture at both high and low multiplicities of infection (MOI). We then performed RNA-seq at regular time intervals, reconstructed DVGs, and analyzed their accumulation dynamics. Our findings indicate that DVGs evolved to exhibit greater diversity and abundance, with deletions and insertions being the most common types. Notably, some high MOI deletions showed very limited temporary existence, while others became prevalent over time. We observed differences in DVG abundance between high and low MOI conditions in HCoV-OC43 samples. The size distribution of HCoV-OC43 genomes with deletions differed between high and low MOI passages. In low MOI lineages, short and long DVGs were the most common, with an additional cluster in high MOI lineages which became more prevalent along evolutionary time. MHV also showed variations in DVG size distribution at different MOI conditions, though they were less pronounced compared to HCoV-OC43, suggesting a more random distribution of DVG sizes. We identified hotspot regions for deletions that evolved at a high MOI, primarily within cistrons encoding structural and accessory proteins. In conclusion, our study illustrates the widespread formation of DVGs during betacoronavirus evolution, influenced by MOI and cell- and virus-specific factors.

Bacterial and viral co-infections in aquaculture under climate warming: co-evolutionary implications, diagnosis, and treatment

Climate change and the associated environmental temperature fluctuations are contributing to increases in the frequency and severity of disease outbreaks in both wild and farmed aquatic species. This has a significant impact on biodiversity and also puts global food production systems, such as aquaculture, at risk. Most infections are the result of complex interactions between multiple pathogens, and understanding these interactions and their co-evolutionary mechanisms is crucial for developing effective diagnosis and control strategies. In this review, we discuss current knowledge on bacteria-bacteria, virus-virus, and bacterial and viral co-infections in aquaculture as well as their co-evolution in the context of global warming. We also propose a framework and different novel methods (e.g. advanced molecular tools such as digital PCR and next-generation sequencing) to (1) precisely identify overlooked co-infections, (2) gain an understanding of the co-infection dynamics and mechanisms by knowing species interactions, and (3) facilitate the development multi-pathogen preventive measures such as polyvalent vaccines. As aquaculture disease outbreaks are forecasted to increase both due to the intensification of practices to meet the protein demand of the increasing global population and as a result of global warming, understanding and treating co-infections in aquatic species has important implications for global food security and the economy.

Geographical Expansion of Avian Metapneumovirus Subtype B: First Detection and Molecular Characterization of Avian Metapneumovirus Subtype B in US Poultry

Avian metapneumovirus (aMPV), classified within the Pneumoviridae family, wreaks havoc on poultry health. It typically causes upper respiratory tract and reproductive tract infections, mainly in turkeys, chickens, and ducks. Four subtypes of AMPV (A, B, C, D) and two unclassified subtypes have been identified, of which subtypes A and B are widely distributed across the world. In January 2024, an outbreak of severe respiratory disease occurred on turkey and chicken farms across different states in the US. Metagenomics sequencing of selected tissue and swab samples confirmed the presence of aMPV subtype B. Subsequently, all samples were screened using an aMPV subtype A and B multiplex real-time RT-PCR kit. Of the 221 farms, 124 (56%) were found to be positive for aMPV-B. All samples were negative for subtype A. Six whole genomes were assembled, five from turkeys and one from chickens; all six assembled genomes showed 99.29 to 99.98% nucleotide identity, indicating a clonal expansion event for aMPV-B within the country. In addition, all six sequences showed 97.74 to 98.58% nucleotide identity with previously reported subtype B sequences, e.g., VCO3/60616, Hungary/657/4, and BR/1890/E1/19. In comparison to these two reference strains, the study sequences showed unique 49-62 amino acid changes across the genome, with maximum changes in glycoprotein (G). One unique AA change from T (Threonine) to I (Isoleucine) at position 153 in G protein was reported only in the chicken aMPV sequence, which differentiated it from turkey sequences. The twelve unique AA changes along with change in polarity of the G protein may indicate that these unique changes played a role in the adaptation of this virus in the US poultry. This is the first documented report of aMPV subtype B in US poultry, highlighting the need for further investigations into its genotypic characterization, pathogenesis, and evolutionary dynamics.

Mottle: Accurate pairwise substitution distance at high divergence through the exploitation of short-read mappers and gradient descent

Current tools for estimating the substitution distance between two related sequences struggle to remain accurate at a high divergence. Difficulties at distant homologies, such as false seeding and over-alignment, create a high barrier for the development of a stable estimator. This is especially true for viral genomes, which carry a high rate of mutation, small size, and sparse taxonomy. Developing an accurate substitution distance measure would help to elucidate the relationship between highly divergent sequences, interrogate their evolutionary history, and better facilitate the discovery of new viral genomes. To tackle these problems, we propose an approach that uses short-read mappers to create whole-genome maps, and gradient descent to isolate the homologous fraction and calculate the final distance value. We implement this approach as Mottle. With the use of simulated and biological sequences, Mottle was able to remain stable to 0.66-0.96 substitutions per base pair and identify viral outgroup genomes with 95% accuracy at the family-order level. Our results indicate that Mottle performs as well as existing programs in identifying taxonomic relationships, with more accurate numerical estimation of genomic distance over greater divergences. By contrast, one limitation is a reduced numerical accuracy at low divergences, and on genomes where insertions and deletions are uncommon, when compared to alternative approaches. We propose that Mottle may therefore be of particular interest in the study of viruses, viral relationships, and notably for viral discovery platforms, helping in benchmarking of homology search tools and defining the limits of taxonomic classification methods. The code for Mottle is available at https://github.com/tphoward/Mottle_Repo.

Increasing evolution, prevalence, and outbreaks for rift valley fever virus in the process of breaking geographical barriers

BACKGROUND

Rift valley fever (RVF) is listed as one of prioritized diseases by WHO. This study aims to describe RVF virus' landscape distribution globally, and to insight dynamics change of its evolution, prevalence, and outbreaks in the process of breaking geographical barriers.

METHODS

A systematic literature review and meta-analyses was conducted to estimate RVF prevalence by hosts using a random-effect model. Molecular clock-based phylogenetic analyses were performed to estimate RVF virus nucleotide substitution rates using nucleotide sequences in NCBI database. RVF virus prevalence, nucleotide substitution rates, and outbreaks were compared before and after breaking geographical barriers twice, respectively.

RESULTS

RVF virus was reported from 26 kinds of hosts covering 48 countries from 1930 to 2022. Since RVF broke geographical barriers, (1) nucleotide substitution rates significantly increased after firstly spreading out of Africa in 2000, (2) prevalence in humans significantly increased from 1.92 % (95 % CI: 0.86-3.25 %) to 3.03 % (95 % CI: 2.09-4.12 %) after it broke Sahara Desert geographical barriers in 1977, and to 5.24 % (95 % CI: 3.81-6.82 %) after 2000, (3) RVF outbreaks in humans and the number of wildlife hosts presented increasing trends. RVF virus spillover may exist between bats and humans, and accelerate viral substitution rates in humans. During outbreaks, the RVF virus substitution rates accelerated in humans. 60.00 % RVF outbreaks occurred 0-2 months after floods and (or) heavy rainfall.

CONCLUSION

RVF has the increasing risk to cause pandemics, and global collaboration on "One Health" is needed to prevent potential pandemics.

Flavonoids derived from medicinal plants as a COVID-19 treatment

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes COVID-19 disease. Through its viral spike (S) protein, the virus enters and infects epithelial cells by utilizing angiotensin-converting enzyme 2 as a host cell's receptor protein. The COVID-19 pandemic had a profound impact on global public health and economies. Although various effective vaccinations and medications are now available to prevent and treat COVID-19, natural compounds derived from medicinal plants, particularly flavonoids, demonstrated therapeutic potential to treat COVID-19 disease. Flavonoids exhibit dual antiviral mechanisms: direct interference with viral invasion and inhibition of replication. Specifically, they target key viral molecules, particularly viral proteases, involved in infection. These compounds showcase significant immunomodulatory and anti-inflammatory properties, effectively inhibiting various inflammatory cytokines. Additionally, emerging evidence supports the potential of flavonoids to mitigate the progression of COVID-19 in individuals with obesity by positively influencing lipid metabolism. This review aims to elucidate the molecular structure of SARS-CoV-2 and the underlying mechanism of action of flavonoids on the virus. This study evaluates the potential anti-SARS-CoV-2 properties exhibited by flavonoid compounds, with a specific interest in their structure and mechanisms of action, as therapeutic applications for the prevention and treatment of COVID-19. Nevertheless, a significant portion of existing knowledge is based on theoretical frameworks and findings derived from in vitro investigations. Further research is required to better assess the effectiveness of flavonoids in combating SARS-CoV-2, with a particular emphasis on in vivo and clinical investigations.

Genetic Variation of Turnip Yellows Virus in Arable and Vegetable Brassica Crops, Perennial Wild Brassicas, and Aphid Vectors Collected from the Plants

Turnip yellows virus (TuYV; Polerovirus, Solemoviridae) infects and causes yield losses in a range of economically important crop species, particularly the Brassicaceae. It is persistently transmitted by several aphid species and is difficult to control. Although the incidence and genetic diversity of TuYV has been extensively investigated in recent years, little is known about how the diversity within host plants relates to that in its vectors. Arable oilseed rape (Brassica napus) and vegetable brassica plants (Brassica oleracea), wild cabbage (B. oleracea), and aphids present on these plants were sampled in the field in three regions of the United Kingdom. High levels of TuYV (82 to 97%) were detected in plants in all three regions following enzyme-linked immunosorbent assays. TuYV was detected by reverse transcription polymerase chain reaction in Brevicoryne brassicae aphids collected from plants, and TuYV sequences were obtained. Two TuYV open reading frames, ORF0 and ORF3, were partially sequenced from 15 plants, and from one aphid collected from each plant. Comparative analyses between TuYV sequences from host plants and B. brassicae collected from respective plants revealed differences between some ORF0 sequences, which possibly indicated that at least two of the aphids might not have been carrying the same TuYV isolates as those present in their host plants. Maximum likelihood phylogenetic analyses including published, the new TuYV sequences described above, 101 previously unpublished sequences of TuYV from oilseed rape in the United Kingdom, and 13 also previously unpublished sequences of TuYV from oilseed rape in Europe and China revealed three distinct major clades for ORF0 and one for ORF3, with some distinct subclades. Some clustering was related to geographic origin. Explanations for TuYV sequence differences between plants and the aphids present on respective plants and implications for the epidemiology and control of TuYV are discussed.

Portable, single nucleotide polymorphism-specific duplex assay for virus surveillance in wastewater

The Argonaute protein from the archaeon Pyrococcus furiosus (PfAgo) is a DNA-guided nuclease that targets DNA with any sequence. We designed a virus detection assay in which the PfAgo enzyme cleaves the reporter probe, thus generating fluorescent signals when amplicons from a reverse transcriptase loop-mediated isothermal amplification (RT-LAMP) assay contain target sequences. We confirmed that the RT-LAMP-PfAgo assay for the SARS-CoV-2 Delta variant produced significantly higher fluorescent signals (p < 0.001) when a single nucleotide polymorphism (SNP), exclusive to the Delta variant, was present, compared to the samples without the SNP. Additionally, the duplex assay for Pepper mild mottle virus (PMMOV) and SARS-CoV-2 detection produced specific fluorescent signals (FAM or ROX) only when the corresponding sequences were present. Furthermore, the RT-LAMP-PfAgo assay does not require dilution to reduce the impact of environmental inhibitors. The limit of detection of the PMMOV assay, determined with 30 wastewater samples, was 28 gc/μL, with a 95 % confidence interval of [11,103]. Finally, using a point-of-use device, the RT-LAMP-PfAgo assay successfully detected PMMOV in wastewater samples. Based on our findings, we conclude that the RT-LAMP-PfAgo assay can be used as a portable, SNP-specific duplex assay, which will significantly improve virus surveillance in wastewater.

SARS-CoV-2 and innate immunity: the good, the bad, and the "goldilocks"

An ancient conflict between hosts and pathogens has driven the innate and adaptive arms of immunity. Knowledge about this interplay can not only help us identify biological mechanisms but also reveal pathogen vulnerabilities that can be leveraged therapeutically. The humoral response to SARS-CoV-2 infection has been the focus of intense research, and the role of the innate immune system has received significantly less attention. Here, we review current knowledge of the innate immune response to SARS-CoV-2 infection and the various means SARS-CoV-2 employs to evade innate defense systems. We also consider the role of innate immunity in SARS-CoV-2 vaccines and in the phenomenon of long COVID.

Virulence evaluation of Israeli Marek's disease virus isolates from commercial poultry using their meq gene sequence

Fifty-seven Gallid alphaherpesvirus 2 (GaHV-2) isolates, collected during a 30-year period (1990-2019) from commercial poultry flocks affected by Marek's disease (MD), were molecularly characterised. The GaHV-2 meq gene was amplified and sequenced to evaluate the virus virulence, based on the number of PPPPs within the proline-rich repeats (PRRs) of its transactivation domain. The present illustration of virus virulence evaluation on a large scale of field virus isolates by molecular analysis exemplifies the practical benefit and usefulness of the molecular marker in commercial GaVH-2 isolates. The alternative assay of GaVH-2 virulence pathotyping is the classical Gold Standard ADOL method, which is difficult and impossible to employ on a large scale using the Specific Pathogen Free (SPF) chicks of the ADOL strains kept in isolators for two months. The phylogenetic analysis performed in the present study showed that the meq gene amino acid sequences of the 57 Israeli strains divide into 16 phylogenetic branches. The virulence evaluation was performed in comparison with 36 GaHV-2 prototype strains, previously characterised by the in vivo Gold Standard ADOL assay. The results obtained revealed that the GaHV-2 strains circulating in Israel have evolved into a higher virulence potential during the years, as the four-proline stretches number in the meq gene decreased over the investigated period, typically of very virulent virus prototypes. The present study supports the meq gene molecular markers for the assessment of field GaVH-2 strains virulence.

Machine learning-based approach KEVOLVE efficiently identifies SARS-CoV-2 variant-specific genomic signatures

Machine learning was shown to be effective at identifying distinctive genomic signatures among viral sequences. These signatures are defined as pervasive motifs in the viral genome that allow discrimination between species or variants. In the context of SARS-CoV-2, the identification of these signatures can assist in taxonomic and phylogenetic studies, improve in the recognition and definition of emerging variants, and aid in the characterization of functional properties of polymorphic gene products. In this paper, we assess KEVOLVE, an approach based on a genetic algorithm with a machine-learning kernel, to identify multiple genomic signatures based on minimal sets of k-mers. In a comparative study, in which we analyzed large SARS-CoV-2 genome dataset, KEVOLVE was more effective at identifying variant-discriminative signatures than several gold-standard statistical tools. Subsequently, these signatures were characterized using a new extension of KEVOLVE (KANALYZER) to highlight variations of the discriminative signatures among different classes of variants, their genomic location, and the mutations involved. The majority of identified signatures were associated with known mutations among the different variants, in terms of functional and pathological impact based on available literature. Here we showed that KEVOLVE is a robust machine learning approach to identify discriminative signatures among SARS-CoV-2 variants, which are frequently also biologically relevant, while bypassing multiple sequence alignments. The source code of the method and additional resources are available at: https://github.com/bioinfoUQAM/KEVOLVE.

Isolation, characterization, and phylogenetic analysis of two new porcine parvovirus 1 isolates from Northern China

Porcine parvovirus (PPV) is a pathogen of infectious reproductive disease, which can cause stillbirth, mummification, embryo death, and infertility (SMEDI) syndrome in pigs. The objective of this study was to gain new insights into the evolution and phylogeny of the PPV1 genome. In this study, we isolated two new PPV1 (HLJ202108-Y and SDLC202109) from northern China and sequenced their whole genomes. The new isolates were found to have three amino acid substitutions (K195R, K562R, and S578P) in nonstructural protein 1. The VP2 amino acid site contained nine nonsynonymous substitutions, including six substitutions of the Kresse strain corresponding to the NADL-2 strain and three substitutions of A414S, S436T, and N555K. Genetic evolution analysis was conducted on 107 reference sequences available in the GenBank database, and 4-5 PPV1 taxa were defined. The new isolates were in the same phylogenetic cluster as strain 27a. The changes in the cluster, specifically marker amino acids, and their potential role in enhancing pathogenicity are discussed in this study. Furthermore, the evolutionary tree map results showed that the strains in China were evolving in two directions: one was becoming increasingly similar to early NADL-2 strains, while the other was evolving toward 27a-like strains. We also compared the proliferation ability of the isolated strains in susceptible cells by analyzing the multistep growth curves. The results showed that the virulence titer of the mutant strain was high. In summary, this study introduced the latest changes in PPV and discussed the virus characteristics that were considered to affect virulence.

Generation of SARS-CoV-2 spike receptor binding domain mutants and functional screening for immune evaders using a novel lentivirus-based system

The emergence of rapid and continuous mutations of severe acute respiratory syndrome 2 (SARS-CoV-2) spike glycoprotein that increased with the Omicron variant points out the necessity to anticipate such mutations for conceiving specific and adaptable therapies to avoid another pandemic. The crucial target for the antibody treatment and vaccine design is the receptor binding domain (RBD) of the SARS-CoV-2 spike. It is also the site where the virus has shown its high ability to mutate and consequently escape immune response. We developed a robust and simple method for generating a large number of functional SARS-CoV-2 spike RBD mutants by error-prone PCR and a novel nonreplicative lentivirus-based system. We prepared anti-RBD wild type (WT) polyclonal antibodies and used them to screen and select for mutant libraries that escape inhibition of virion entry into recipient cells expressing human angiotensin-converting enzyme 2 and transmembrane serine protease 2. We isolated, cloned, and sequenced six mutants totally bearing nine mutation sites. Eight mutations were found in successive WT variants, including Omicron and other recombinants, whereas one is novel. These results, together with the detailed functional analyses of two mutants provided the proof of concept for our approach.

Quantification of Virion-Sense and Complementary-Sense DNA Strands of Circular Single-Stranded DNA Viruses

Circular ssDNA viruses are ubiquitous and can be found in both prokaryotes and eukaryotes. To understand the interaction of ssDNA viruses with their hosts, it is important to characterize the dynamics of viral sense (VS) and complementary-sense (CS) viral strands during the infection process. Here, we present a simple and rapid protocol that allows sensitive and accurate determination of the VS and CS strands generated during viral infection.The method consists of a two-step qPCR in which the first step uses a strand-specific (CS or VS) labeled primer and T4 DNA polymerase that lacks strand displacement activity and makes a single copy per VS or CS strand. Next, the T4 DNA polymerase and unincorporated oligonucleotides are removed by a silica membrane spin column. Finally, the purified VS or CS strands are quantified by qPCR in a second step in which amplification uses a tag primer and a specific primer. Absolute quantification of VS and CS strands is obtained by extrapolating the Cq data to a standard curve of ssDNA, which can be generated by phagemid expression. Quantification of VS and CS strands of two geminiviruses in infections of Solanum lycopersicum (tomato) and Nicotiana benthamiana plants using this method is shown.

Bioinformatic Approaches for Comparative Analysis of Viruses

The field of viral genomic studies has experienced an unprecedented increase in data volume. New strains of known viruses are constantly being added to the GenBank database and so are completely new species with little or no resemblance to our databases of sequences. In addition to this, metagenomic techniques have the potential to further increase the number and rate of sequenced genomes. Besides, it is important to consider that viruses have a set of unique features that often break down molecular biology dogmas, e.g., the flux of information from RNA to DNA in retroviruses and the use of RNA molecules as genomes. As a result, extracting meaningful information from viral genomes remains a challenge and standard methods for comparing the unknown and our databases of characterized sequences may need adaptations. Thus, several bioinformatic approaches and tools have been created to address the challenge of analyzing viral data. This chapter offers descriptions and protocols of some of the most important bioinformatic techniques for comparative analysis of viruses. The authors also provide comments and discussion on how viruses' unique features can affect standard analyses and how to overcome some of the major sources of problems. Protocols and topics emphasize online tools (which are more accessible to users) and give the real experience of what most bioinformaticians do in day-by-day work with command-line pipelines. The topics discussed include (1) clustering related genomes, (2) whole genome multiple sequence alignments for small RNA viruses, (3) protein alignment for marker genes and species affiliation, (4) variant calling and annotation, and (5) virome analyses and pathogen identification.

Taxonomic Classification of Geminiviruses Based on Pairwise Sequence Comparisons

Geminiviridae is the largest and one of the most diverse families of plant viruses, comprising 14 genera demarcated based on host range, type of insect vector, and phylogenetic relationships. The use of unbiased, whole-genome multiple displacement amplification techniques coupled with high-throughput sequencing has greatly expanded our knowledge of geminivirus diversity over the last two decades. As a result, a large number of new species have been described in recent years. Species demarcation criteria in the family are entirely based on sequence comparisons, but the specific cutoff values vary for each genus. The purpose of this chapter is to provide a step-by-step pipeline to classify new species in the family Geminiviridae.

Genotyping and phylogeography of infectious bronchitis virus isolates from Pakistan show unique linkage to GI-24 lineage

Infectious bronchitis virus (IBV) is prevalent in Pakistan causing enormous economic losses. To date no clear data are available on circulating genotypes and phylogeographic spread of the virus. Hence current study assessed these parameters for all available IBV Pakistani isolates, based on the 9 new sequences, with respect to other Asian and non-Asian countries. Results indicated that all Pakistani isolates belonged to genotype I (GI), with more than half of them (16/27) belonging to the GI-24 lineage, against which no vaccine is available. Three possible introduction events of the GI-13 IBV lineage into Pakistan, based on the estimated IBV population using isolates from this study, were observed possibly from Afghanistan, China, and/or Egypt. These events were further analyzed on the S1 amino acid level which showed unique alterations (S250H, T270K, and Q298S) in 1 isolate (IBV4, GI-13) when compared to GI-1 lineage. Both GI-1 and GI-13 Pakistani strains showed close homology with homologous vaccine strains that are used in Pakistan. For GI-24 strains, none of the used vaccines showed substantial homology, necessitating the need for further exploration of this lineage and vaccine design. In addition, our findings highlight the importance of genomic surveillance to support phylogeographical studies on IBV in genotyping and molecular epidemiology.

Parallel evolution and enhanced virulence upon in vivo passage of an RNA virus in Drosophila melanogaster

Virus evolution is strongly affected by antagonistic co-evolution of virus and host. Host immunity positively selects for viruses that evade the immune response, which in turn may drive counter-adaptations in host immune genes. We investigated how host immune pressure shapes virus populations, using the fruit fly Drosophila melanogaster and its natural pathogen Drosophila C virus (DCV), as a model. We performed an experimental evolution study in which DCV was serially passaged for ten generations in three fly genotypes differing in their antiviral RNAi response: wild-type flies and flies in which the endonuclease gene Dicer-2 was either overexpressed or inactivated. All evolved virus populations replicated more efficiently in vivo and were more virulent than the parental stock. The number of polymorphisms increased in all three host genotypes with passage number, which was most pronounced in Dicer-2 knockout flies. Mutational analysis showed strong parallel evolution, as mutations accumulated in a specific region of the VP3 capsid protein in every lineage in a host genotype-independent manner. The parental tyrosine at position ninety-five of VP3 was substituted with either one of five different amino acids in fourteen out of fifteen lineages. However, no consistent amino acid changes were observed in the viral RNAi suppressor gene 1A, nor elsewhere in the genome in any of the host backgrounds. Our study indicates that the RNAi response restricts the sequence space that can be explored by viral populations. Moreover, our study illustrates how evolution towards higher virulence can be a highly reproducible, yet unpredictable process.

Flexible changes to the Heliothis virescens ascovirus 3h (HvAV-3h) virion components affect pathogenicity against different host larvae species

IMPORTANCE

Different pathogenic processes of a virus in different hosts are related to the host individual differences, which makes the virus undergoes different survival pressures. Here, we found that the virions of an insect virus, Heliothis virescens ascovirus 3h (HvAV-3h), had different protein composition when they were purified from different host larval species. These "adaptive changes" of the virions were analyzed in detail in this study, which mainly included the differences of the protein composition of virions and the differences in affinity between virions and different host proteins. The results of this study revealed the flexible changes of viruses to help themselves adapt to different hosts. Also, these interesting findings can provide new insights to improve our understanding of virus adaptability and virulence differentiation caused by the adaptation process.

The Novel Porcine Parvoviruses: Current State of Knowledge and Their Possible Implications in Clinical Syndromes in Pigs

Parvoviruses (PVs) affect various animal species causing different diseases. To date, eight different porcine parvoviruses (PPV1 through PPV8) are recognized in the swine population, all of which are distributed among subfamilies and genera of the Parvoviridae family. PPV1 is the oldest and is recognized as the primary agent of SMEDI, while the rest of the PPVs (PPV2 through PPV8) are called novel PPVs (nPPVs). The pathogenesis of nPPVs is still undefined, and whether these viruses are putative disease agents is unknown. Structurally, the PPVs are very similar; the differences occur mainly at the level of their genomes (ssDNA), where there is variation in the number and location of the coding genes. Additionally, it is considered that the genome of PVs has mutation rates similar to those of ssRNA viruses, that is, in the order of 10-5-10-4 nucleotide/substitution/year. These mutations manifest mainly in the VP protein, constituting the viral capsid, affecting virulence, tropism, and viral antigenicity. For nPPVs, mutation rates have already been established that are similar to those already described; however, within this group of viruses, the highest mutation rate has been reported for PPV7. In addition to the mutations, recombinations are also reported, mainly in PPV2, PPV3, and PPV7; these have been found between strains of domestic pigs and wild boars and in a more significant proportion in VP sequences. Regarding affinity for cell types, nPPVs have been detected with variable prevalence in different types of organs and tissues; this has led to the suggestion that they have a broad tropism, although proportionally more have been found in lung and lymphoid tissue such as spleen, tonsils, and lymph nodes. Regarding their epidemiology, nPPVs are present on all continents (except PPV8, only in Asia), and within pig farms, the highest prevalences detecting viral genomes have been seen in the fattener and finishing groups. The relationship between nPPVs and clinical manifestations has been complicated to establish. However, there is already some evidence that establishes associations. One of them is PPV2 with porcine respiratory disease complex (PRDC), where causality tests (PCR, ISH, and histopathology) lead to proposing the PPV2 virus as a possible agent involved in this syndrome. With the other nPPVs, there is still no clear association with any pathology. These have been detected in different systems (respiratory, reproductive, gastrointestinal, urinary, and nervous), and there is still insufficient evidence to classify them as disease-causing agents. In this regard, nPPVs (except PPV8) have been found to cause porcine reproductive failure (PRF), with the most prevalent being PPV4, PPV6, and PPV7. In the case of PRDC, nPPVs have also been detected, with PPV2 having the highest viral loads in the lungs of affected pigs. Regarding coinfections, nPPVs have been detected in concurrence in healthy and sick pigs, with primary PRDC and PRF viruses such as PCV2, PCV3, and PRRSV. The effect of these coinfections is not apparent; it is unknown whether they favor the replication of the primary agents, the severity of the clinical manifestations, or have no effect. The most significant limitation in the study of nPPVs is that their isolation has been impossible; therefore, there are no studies on their pathogenesis both in vitro and in vivo. For all of the above, it is necessary to propose basic and applied research on nPPVs to establish if they are putative disease agents, establish their effect on coinfections, and measure their impact on swine production.

Equine dermatitis outbreak associated with parapoxvirus

Parapoxviruses (PPV) cause skin and mucous membrane lesions in several animal species, and of the five recognized PPVs, at least three are zoonotic. Equine PPV (EqPPV) is the sixth one initially described in humans in the United States and later in a severely sick horse in Finland in 2013-2015. In 2021-2022, a large-scale pustulo-vesicular pastern dermatitis outbreak occurred in horses all over Finland. This study aimed at analysing the outbreak, identifying and describing the causative agent, describing clinical signs, and searching for risk factors. EqPPV was identified as a probable causative agent and co-infections with several potentially pathogenic and zoonotic bacteria were observed. Histopathologically, suppurative and ulcerative dermatitis was diagnosed. Due to the lack of specific tests for this virus, we developed a novel diagnostic EqPPV-PCR with sensitivity of 10 copies/reaction. Based on a large proportion of the genome sequenced directly from clinical samples, very little variation was detected between the sequences of the case from 2013 and the cases from 2021 to 2022. Based on an epidemiological survey, the main risk factor for pastern dermatitis was having racehorses. Approximately one third of the horses at each affected stable got clinical dermatitis, manifesting as severe skin lesions. Skin lesions were also occasionally reported in humans, indicating potential zoonotic transmission. Case stables commonly reported attendance at race events before acquiring the disease. Survey also identified differences in practises between case and control stables. Taken together, these results enable a better preparedness, diagnostics, and guidelines for future outbreaks.

African swine fever virus - variants on the rise

African swine fever virus (ASFV), a large and complex DNA-virus circulating between soft ticks and indigenous suids in sub-Saharan Africa, has made its way into swine populations from Europe to Asia. This virus, causing a severe haemorrhagic disease (African swine fever) with very high lethality rates in wild boar and domestic pigs, has demonstrated a remarkably high genetic stability for over 10 years. Consequently, analyses into virus evolution and molecular epidemiology often struggled to provide the genetic basis to trace outbreaks while few resources have been dedicated to genomic surveillance on whole-genome level. During its recent incursion into Germany in 2020, ASFV has unexpectedly diverged into five clearly distinguishable linages with at least ten different variants characterized by high-impact mutations never identified before. Noticeably, all new variants share a frameshift mutation in the 3' end of the DNA polymerase PolX gene O174L, suggesting a causative role as possible mutator gene. Although epidemiological modelling supported the influence of increased mutation rates, it remains unknown how fast virus evolution might progress under these circumstances. Moreover, a tailored Sanger sequencing approach allowed us, for the first time, to trace variants with genomic epidemiology to regional clusters. In conclusion, our findings suggest that this new factor has the potential to dramatically influence the course of the ASFV pandemic with unknown outcome. Therefore, our work highlights the importance of genomic surveillance of ASFV on whole-genome level, the need for high-quality sequences and calls for a closer monitoring of future phenotypic changes of ASFV.

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.

A comprehensive review on the global efforts on vaccines and repurposed drugs for combating COVID-19

The recently discovered coronavirus, known as SARS-CoV-2, is a highly contagious and potentially lethal viral infection that was declared a pandemic by the World Health Organization on March 11, 2020. Since the beginning of the pandemic, an unprecedented number of COVID-19 vaccine candidates have been investigated for their potential to manage the pandemic. Herein, we reviewed vaccine development and the associated research effort, both traditional and forward-looking, to demonstrate the advantages and disadvantages of their technology, in addition to their efficacy limitations against mutant SARS-CoV-2. Moreover, we report repurposed drug discovery, which mainly focuses on virus-based and host-based targets, as well as their inhibitors. SARS-CoV-2 targets include the main protease (Mpro), and RNA-dependent RNA-polymerase (RdRp), which are the most well-studied and conserved across coronaviruses, enabling the development of broad-spectrum inhibitors of these enzymes.

Spike Gene Analysis and Prevalence of Porcine Epidemic Diarrhea Virus from Pigs in South Korea: 2013-2022

From late 2013-2022, 1131 cases of porcine epidemic diarrhea (PED) were reported to the Korean Animal Health Integrated System (KAHIS). There were four major outbreaks from winter to spring (2013-2014, 2017-2018, 2018-2019, and 2021-2022), with the main outbreaks occurring in Chungnam (CN), Jeonbuk (JB), and Jeju (JJ). Analysis of the complete spike (S) gene of 140/1131 KAHIS PEDV cases nationwide confirmed that 139 belonged to the G2b genotype and 1 to the G2a genotype. Among them, two strains (K17GG1 and K17GB3) were similar to an S INDEL isolated in the United States (strain OH851), and 12 strains had deletions (nucleotides (nt) 3-99) or insertions (12 nt) within the S gene. PEDVs in JJ formed a regionally independent cluster. The substitution rates (substitutions/site/year) were as follows: 1.5952 × 10-3 in CN, 1.8065 × 10-3 in JB, and 1.5113 × 10-3 in JJ. A Bayesian skyline plot showed that the effective population size of PEDs in JJ fell from 2013-2022, whereas in CN and JB it was maintained. Genotyping of 340 Korean PEDV strains, including the 140 PEDVs in this study and 200 Korean reference strains from GenBank, revealed that only the highly pathogenic non-INDEL type (G2b) was dominant from 2020 onwards. Therefore, it is predicted that the incidence of PED will be maintained by the G2b (non-INDEL) genotype.

Molecular Epidemiology of Rotavirus A in Calves: Evolutionary Analysis of a Bovine G8P[11] Strain and Spatio-Temporal Dynamics of G6 Lineages in the Americas

Rotavirus A (RVA) causes diarrhea in calves and frequently possesses the G6 and P[5]/P[11] genotypes, whereas G8 is less common. We aimed to compare RVA infections and G/P genotypes in beef and dairy calves from major livestock regions of Argentina, elucidate the evolutionary origin of a G8 strain and analyze the G8 lineages, infer the phylogenetic relationship of RVA field strains, and investigate the evolution and spatio-temporal dynamics of the main G6 lineages in American countries. Fecal samples (n = 422) from diarrheic (beef, 104; dairy, 137) and non-diarrheic (beef, 78; dairy, 103) calves were analyzed by ELISA and semi-nested multiplex RT-PCR. Sequencing, phylogenetic, phylodynamic, and phylogeographic analyses were performed. RVA infections were more frequent in beef (22.0%) than in dairy (14.2%) calves. Prevalent genotypes and G6 lineages were G6(IV)P[5] in beef (90.9%) and G6(III)P[11] (41.2%) or mixed genotypes (23.5%) in dairy calves. The only G8 strain was phylogenetically related to bovine and artiodactyl bovine-like strains. Re-analyses inside the G8 genotype identified G8(I) to G8(VIII) lineages. Of all G6 strains characterized, the G6(IV)P[5](I) strains from "Cuenca del Salado" (Argentina) and Uruguay clustered together. According to farm location, a clustering pattern for G6(IV)P[5] strains of beef farms was observed. Both G6 lineage strains together revealed an evolutionary rate of 1.24 × 10-3 substitutions/site/year, and the time to the most recent common ancestor was dated in 1853. The most probable ancestral locations were Argentina in 1981 for G6(III) strains and the USA in 1940 for G6(IV) strains. The highest migration rates for both G6 lineages together were from Argentina to Brazil and Uruguay. Altogether, the epidemiology, genetic diversity, and phylogeny of RVA in calves can differ according to the production system and farm location. We provide novel knowledge about the evolutionary origin of a bovine G8P[11] strain. Finally, bovine G6 strains from American countries would have originated in the USA nearly a century before its first description.

Cobalt-assisted route to rare carbocyclic C-ribonucleosides

(Re)emerging RNA viruses have been major threats to public health in the past years, and from the few drugs available, nucleoside analogues are still at the cornerstone of the antiviral therapy. Among them, the synthesis of carbocyclic C-nucleosides is suffering from long syntheses and poor yields. Herein we report a concise stereoselective synthesis of rare carbocyclic C-nucleosides (11a-l) bearing non-canonical nucleobases through a cobalt-assisted-route as key step starting from the optically pure (-)-cyclopentenone 1. This approach paves the route for novel carbocyclic C-nucleoside discovery.

Distinct pathways of adaptive evolution in Cryptococcus neoformans reveal a mutation in adenylyl cyclase with trade-offs for pathogenicity

Pathogenic fungi populate a wide range of environments and infect a diversity of host species. Despite this substantial biological flexibility, the impact of interactions between fungi and their hosts on the evolution of pathogenicity remains unclear. We studied how repeated interactions between the fungus Cryptococcus neoformans and relevant environmental and mammalian host cells-amoeba and mouse macrophages-shape the evolution of this model fungal pathogen. First, using a collection of clinical and environmental isolates of C. neoformans, we characterized a range of survival phenotypes for these strains when exposed to host cells of different species. We then performed serial passages of an environmentally isolated C. neoformans strain through either amoeba or macrophages for ∼75 generations to observe how these interactions select for improved replication within hosts. In one adapted population, we identified a single point mutation in the adenylyl cyclase gene, CAC1, that swept to fixation and confers a strong competitive advantage for growth inside macrophages. Strikingly, this growth advantage in macrophages is inversely correlated with disease severity during mouse infections, suggesting that adaptation to specific host niches can markedly reduce the pathogenicity of these fungi. These results raise intriguing questions about the influence of cyclic AMP (cAMP) signaling on pathogenicity and highlight the role of seemingly small adaptive changes in promoting fundamental shifts in the intracellular behavior and virulence of these important human pathogens.

A Proofreading Mutation with an Allosteric Effect Allows a Cluster of SARS-CoV-2 Viruses to Rapidly Evolve

The RNA-dependent RNA polymerase of the severe acute respiratory syndrome coronavirus 2 virus is error prone, with errors being corrected by the exonuclease (NSP14) proofreading mechanism. However, the mutagenesis and subsequent evolutionary trajectory of the virus is mediated by the delicate interplay of replicase fidelity and environmental pressures. Here, we have shown that a single, distal mutation (F60S) in NSP14 can have a profound impact upon proofreading with an increased accumulation of mutations and elevated evolutionary rate being observed. Understanding the implications of these changes is crucial, as these underlying mutational processes may have important implications for understanding the population-wide evolution of the virus. This study underscores the urgent need for continued research into the replicative mechanisms of this virus to combat its continued impact on global health, through the re-emergence of immuno-evasive variants.

New Potential Role of European Mouflon (Ovis aries musimon) in the Epidemiology of Border Disease in the Pyrenees

Since 2001, high-mortality outbreaks of border disease (BD) have negatively affected populations of Pyrenean chamois (Rupicapra pyrenaica pyrenaica). Studies in the affected areas determined that sympatric wild ruminants did not seem to have an epidemiologic role in the circulation of border disease virus (BDV). However, the recent increase in European mouflon (Ovis aries musimon) densities might enhance the risk of pathogen transmission among chamois and mouflons. We conducted a serologic and virologic investigation of BDV in European mouflon from the Spanish Pyrenees, with the aim of determining potential changes in the role of this species in BDV epidemiology. From 2018 to 2022, we detected antibodies against BDV in 31/185 (16.7%) animals but did not detect BDV RNA in any spleen sample (0/65). These results indicate that BDV infection is occurring in these mouflon populations to a greater extent than previously described, which could shift the current understanding of BD epidemiology in the Pyrenees and cause an unpredictable effect on both chamois and mouflon populations. Further studies on the molecular identification of BDV in mouflon and chamois are required to better understand the contribution of mouflon in the epidemiology of BD.

Viral Prevalence and Genomic Xenology in the Coevolution of HzNV-2 (Nudiviridae) with Host Helicoverpa zea (Lepidoptera: Noctuidae)

Insect viruses have been described from numerous lineages, yet patterns of genetic exchange and viral prevalence, which are essential to understanding host-virus coevolution, are rarely studied. In Helicoverpa zea, the virus HzNV-2 can cause deformity of male and female genitalia, resulting in sterility. Using ddPCR, we found that male H. zea with malformed genitalia (agonadal) contained high levels of HzNV-2 DNA, confirming previous work. HzNV-2 was found to be prevalent throughout the United States, at more than twice the rate of the baculovirus HaSNPV, and that it contained several host-acquired DNA sequences. HzNV-2 possesses four recently endogenized lepidopteran genes and several more distantly related genes, including one gene with a bacteria-like sequence found in both host and virus. Among the recently acquired genes is cytosolic serine hydroxymethyltransferase (cSHMT). In nearly all tested H. zea, cSHMT contained a 200 bp transposable element (TE) that was not found in cSHMT of the sister species H. armigera. No other virus has been found with host cSHMT, and the study of this shared copy, including possible interactions, may yield new insights into the function of this gene with possible applications to insect biological control, and gene editing.

Viral Evolution Shaped by Host Proteostasis Networks

Understanding the factors that shape viral evolution is critical for developing effective antiviral strategies, accurately predicting viral evolution, and preventing pandemics. One fundamental determinant of viral evolution is the interplay between viral protein biophysics and the host machineries that regulate protein folding and quality control. Most adaptive mutations in viruses are biophysically deleterious, resulting in a viral protein product with folding defects. In cells, protein folding is assisted by a dynamic system of chaperones and quality control processes known as the proteostasis network. Host proteostasis networks can determine the fates of viral proteins with biophysical defects, either by assisting with folding or by targeting them for degradation. In this review, we discuss and analyze new discoveries revealing that host proteostasis factors can profoundly shape the sequence space accessible to evolving viral proteins. We also discuss the many opportunities for research progress proffered by the proteostasis perspective on viral evolution and adaptation.

Trends and Challenges in the Surveillance and Control of Avian Metapneumovirus

Among the respiratory pathogens of birds, the Avian Metapneumovirus (aMPV) is one of the most relevant, as it is responsible for causing infections of the upper respiratory tract and may induce respiratory syndromes. aMPV is capable of affecting the reproductive system of birds, directly impacting shell quality and decreasing egg production. Consequently, this infection can cause disorders related to animal welfare and zootechnical losses. The first cases of respiratory syndromes caused by aMPV were described in the 1970s, and today six subtypes (A, B, C, D, and two more new subtypes) have been identified and are widespread in all chicken and turkey-producing countries in the world, causing enormous economic losses for the poultry industry. Conventionally, immunological techniques are used to demonstrate aMPV infection in poultry, however, the identification of aMPV through molecular techniques helped in establishing the traceability of the virus. This review compiles data on the main aMPV subtypes present in different countries; aMPV and bacteria co-infection; vaccination against aMPV and viral selective pressure, highlighting the strategies used to prevent and control respiratory disease; and addresses tools for viral diagnosis and virus genome studies aiming at improving and streamlining pathogen detection and corroborating the development of new vaccines that can effectively protect herds, preventing viral escapes.

A-to-I RNA editing by ADAR and its therapeutic applications: From viral infections to cancer immunotherapy

ADAR deaminases catalyze adenosine-to-inosine (A-to-I) editing on double-stranded RNA (dsRNA) substrates that regulate an umbrella of biological processes. One of the two catalytically active ADAR enzymes, ADAR1, plays a major role in innate immune responses by suppression of RNA sensing pathways which are orchestrated through the ADAR1-dsRNA-MDA5 axis. Unedited immunogenic dsRNA substrates are potent ligands for the cellular sensor MDA5. Upon activation, MDA5 leads to the induction of interferons and expression of hundreds of interferon-stimulated genes with potent antiviral activity. In this way, ADAR1 acts as a gatekeeper of the RNA sensing pathway by striking a fine balance between innate antiviral responses and prevention of autoimmunity. Reduced editing of immunogenic dsRNA by ADAR1 is strongly linked to the development of common autoimmune and inflammatory diseases. In viral infections, ADAR1 exhibits both antiviral and proviral effects. This is modulated by both editing-dependent and editing-independent functions, such as PKR antagonism. Several A-to-I RNA editing events have been identified in viruses, including in the insidious viral pathogen, SARS-CoV-2 which regulates viral fitness and infectivity, and could play a role in shaping viral evolution. Furthermore, ADAR1 is an attractive target for immuno-oncology therapy. Overexpression of ADAR1 and increased dsRNA editing have been observed in several human cancers. Silencing ADAR1, especially in cancers that are refractory to immune checkpoint inhibitors, is a promising therapeutic strategy for cancer immunotherapy in conjunction with epigenetic therapy. The mechanistic understanding of dsRNA editing by ADAR1 and dsRNA sensing by MDA5 and PKR holds great potential for therapeutic applications. This article is categorized under: RNA Processing > RNA Editing and Modification RNA in Disease and Development > RNA in Disease.

Genomic evolution of bovine viral diarrhea virus based on complete genome and individual gene analyses

Bovine viral diarrhea virus (BVDV) genome consists of a single-stranded, positive-sense RNA with high genetic diversity. In the last years, significant progress has been achieved in BVDV knowledge evolution through phylodynamic analysis based on the partial 5'UTR sequences, whereas a few studies have used other genes or the complete coding sequence (CDS). However, no research has evaluated and compared BVDV evolutionary history based on the complete genome (CG), CDS, and individual genes. In this study, phylodynamic analyses were carried out with BVDV-1 (Pestivirus A) and BVDV-2 (Pestivirus B) CG sequences available on the GenBank database and each genomic region: CDS, UTRs, and individual genes. In comparison to the CG, the estimations for both BVDV species varied according to the dataset used, pointing out the importance of considering the analyzed genomic region when concluding. This study may provide new insight into BVDV evolution history while highlighting the need to increase the available BVDV CG sequences to perform more comprehensive phylodynamic studies in the future.

Optimized quantification of intra-host viral diversity in SARS-CoV-2 and influenza virus sequence data

High error rates of viral RNA-dependent RNA polymerases lead to diverse intra-host viral populations during infection. Errors made during replication that are not strongly deleterious to the virus can lead to the generation of minority variants. However, accurate detection of minority variants in viral sequence data is complicated by errors introduced during sample preparation and data analysis. We used synthetic RNA controls and simulated data to test seven variant-calling tools across a range of allele frequencies and simulated coverages. We show that choice of variant caller and use of replicate sequencing have the most significant impact on single-nucleotide variant (SNV) discovery and demonstrate how both allele frequency and coverage thresholds impact both false discovery and false-negative rates. When replicates are not available, using a combination of multiple callers with more stringent cutoffs is recommended. We use these parameters to find minority variants in sequencing data from SARS-CoV-2 clinical specimens and provide guidance for studies of intra-host viral diversity using either single replicate data or data from technical replicates. Our study provides a framework for rigorous assessment of technical factors that impact SNV identification in viral samples and establishes heuristics that will inform and improve future studies of intra-host variation, viral diversity, and viral evolution. IMPORTANCE When viruses replicate inside a host cell, the virus replication machinery makes mistakes. Over time, these mistakes create mutations that result in a diverse population of viruses inside the host. Mutations that are neither lethal to the virus nor strongly beneficial can lead to minority variants that are minor members of the virus population. However, preparing samples for sequencing can also introduce errors that resemble minority variants, resulting in the inclusion of false-positive data if not filtered correctly. In this study, we aimed to determine the best methods for identification and quantification of these minority variants by testing the performance of seven commonly used variant-calling tools. We used simulated and synthetic data to test their performance against a true set of variants and then used these studies to inform variant identification in data from SARS-CoV-2 clinical specimens. Together, analyses of our data provide extensive guidance for future studies of viral diversity and evolution.

Porcine Reproductive and Respiratory Syndrome (PRRSV2) Viral Diversity within a Farrow-to-Wean Farm Cohort Study

Describing PRRSV whole-genome viral diversity data over time within the host and within-farm is crucial for a better understanding of viral evolution and its implications. A cohort study was conducted at one naïve farrow-to-wean farm reporting a PRRSV outbreak. All piglets 3-5 days of age (DOA) born to mass-exposed sows through live virus inoculation with the recently introduced wild-type virus two weeks prior were sampled and followed up at 17-19 DOA. Samples from 127 piglets were individually tested for PRRSV by RT-PCR and 100 sequences were generated using Oxford Nanopore Technologies chemistry. Female piglets had significantly higher median Ct values than males (15.5 vs. 13.7, Kruskal-Wallis p < 0.001) at 3-5 DOA. A 52.8% mortality between sampling points was found, and the odds of dying by 17-19 DOA decreased with every one unit increase in Ct values at 3-5 DOA (OR = 0.76, 95% CI 0.61-0.94, p = 0.01). Although the within-pig percent nucleotide identity was overall high (99.7%) between 3-5 DOA and 17-19 DOA samples, ORFs 4 and 5a showed much lower identities (97.26% and 98.53%, respectively). When looking solely at ORF5, 62% of the sequences were identical to the 3-5 DOA consensus. Ten and eight regions showed increased nucleotide and amino acid genetic diversity, respectively, all found throughout ORFs 2a/2b, 4, 5a/5, 6, and 7.