Mechanisms and consequences of positive-strand RNA virus recombination

A review of the recombination events, mechanisms and consequences of Coxsackievirus A6

Hand, foot, and mouth disease (HFMD) is one of the most common class C infectious diseases, posing a serious threat to public health worldwide. Enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16) have been regarded as the major pathogenic agents of HFMD; however, since an outbreak caused by coxsackievirus A6 (CV-A6) in France in 2008, CV-A6 has gradually become the predominant pathogen in many regions. CV-A6 infects not only children but also adults, and causes atypical clinical symptoms such as a more generalized rash, eczema herpeticum, high fever, and onychomadesis, which are different from the symptoms associated with EV-A71 and CV-A16. Importantly, the rate of genetic recombination of CV-A6 is high, which can lead to changes in virulence and the rapid evolution of other characteristics, thus posing a serious threat to public health. To date, no specific vaccines or therapeutics have been approved for CV-A6 prevention or treatment, hence it is essential to fully understand the relationship between recombination and evolution of this virus. Here, we systematically review the genetic recombination events of CV-A6 that have occurred worldwide and explore how these events have promoted virus evolution, thus providing important information regarding future HFMD surveillance and prevention.

A novel VP1-based enzyme-linked immunosorbent assay revealed widespread Enterovirus G infections in Guangxi, China

Enterovirus G (EV-G) has recently been shown to affect weight gain and cause neurological symptoms in piglets. However, the serological investigation of EV-G is limited. In this study, we developed a novel serological detection method based on the structural protein, VP1 of EV-G. The intra-assay and inter-assay coefficient variations were 3.2-8.9% and 2.6-8.0%, respectively. There was no cross-reaction of the VP1-based enzyme-linked immunosorbent assay (ELISA) with antisera against the other known porcine viruses. In addition, a comparison was made with other methods including the developed indirect ELISAs based on VP2 and VP3 proteins and western blot (WB) analysis, which demonstrated the reliability of the novel method. Using the VP1-based ELISA, we carried out the first seroepidemiological survey of EV-G in China by testing 1041 serum samples collected from different pig farms in Guangxi from 2019 to 2021. Our results showed that 68.78% of the serum samples and 100% of the pig farms were positive for EV-G, with a relatively high incidence of seropositivity in pigs of different ages. This was specifically evident in fattening pigs and sows, which may suggest that the piglets have experienced an infection with EV-G during their growth process. Our data provide the first serological evidence of EV-G infections in pigs from China and reveal the widespread presence of EV-G infections in Guangxi, China.

Genetic Variation and Evolutionary Analysis of Eggplant Mottled Dwarf Virus Isolates from Spain

The genetic variation and population structure of gene N (nucleocapsid) and part of gene L (replicase) from 13 eggplant mottle dwarf virus (EMDV) isolates from Spain were evaluated and compared with sequences of EMDV isolates from other countries retrieved from GenBank. Phylogenetic inference of part of gene L showed three main clades, one containing an EMDV isolate from Australia and the other two containing isolates from Iran and Europe, as well as four subclades. EMDV isolates from Spain were genetically very similar and grouped in a subclade together with one isolate from Germany and one from the UK. No new recombination events were detected in addition to one recombination previously reported, suggesting that recombination is rare for EMDV. The comparison of synonymous and non-synonymous rates showed that negative selection played an important role, and only two codons were under positive selection. Genetic differentiation (Fst test), phylogenetic and nucleotide diversity analyses suggest a unique introduction of EMDV to Spain and low gene flow with other countries. In contrast, Greece and Italy showed diverse populations with high gene flow between both.

Recombination between coronaviruses and synthetic RNAs and biorisk implications motivated by a SARS-CoV-2 FCS origin controversy

The urgent need for improved policy, regulation, and oversight of research with potential pandemic pathogens (PPPs) has been widely acknowledged. A 2022 article in Frontiers in Virology raises questions, reporting on a 100% sequence homology between the SARS-CoV-2 furin cleavage site (FCS) and the negative strand of a 2017 patented sequence. Even though Ambati and collaborators suspect a possible inadvertent or intentional cause leading to the FCS insert, the related underpinnings have not been studied from the perspective of potential biorisk policy gaps. A commentary on their article contests the low coincidence likelihood that was calculated by Ambati et al., arguing that the sequence match could have been a chance occurrence alone. Additionally, it has been suggested that the odds of the recombination event may be low. These considerations seem to have put many speculations related to any implied viral beginnings, notably from a research setting likely outside the Wuhan Institute of Virology, to rest. However, potential implications for future disasters in terms of biosafety and biosecurity have not been addressed. To demonstrate the feasibility of the Ambati et al. postulate, a theoretical framework is developed that substantially extends the research orientations implicated by these authors and the related patent. It is argued that specific experimental conditions, in combination, could significantly increase the implied recombination profile between coronaviruses and synthetic RNAs. Consequently, this article scrutinizes these largely unrecognized vulnerabilities to discuss implications across the spectrum of the biological risk landscape, with special attention to a potential "crime harvest." Focusing on insufficiently understood features of interaction between the natural and man-made world, vulnerabilities related to contaminants, camouflaging, and various misuse potentials fostered by the digitization and computerization of synthetic biology, it highlights novel biorisk gaps not covered by existing PPP policy. Even though this work does not aim to provide proof of the viral origin, it will make the point that, in theory, a convergence of under-appreciated lab experiments and technologies could have led to the SARS-CoV-2 FCS insert, which analogously could be exploited by various threat actors for the clandestine genesis of similar or even worse pathogens.

Evaluation of Potential In Vitro Recombination Events in Codon Deoptimized FMDV Strains

Codon deoptimization (CD) has been recently used as a possible strategy to derive foot-and-mouth disease (FMD) live-attenuated vaccine (LAV) candidates containing DIVA markers. However, reversion to virulence, or loss of DIVA, from possible recombination with wild-type (WT) strains has yet to be analyzed. An in vitro assay was developed to quantitate the levels of recombination between WT and a prospective A24-P2P3 partially deoptimized LAV candidate. By using two genetically engineered non-infectious RNA templates, we demonstrate that recombination can occur within non-deoptimized viral genomic regions (i.e., 3′end of P3 region). The sequencing of single plaque recombinants revealed a variety of genome compositions, including full-length WT sequences at the consensus level and deoptimized sequences at the sub-consensus/consensus level within the 3′end of the P3 region. Notably, after further passage, two recombinants that contained deoptimized sequences evolved to WT. Overall, recombinants featuring large stretches of CD or DIVA markers were less fit than WT viruses. Our results indicate that the developed assay is a powerful tool to evaluate the recombination of FMDV genomes in vitro and should contribute to the improved design of FMDV codon deoptimized LAV candidates.

Patterns and Temporal Dynamics of Natural Recombination in Noroviruses

Noroviruses infect a wide range of mammals and are the major cause of gastroenteritis in humans. Recombination at the junction of ORF1 encoding nonstructural proteins and ORF2 encoding major capsid protein VP1 is a well-known feature of noroviruses. Using all available complete norovirus sequences, we systematically analyzed patterns of natural recombination in the genus Norovirus both throughout the genome and across the genogroups. Recombination events between nonstructural (ORF1) and structural genomic regions (ORF2 and ORF3) were found in all analyzed genogroups of noroviruses, although recombination was most prominent between members of GII, the most common genogroup that infects humans. The half-life times of recombinant forms (clades without evidence of recombination) of human GI and GII noroviruses were 10.4 and 8.4-11.3 years, respectively. There was evidence of many recent recombination events, and most noroviruses that differed by more than 18% of nucleotide sequence were recombinant relative to each other. However, there were no distinct recombination events between viruses that differed by over 42% in ORF2/3, consistent with the absence of systematic recombination between different genogroups. The few inter-genogroup recombination events most likely occurred between ancient viruses before they diverged into contemporary genogroups. The recombination events within ORF1 or between ORF2/3 were generally rare. Thus, noroviruses routinely exchange full structural and nonstructural blocks of the genome, providing a modular evolution.

Incipient Parallel Evolution of SARS-CoV-2 Deltacron Variant in South Brazil

With the coexistence of multiple lineages and increased international travel, recombination and gene flow are likely to become increasingly important in the adaptive evolution of SARS-CoV-2. These processes could result in genetic introgression and the incipient parallel evolution of multiple recombinant lineages. However, identifying recombinant lineages is challenging, and the true extent of recombinant evolution in SARS-CoV-2 may be underestimated. This study describes the first SARS-CoV-2 Deltacron recombinant case identified in Brazil. We demonstrate that the recombination breakpoint is at the beginning of the Spike gene. The 5' genome portion (circa 22 kb) resembles the AY.101 (Delta), and the 3' genome portion (circa 8 kb nucleotides) is most similar to the BA.1.1 (Omicron). Furthermore, evolutionary genomic analyses indicate that the new strain emerged after a single recombination event between lineages of diverse geographical locations in December 2021 in South Brazil. This Deltacron, AYBA-RS, is one of the dozens of recombinants described in 2022. The submission of only four sequences in the GISAID database suggests that this lineage had a minor epidemiological impact. However, the recent emergence of this and other Deltacron recombinant lineages (XD, XF, and XS) suggests that gene flow and recombination may play an increasingly important role in the COVID-19 pandemic. We explain the evolutionary and population genetic theory that supports this assertion, concluding that this stresses the need for continued genomic surveillance. This monitoring is vital for countries where multiple variants are present, as well as for countries that receive significant inbound international travel.

Evolution of cucurbit-infecting tobamoviruses: Recombination and codon usage bias

Currently, there are seven cucurbit-infecting tobamoviruses comprising cucumber green mottle mosaic virus (CGMMV), Kyuri green mottle mosaic virus (KGMMV), cucumber fruit mottle mosaic virus (CFMMV), zucchini green mottle mosaic virus (ZGMMV), cucumber mottle virus (CMoV), watermelon green mottle mosaic virus (WGMMV), and Trichosanthes mottle mosaic virus (TrMMV). To gain more insights into their evolution, recombination analyses were conducted. Four CGMMV isolates and one KGMMV isolate were suggested to be recombinants. And there was an interspecies recombination event between CGMMV and ZGMMV. Phylogenetic incongruence was also observed for CGMMV and KGMMV. A probable ancestral pattern was inferred for the gene junction region between RdRp and MP. Codon usage bias analysis revealed that the viral genes had additional influence independent of compositional constraint. In codon preference, the seven viruses were both similar to and different from the host cucumber (Cucumis sativus). Moreover, the viruses were not deficient in CpG and UpA dinucleotides.

Genetic diversity and phylogenetic characteristics of viruses in lily plants in Beijing

Lily (Lilium) is an important bulbous perennial herb that is frequently infected by one or more viruses. To investigate the diversity of lily viruses, lilies with virus-like symptoms in Beijing were collected to perform small RNA deep sequencing. Then, the 12 complete and six nearly full-length viral genomes, including six known viruses and two novel viruses were determined. Based on sequence and phylogenetic analyses, two novel viruses were considered to be members of the genera Alphaendornavirus (Endornaviridae) and Polerovirus (Solemoviridae). These two novel viruses were provisionally named lily-associated alphaendornavirus 1 (LaEV-1) and lily-associated polerovirus 1 (LaPV-1). Based on sequence, phylogenetic and recombination analyses, strawberry latent ringspot virus (SLRSV) in the genus Stralarivirus (Secoviridae) was identified for the first time in China, and shown to exhibit the highest nucleotide (nt) diversity among the available full-length SLRSV genome sequences, with the highest identities of 79.5% for RNA1 and 80.9% for RNA2. Interestingly, the protease cofactor region in RNA1 was 752 aa in length, whereas those of the other 27 characterized isolates ranged from 700-719 aa in length. The genome sequences of lily virus A (Potyvirus), lily virus X (Potexvirus), and plantago asiatica mosaic virus (Potexvirus) exhibited varying degrees of sequence diversity at the nucleotide level compared with their corresponding characterized isolates. In addition, plantago asiatica mosaic virus (PlAMV) tended to cluster on a host species-basis. One identified lily mottle virus (Potyvirus) isolate was detected as a recombinant, and which clustered in a different group with four other isolates. Seven identified lily symptomless virus (Carlavirus) isolates, including one recombinant, were clustered into three clades. Our results revealed the genetic diversity of lily-infecting viruses, and sequence insertion, host species and recombination are factors that likely contribute to this diversity. Collectively, our results provide useful information regarding the control of viral disease in lily.

A Glimpse on the Evolution of RNA Viruses: Implications and Lessons from SARS-CoV-2

RNA viruses are characterised by extremely high genetic variability due to fast replication, large population size, low fidelity, and (usually) a lack of proofreading mechanisms of RNA polymerases leading to high mutation rates. Furthermore, viral recombination and reassortment may act as a significant evolutionary force among viruses contributing to greater genetic diversity than obtainable by mutation alone. The above-mentioned properties allow for the rapid evolution of RNA viruses, which may result in difficulties in viral eradication, changes in virulence and pathogenicity, and lead to events such as cross-species transmissions, which are matters of great interest in the light of current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemics. In this review, we aim to explore the molecular mechanisms of the variability of viral RNA genomes, emphasising the evolutionary trajectory of SARS-CoV-2 and its variants. Furthermore, the causes and consequences of coronavirus variation are explored, along with theories on the origin of human coronaviruses and features of emergent RNA viruses in general. Finally, we summarise the current knowledge on the circulating variants of concern and highlight the many unknowns regarding SARS-CoV-2 pathogenesis.

Heterologous RNA Recombination in the Cystoviruses φ6 and φ8: A Mechanism of Viral Variation and Genome Repair

Recombination and mutation of viral genomes represent major mechanisms for viral evolution and, in many cases, moderate pathogenicity. Segmented genome viruses frequently undergo reassortment of the genome via multiple infection of host organisms, with influenza and reoviruses being well-known examples. Specifically, major genomic shifts mediated by reassortment are responsible for radical changes in the influenza antigenic determinants that can result in pandemics requiring rapid preventative responses by vaccine modifications. In contrast, smaller mutational changes brought about by the error-prone viral RNA polymerases that, for the most part, lack a replication base mispairing editing function produce small mutational changes in the RNA genome during replication. Referring again to the influenza example, the accumulated mutations-known as drift-require yearly vaccine updating and rapid worldwide distribution of each new formulation. Coronaviruses with a large positive-sense RNA genome have long been known to undergo intramolecular recombination likely mediated by copy choice of the RNA template by the viral RNA polymerase in addition to the polymerase-based mutations. The current SARS-CoV-2 origin debate underscores the importance of understanding the plasticity of viral genomes, particularly the mechanisms responsible for intramolecular recombination. This review describes the use of the cystovirus bacteriophage as an experimental model for recombination studies in a controlled manner, resulting in the development of a model for intramolecular RNA genome alterations. The review relates the sequence of experimental studies from the laboratory of Leonard Mindich, PhD at the Public Health Research Institute-then in New York City-and covers a period of approximately 12 years. Hence, this is a historical scientific review of research that has the greatest relevance to current studies of emerging RNA virus pathogens.

A 21L/BA.2-21K/BA.1 "MixOmicron" SARS-CoV-2 hybrid undetected by qPCR that screen for variant in routine diagnosis

Among the multiple SARS-CoV-2 variants identified since summer 2020, several have co-circulated, creating opportunities for coinfections and potentially genetic recombinations that are common in coronaviruses. Viral recombinants are indeed beginning to be reported more frequently. Here, we describe a new SARS-CoV-2 recombinant genome that is mostly that of a Omicron 21L/BA.2 variant but with a 3' tip originating from a Omicron 21K/BA.1 variant. Two such genomes were obtained in our institute from adults sampled in February 2022 in university hospitals of Marseille, southern France, by next-generation sequencing carried out with the Illumina or Nanopore technologies. The recombination site was located between nucleotides 26,858-27,382. In the two genomic assemblies, mean sequencing depth at mutation-harboring positions was 271 and 1362 reads and mean prevalence of the majoritary nucleotide was 99.3 ± 2.2% and 98.8 ± 1.6%, respectively. Phylogeny generated trees with slightly different topologies according to whether genomes analyzed were depleted or not of the 3' tip. This 3' terminal end brought in the Omicron 21L/BA.2 genome a short transposable element of 41 nucleotides named S2m that is present in most SARS-CoV-2 except a few variants among which the Omicron 21L/BA.2 variant and may be involved in virulence. Importantly, this recombinant is not detected by currently used qPCR that screen for variants in routine diagnosis. The present observation emphasizes the need to survey closely the genetic pathways of SARS-CoV-2 variability by whole genome sequencing, and it could contribute to gain a better understanding of factors that lead to observed differences between epidemic potentials of the different variants.

Adaptation of HIV-1/HIV-2 Chimeras with Defects in Genome Packaging and Viral Replication

Frequent recombination is a hallmark of retrovirus replication. In rare cases, recombination occurs between distantly related retroviruses, generating novel viruses that may significantly impact viral evolution and public health. These recombinants may initially have substantial replication defects due to impaired interactions between proteins and/or nucleic acids from the two parental viruses. However, given the high mutation rates of retroviruses, these recombinants may be able to evolve improved compatibility of these viral elements. To test this hypothesis, we examined the adaptation of chimeras between two distantly related human pathogens: HIV-1 and HIV-2. We constructed HIV-1-based chimeras containing the HIV-2 nucleocapsid (NC) domain of Gag or the two zinc fingers of HIV-2 NC, which are critical for specific recognition of viral RNA. These chimeras exhibited significant defects in RNA genome packaging and replication kinetics in T cells. However, in some experiments, the chimeric viruses replicated with faster kinetics when repassaged, indicating that viral adaptation had occurred. Sequence analysis revealed the acquisition of a single amino acid substitution, S18L, in the first zinc finger of HIV-2 NC. This substitution, which represents a switch from a conserved HIV-2 residue to a conserved HIV-1 residue at this position, partially rescued RNA packaging and replication kinetics. Further analysis revealed that the combination of two substitutions in HIV-2 NC, W10F and S18L, almost completely restored RNA packaging and replication kinetics. Our study demonstrates that chimeras of distantly related retroviruses can adapt and significantly enhance their replication by acquiring a single substitution. IMPORTANCE Novel retroviruses can emerge from recombination between distantly related retroviruses. Most notably, HIV-1 originated from zoonotic transmission of a novel recombinant (SIV_cpz) into humans. Newly generated recombinants may initially have significant replication defects due to impaired interactions between viral proteins and/or nucleic acids, such as between cis- and trans-acting elements from the two parental viruses. However, provided that the recombinants retain some ability to replicate, they may be able to adapt and repair the defective interactions. Here, we used HIV-1 and HIV-2 Gag chimeras as a model system for studying the adaptation of recombinant viruses. We found that only two substitutions in the HIV-2 NC domain, W10F and S18L, were required to almost fully restore RNA genome packaging and replication kinetics. These results illustrate the extremely flexible nature of retroviruses and highlight the possible emergence of novel recombinants in the future that could pose a significant threat to public health.

Revealing the uncharacterised diversity of amphibian and reptile viruses

Amphibians and non-avian reptiles represent a significant proportion of terrestrial vertebrates, however knowledge of their viruses is not proportional to their abundance. Many amphibians and reptiles have strict habitual environments and localised populations and are vulnerable to viral outbreaks and potential elimination as a result. We sought to identify viruses that were hidden in amphibian and reptile metatranscriptomic data by screening 235 RNA-sequencing datasets from a 122 species covering 25 countries. We identified 26 novel viruses and eight previously characterised viruses from fifteen different viral families. Twenty-five viruses had RNA genomes with identity to Arteriviridae, Tobaniviridae, Hantaviridae, Rhabdoviridae, Astroviridae, Arenaviridae, Hepeviridae, Picornaviridae, Orthomyxoviridae, Reoviridae, Flaviviridae and Caliciviridae. In addition to RNA viruses, we also screened datasets for DNA viral transcripts, which are commonly excluded from transcriptomic analysis. We identified ten DNA viruses with identity to Papillomaviridae, Parvoviridae, Circoviridae and Adomaviridae. With the addition of these viruses, we expand the global amphibian and reptile virome and identify new potentially pathogenic viruses that could challenge populations. We speculate that amphibian viruses often have simpler genomes than those in amniotes, as in the case of the Secondpapillomavirinae and Orthomyxoviridae viruses identified in this study. In addition, we find evidence of inter-family recombination in RNA viruses, and we also identify new members of the recombinant Adomaviridae family. Overall, we provide insights into the uncharacterised diversity of amphibian and reptile viruses with the aim of improving population management, treatment and conservation into the future.

Molecular Characterization of Whole-Genome SARS-CoV-2 from the First Suspected Cases of the XE Variant in the Lazio Region, Italy

We report two cases of SARS-CoV-2 recombinant variant XE detected in nasopharyngeal swabs (NPS) of hospitalized patients with no evident epidemiological link in Lazio, Central Italy. Whole-Genome Sequencing (WGS) performed on an Ion Torrent GSS5 platform according to Italian flash surveys showed genomes corresponding to the PANGOLIN unclassified lineage and the Nextclade XE clade. Further analyses were then carried out to investigate more deeply the genetic characteristics of these XE-like sequences. When phylogenetic trees, by using IQ-TREE, were built splitting the genome into two regions according to the putative XE recombination site, the upstream and downstream regions were seen to be clustered near BA.1 and BA.2 sequences, respectively. However, our XE-like sequences clustered separately, with a significant bootstrap, from the classified European and Italian XE strains, although the recombination site between BA.1 and BA.2 was identified at the nucleotide site 11556 by RDP4 software, consistent with the putative XE breakpoint. These findings show the risk of the introduction of novel recombinant variants of SARS-CoV-2 and the existence of XE-like strains, phylogenetically separated, that could make their exact taxonomy difficult. It follows the need for continued SARS-CoV-2 surveillance by WGS.

Co-existence and co-infection of influenza A viruses and coronaviruses: Public health challenges

Since the 20^th century, humans have lived through five pandemics caused by influenza A viruses (IAVs) (H1N1/1918, H2N2/1957, H3N2/1968, and H1N1/2009) and the coronavirus (CoV) severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). IAVs and CoVs both have broad host ranges and share multiple hosts. Virus co-circulation and even co-infections facilitate genetic reassortment among IAVs and recombination among CoVs, further altering virus evolution dynamics and generating novel variants with increased cross-species transmission risk. Moreover, SARS-CoV-2 may maintain long-term circulation in humans as seasonal IAVs. Co-existence and co-infection of both viruses in humans could alter disease transmission patterns and aggravate disease burden. Herein, we demonstrate how virus-host ecology correlates with the co-existence and co-infection of IAVs and/or CoVs, further affecting virus evolution and disease dynamics and burden, calling for active virus surveillance and countermeasures for future public health challenges.

Culture and identification of a "Deltamicron" SARS-CoV-2 in a three cases cluster in southern France

Multiple SARS-CoV-2 variants have successively, or concomitantly spread worldwide since the summer of 2020. A few co-infections with different variants were reported and genetic recombinations, common among coronaviruses, were reported or suspected based on co-detection of signature mutations of different variants in a given genome. Here we report three infections in southern France with a Delta 21J_AY.4-Omicron 21K/BA.1 "Deltamicron" recombinant. The hybrid genome harbors signature mutations of the two lineages, supported by a mean sequencing depth of 1163-1421 reads and a mean nucleotide diversity of 0.1%-0.6%. It is composed of the near full-length spike gene (from codons 156-179) of an Omicron 21K/BA.1 variant in a Delta 21J/AY.4 lineage backbone. Importantly, we cultured an isolate of this recombinant and sequenced its genome. It was observed by scanning electron microscopy. As it is misidentified with current variant screening quantitative polymerase chain reaction (qPCR), we designed and implemented for routine diagnosis a specific duplex qPCR. Finally, structural analysis of the recombinant spike suggested its hybrid content could optimize viral binding to the host cell membrane. These findings prompt further studies of the virological, epidemiological, and clinical features of this recombinant.

Sequential Appearance and Isolation of a SARS-CoV-2 Recombinant between Two Major SARS-CoV-2 Variants in a Chronically Infected Immunocompromised Patient

Genetic recombination is a major evolutionary mechanism among RNA viruses, and it is common in coronaviruses, including those infecting humans. A few SARS-CoV-2 recombinants have been reported to date whose genome harbored combinations of mutations from different mutants or variants, but only a single patient's sample was analyzed, and the virus was not isolated. Here, we report the gradual emergence of a hybrid genome of B.1.160 and Alpha variants in a lymphoma patient chronically infected for 14 months, and we isolated the recombinant virus. The hybrid genome was obtained by next-generation sequencing, and the recombination sites were confirmed by PCR. This consisted of a parental B.1.160 backbone interspersed with two fragments, including the spike gene, from an Alpha variant. An analysis of seven sequential samples from the patient decoded the recombination steps, including the initial infection with a B.1.160 variant, then a concurrent infection with this variant and an Alpha variant, the generation of hybrid genomes, and eventually the emergence of a predominant recombinant virus isolated at the end of the patient's follow-up. This case exemplifies the recombination process of SARS-CoV-2 in real life, and it calls for intensifying the genomic surveillance in patients coinfected with different SARS-CoV-2 variants, and more generally with several RNA viruses, as this may lead to the appearance of new viruses.

Recombination in Coronaviruses, with a Focus on SARS-CoV-2

Recombination is a common evolutionary tool for RNA viruses, and coronaviruses are no exception. We review here the evidence for recombination in SARS-CoV-2 and reconcile nomenclature for recombinants, discuss their origin and fitness, and speculate how recombinants could make a difference in the future of the COVID-19 pandemics.

Crowning Touches in Positive-Strand RNA Virus Genome Replication Complex Structure and Function

Positive-strand RNA viruses, the largest genetic class of eukaryotic viruses, include coronaviruses and many other established and emerging pathogens. A major target for understanding and controlling these viruses is their genome replication, which occurs in virus-induced membrane vesicles that organize replication steps and protect double-stranded RNA intermediates from innate immune recognition. The structure of these complexes has been greatly illuminated by recent cryo-electron microscope tomography studies with several viruses. One key finding in diverse systems is the organization of crucial viral RNA replication factors in multimeric rings or crowns that among other functions serve as exit channels gating release of progeny genomes to the cytosol for translation and encapsidation. Emerging results suggest that these crowns serve additional important purposes in replication complex assembly, function, and interaction with downstream processes such as encapsidation. The findings provide insights into viral function and evolution and new bases for understanding, controlling, and engineering positive-strand RNA viruses. Expected final online publication date for the Annual Review of Virology, Volume 9 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Recombination in Positive-Strand RNA Viruses

RNA recombination is a major driver of genetic shifts tightly linked to the evolution of RNA viruses. Genomic recombination contributes substantially to the emergence of new viral lineages, expansion in host tropism, adaptations to new environments, and virulence and pathogenesis. Here, we review some of the recent progress that has advanced our understanding of recombination in positive-strand RNA viruses, including recombination triggers and the mechanisms behind them. The study of RNA recombination aids in predicting the probability and outcome of viral recombination events, and in the design of viruses with reduced recombination frequency as candidates for the development of live attenuated vaccines. Surveillance of viral recombination should remain a priority in the detection of emergent viral strains, a goal that can only be accomplished by expanding our understanding of how these events are triggered and regulated.

Generated Randomly and Selected Functionally? The Nature of Enterovirus Recombination

Genetic recombination in RNA viruses is an important evolutionary mechanism. It contributes to population diversity, host/tissue adaptation, and compromises vaccine efficacy. Both the molecular mechanism and initial products of recombination are relatively poorly understood. We used an established poliovirus-based in vitro recombination assay to investigate the roles of sequence identity and RNA structure, implicated or inferred from an analysis of circulating recombinant viruses, in the process. In addition, we used next-generation sequencing to investigate the early products of recombination after cellular coinfection with different poliovirus serotypes. In independent studies, we find no evidence for a role for RNA identity or structure in determining recombination junctions location. Instead, genome function and fitness are of greater importance in determining the identity of recombinant progeny. These studies provide further insights into this important evolutionary mechanism and emphasize the critical nature of the selection process on a mixed virus population.

Genomic Bootstrap Barcodes and Their Application to Study the Evolution of Sarbecoviruses

Recombination creates mosaic genomes containing regions with mixed ancestry, and the accumulation of such events over time can complicate greatly many aspects of evolutionary inference. Here, we developed a sliding window bootstrap (SWB) method to generate genomic bootstrap (GB) barcodes to highlight the regions supporting phylogenetic relationships. The method was applied to an alignment of 56 sarbecoviruses, including SARS-CoV and SARS-CoV-2, responsible for the SARS epidemic and COVID-19 pandemic, respectively. The SWB analyses were also used to construct a consensus tree showing the most reliable relationships and better interpret hidden phylogenetic signals. Our results revealed that most relationships were supported by just a few genomic regions and confirmed that three divergent lineages could be found in bats from Yunnan: SCoVrC, which groups SARS-CoV related coronaviruses from China; SCoV2rC, which includes SARS-CoV-2 related coronaviruses from Southeast Asia and Yunnan; and YunSar, which contains a few highly divergent viruses recently described in Yunnan. The GB barcodes showed evidence for ancient recombination between SCoV2rC and YunSar genomes, as well as more recent recombination events between SCoVrC and SCoV2rC genomes. The recombination and phylogeographic patterns suggest a strong host-dependent selection of the viral RNA-dependent RNA polymerase. In addition, SARS-CoV-2 appears as a mosaic genome composed of regions sharing recent ancestry with three bat SCoV2rCs from Yunnan (RmYN02, RpYN06, and RaTG13) or related to more ancient ancestors in bats from Yunnan and Southeast Asia. Finally, our results suggest that viral circular RNAs may be key molecules for the mechanism of recombination.

WGS- versus ORF5-Based Typing of PRRSV: A Belgian Case Study

Porcine reproductive and respiratory syndrome virus (PRRSV) is the causative agent of one of the most widespread and economically devastating diseases in the swine industry. Typing circulating PRRSV strains by means of sequencing is crucial for developing adequate control strategies. Most genetic studies only target the highly variable open reading frame (ORF) 5, for which an extensive database is available. In this study, we performed whole-genome sequencing (WGS) on a collection of 124 PRRSV-1 positive serum samples that were collected over a 5-year period (2015–2019) in Belgium. Our results show that (nearly) complete PRRSV genomes can be obtained directly from serum samples with a high success rate. Analysis of the coding regions confirmed the exceptionally high genetic diversity, even among Belgian PRRSV-1 strains. To gain more insight into the added value of WGS, we performed phylogenetic cluster analyses on separate ORF datasets as well as on a single, concatenated dataset (CDS) containing all ORFs. A comparison between the CDS and ORF clustering schemes revealed numerous discrepancies. To explain these differences, we performed a large-scale recombination analysis, which allowed us to identify a large number of potential recombination events that were scattered across the genome. As PRRSV does not contain typical recombination hot-spots, typing PRRSV strains based on a single ORF is not recommended. Although the typing accuracy can be improved by including multiple regions, our results show that the full genetic diversity among PRRSV strains can only be captured by analysing (nearly) complete genomes. Finally, we also identified several vaccine-derived recombinant strains, which once more raises the question of the safety of these vaccines.

First come, first served: superinfection exclusion in Deformed wing virus is dependent upon sequence identity and not the order of virus acquisition

Deformed wing virus (DWV) is the most important globally distributed pathogen of honey bees and, when vectored by the ectoparasite Varroa destructor, is associated with high levels of colony losses. Divergent DWV types may differ in their pathogenicity and are reported to exhibit superinfection exclusion upon sequential infections, an inevitability in a Varroa-infested colony. We used a reverse genetic approach to investigate competition and interactions between genetically distinct or related virus strains, analysing viral load over time, tissue distribution with reporter gene-expressing viruses and recombination between virus variants. Transient competition occurred irrespective of the order of virus acquisition, indicating no directionality or dominance. Over longer periods, the ability to compete with a pre-existing infection correlated with the genetic divergence of the inoculae. Genetic recombination was observed throughout the DWV genome with recombinants accounting for ~2% of the population as determined by deep sequencing. We propose that superinfection exclusion, if it occurs at all, is a consequence of a cross-reactive RNAi response to the viruses involved, explaining the lack of dominance of one virus type over another. A better understanding of the consequences of dual- and superinfection will inform development of cross-protective honey bee vaccines and landscape-scale DWV transmission and evolution.

Application Route and Immune Status of the Host Determine Safety and Oncolytic Activity of Oncolytic Coxsackievirus B3 Variant PD-H

The coxsackievirus B3 strain PD-0 has been proposed as a new oncolytic virus for the treatment of colorectal carcinoma. Here, we generated a cDNA clone of PD-0 and analyzed the virus PD-H, newly generated from this cDNA, in xenografted and syngenic models of colorectal cancer. Replication and cytotoxic assays revealed that PD-H replicated and lysed colorectal carcinoma cell lines in vitro as well as PD-0. Intratumoral injection of PD-H into subcutaneous DLD-1 tumors in nude mice resulted in strong inhibition of tumor growth and significantly prolonged the survival of the animals, but virus-induced systemic infection was observed in one of the six animals. In a syngenic mouse model of subcutaneously growing Colon-26 tumors, intratumoral administration of PD-H led to a significant reduction of tumor growth, the prolongation of animal survival, the prevention of tumor-induced cachexia, and the elevation of CD3+ and dendritic cells in the tumor microenvironment. No virus-induced side effects were observed. After intraperitoneal application, PD-H induced weak pancreatitis and myocarditis in immunocompetent mice. By equipping the virus with target sites of miR-375, which is specifically expressed in the pancreas, organ infections were prevented. Moreover, employment of this virus in a syngenic mouse model of CT-26 peritoneal carcinomatosis resulted in a significant reduction in tumor growth and an increase in animal survival. The results demonstrate that the immune status of the host, the route of virus application, and the engineering of the virus with target sites of suitable microRNAs are crucial for the use of PD-H as an oncolytic virus.

Imprecise recombinant viruses evolve via a fitness-driven, iterative process of polymerase template-switching events

Recombination is a common feature of many positive-strand RNA viruses, playing an important role in virus evolution. However, to date, there is limited understanding of the mechanisms behind the process. Utilising in vitro assays, we have previously shown that the template-switching event of recombination is a random and ubiquitous process that often leads to recombinant viruses with imprecise genomes containing sequence duplications. Subsequently, a process termed resolution, that has yet to be mechanistically studied, removes these duplicated sequences resulting in a virus population of wild type length genomes. Using defined imprecise recombinant viruses together with Oxford Nanopore and Illumina high throughput next generation sequencing technologies we have investigated the process of resolution. We show that genome resolution involves subsequent rounds of template-switching recombination with viral fitness resulting in the survival of a small subset of recombinant genomes. This alters our previously held understanding that recombination and resolution are independent steps of the process, and instead demonstrates that viruses undergo frequent and continuous recombination events over a prolonged period until the fittest viruses, predominantly those with wild type length genomes, dominate the population.

Viruses with positive-sense RNA genomes, such as poliovirus, have several mechanisms by which they evolve. One of these is the process of recombination involving the large-scale exchange of genetic information. Recombination occurs during replication when the viral polymerase, bound to the nascent RNA chain, switches from copying one genome to another. However, the polymerase does not always accurately switch between the two, resulting in sequence duplications or deletions, and genomes that are referred to as imprecise. Over multiple rounds of replication sequence duplications are lost and genomes are resolved to wild type length, but it is unclear how this occurs. Here we used synthetic polioviruses containing defined sequence duplications to determine that the genome population undergoes repeated rounds of recombination until sequence duplications are lost and viruses with precise, wild type length genomes are selected for. This selection is based on the overall fitness of the virus population, with less fit imprecise viruses evolving more quickly. Our study suggests that recombination is a continual process where virus fitness drives the selection of a small subset of recombinant variants. These data are important for understanding how novel viruses evolve via recombination and how this process can be blocked to prevent novel and dangerous pathogens from arising.

Origin, phylogeny, variability and epitope conservation of SARS-CoV-2 worldwide

The coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses innumerous challenges, like understanding what triggered the emergence of this new human virus, how this RNA virus is evolving or how the variability of viral genome may impact the primary structure of proteins that are targets for vaccine. We analyzed 19471 SARS-CoV-2 genomes available at the GISAID database from all over the world and 3335 genomes of other Coronoviridae family members available at GenBank, collecting SARS-CoV-2 high-quality genomes and distinct Coronoviridae family genomes. Additionally, we analyzed 199,984 spike glycoprotein sequences. Here, we identify a SARS-CoV-2 emerging cluster containing 13 closely related genomes isolated from bat and pangolin that showed evidence of recombination, which may have contributed to the emergence of SARS-CoV-2. The analyzed SARS-CoV-2 genomes presented 9632 single nucleotide variants (SNVs) corresponding to a variant density of 0.3 over the genome, and a clear geographic distribution. SNVs are unevenly distributed throughout the genome and hotspots for mutations were found for the spike gene and ORF 1ab. We describe a set of predicted spike protein epitopes whose variability is negligible. Additionally, all predicted epitopes for the structural E, M and N proteins are highly conserved. The amino acid changes present in the spike glycoprotein of variables of concern (VOCs) comprise between 3.4% and 20.7% of the predicted epitopes of this protein. These results favors the continuous efficacy of the available vaccines targeting the spike protein, and other structural proteins. Multiple epitopes vaccines should sustain vaccine efficacy since at least some of the epitopes present in variability regions of VOCs are conserved and thus recognizable by antibodies.

Genetic Characteristics of Human Parainfluenza Virus Types 1-4 From Patients With Clinical Respiratory Tract Infection in China

Human parainfluenza viruses (HPIV1–4) cause acute respiratory tract infections, thereby impacting human health worldwide. However, there are no current effective antivirals or licensed vaccines for infection prevention. Moreover, sequence information for human parainfluenza viruses (HPIVs) circulating in China is inadequate. Therefore, to shed light on viral genetic diversity and evolution, we collected samples from patients infected with HPIV1–4 in China from 2012 to 2018 to sequence the viruses. We obtained 24 consensus sequences, comprising 1 for HPIV1, 2 for HPIV2, 19 for HPIV3, and 2 for HPIV4A. Phylogenetic analyses classified the 1 HPIV1 into clade 2, and the 2 HPIV4 sequences into cluster 4A. Based on the hemagglutinin-neuraminidase (HN) gene, a new sub-cluster was identified in one of the HPIV2, namely G1c, and the 19 HPIV3 sequences were classified into the genetic lineages of C3f and C3a. The results indicated that HPIV1–4 were co-circulated in China. Further, the lineages of sub-cluster C3 of HPIV3 were co-circulated in China. A recombination analysis indicated that a putative recombination event may have occurred in the HN gene of HPIV3. In the obtained sequences of HPIV3, we found that two amino acid substitution sites (R73K in the F protein of PUMCH14028/2014 and A281V in the HN protein of PUMCH13961/2014) and a negative selection site (amino acid position 398 in the F protein) corresponded to the previously reported neutralization-related sites. Moreover, amino acid substitution site (K108E) corresponded to the negative selection site (amino acid position 108) in the 10 F proteins of HPIV3. However, no amino acid substitution site corresponded to the glycosylation site in the obtained HPIV3 sequences. These results might help in studying virus evolution, developing vaccines, and monitoring HPIV-related respiratory diseases.

High-Throughput Fluorescent Assay for Inhibitor Screening of Proteases from RNA Viruses

Spanish flu, polio epidemics, and the ongoing COVID-19 pandemic are the most profound examples of severe widespread diseases caused by RNA viruses. The coronavirus pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) demands affordable and reliable assays for testing antivirals. To test inhibitors of viral proteases, we have developed an inexpensive high-throughput assay based on fluorescent energy transfer (FRET). We assayed an array of inhibitors for papain-like protease from SARS-CoV-2 and validated it on protease from the tick-borne encephalitis virus to emphasize its versatility. The reaction progress is monitored as loss of FRET signal of the substrate. This robust and reproducible assay can be used for testing the inhibitors in 96- or 384-well plates.

Full-Genome Sequences of Alphacoronaviruses and Astroviruses from Myotis and Pipistrelle Bats in Denmark

Bat species worldwide are receiving increased attention for the discovery of emerging viruses, cross-species transmission, and zoonoses, as well as for characterizing virus infections specific to bats. In a previous study, we investigated the presence of coronaviruses in faecal samples from bats at different locations in Denmark, and made phylogenies based on short, partial ORF1b sequences. In this study, selected samples containing bat coronaviruses from three different bat species were analysed, using a non-targeted approach of next-generation sequencing. From the resulting metagenomics data, we assembled full-genome sequences of seven distinct alphacoronaviruses, three astroviruses, and a polyomavirus, as well as partial genome sequences of rotavirus H and caliciviruses, from the different bat species. Comparisons to published sequences indicate that the bat alphacoronaviruses belong to three different subgenera-i.e., Pedacovirus, Nyctacovirus, and Myotacovirus-that the astroviruses may be new species in the genus Mamastrovirus, and that the polyomavirus could also be a new species, but unassigned to a genus. Furthermore, several viruses of invertebrates-including two Rhopalosiphum padi (aphid) viruses and a Kadipiro virus-present in the faecal material were assembled. Interestingly, this is the first detection in Europe of a Kadipiro virus.

Rapid protein sequence evolution via compensatory frameshift is widespread in RNA virus genomes

Background

RNA viruses possess remarkable evolutionary versatility driven by the high mutability of their genomes. Frameshifting nucleotide insertions or deletions (indels), which cause the premature termination of proteins, are frequently observed in the coding sequences of various viral genomes. When a secondary indel occurs near the primary indel site, the open reading frame can be restored to produce functional proteins, a phenomenon known as the compensatory frameshift.

Results

In this study, we systematically analyzed publicly available viral genome sequences and identified compensatory frameshift events in hundreds of viral protein-coding sequences. Compensatory frameshift events resulted in large-scale amino acid differences between the compensatory frameshift form and the wild type even though their nucleotide sequences were almost identical. Phylogenetic analyses revealed that the evolutionary distance between proteins with and without a compensatory frameshift were significantly overestimated because amino acid mismatches caused by compensatory frameshifts were counted as substitutions. Further, this could cause compensatory frameshift forms to branch in different locations in the protein and nucleotide trees, which may obscure the correct interpretation of phylogenetic relationships between variant viruses.

Conclusions

Our results imply that the compensatory frameshift is one of the mechanisms driving the rapid protein evolution of RNA viruses and potentially assisting their host-range expansion and adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04182-9.

Characterizing Transcriptional Regulatory Sequences in Coronaviruses and Their Role in Recombination

Novel coronaviruses, including SARS-CoV-2, SARS, and MERS, often originate from recombination events. The mechanism of recombination in RNA viruses is template switching. Coronavirus transcription also involves template switching at specific regions, called transcriptional regulatory sequences (TRS). It is hypothesized but not yet verified that TRS sites are prone to recombination events. Here, we developed a tool called SuPER to systematically identify TRS in coronavirus genomes and then investigated whether recombination is more common at TRS. We ran SuPER on 506 coronavirus genomes and identified 465 TRS-L and 3,509 TRS-B. We found that the TRS-L core sequence (CS) and the secondary structure of the leader sequence are generally conserved within coronavirus genera but different between genera. By examining the location of recombination breakpoints with respect to TRS-B CS, we observed that recombination hotspots are more frequently colocated with TRS-B sites than expected.

Structure and Expression of Large (+)RNA Genomes of Viruses of Higher Eukaryotes

Viral positive-sense RNA genomes evolve rapidly due to the high mutation rates during replication and RNA recombination, which allowing the viruses to acquire and modify genes for their adaptation. The size of RNA genome is limited by several factors, including low fidelity of RNA polymerases and packaging constraints. However, the 12-kb size limit is exceeded in the two groups of eukaryotic (+)RNA viruses – animal nidoviruses and plant closteroviruses. These virus groups have several traits in common. Their genomes contain 5′-proximal genes that are expressed via ribosomal frameshifting and encode one or two papain-like protease domains, membrane-binding domain(s), methyltransferase, RNA helicase, and RNA polymerase. In addition, some nidoviruses (i.e., coronaviruses) contain replication-associated domains, such as proofreading exonuclease, putative primase, nucleotidyltransferase, and endonuclease. In both nidoviruses and closteroviruses, the 3′-terminal part of the genome contains genes for structural and accessory proteins expressed via a nested set of coterminal subgenomic RNAs. Coronaviruses and closteroviruses have evolved to form flexuous helically symmetrical nucleocapsids as a mean to resolve packaging constraints. Since phylogenetic reconstructions of the RNA polymerase domains indicate only a marginal relationship between the nidoviruses and closteroviruses, their similar properties likely have evolved convergently, along with the increase in the genome size.

Viral mouse models used to study multiple sclerosis: past and present

Multiple sclerosis (MS) is a common inflammatory demyelinating disease of the central nervous system. Although the etiology of MS is unknown, genetics and environmental factors, such as infections, play a role. Viral infections of mice have been used as model systems to study this demyelinating disease of humans. Three viruses that have long been studied in this capacity are Theiler’s murine encephalomyelitis virus, mouse hepatitis virus, and Semliki Forest virus. This review describes the viruses themselves, the infection process, the disease caused by infection and its accompanying pathology, and the model systems and their usefulness in studying MS.

Development and application of reverse genetic technology for the influenza virus

Influenza virus is a common virus in people's daily lives, and it has certain infectivity in humans and animals. Influenza viruses have the characteristics of a high mutation rate and wide distribution. Reverse genetic technology is primarily used to modify viruses at the DNA level through targeted modification of the virus cDNA. Genetically modified influenza viruses have a unique advantage when researching the transmission and pathogenicity of influenza. With the continuous development of oncolytic viruses in recent years, studies have found that influenza viruses also have certain oncolytic activity. Influenza viruses can specifically recognize tumor cells; activate cytotoxic T cells, NK cells, dendritic cells, etc.; and stimulate the body to produce an immune response, thereby killing tumor cells. This article will review the development and application of influenza virus reverse genetic technology.

Molecular epidemiology of coxsackievirus A16 circulating in children in Beijing, China from 2010 to 2019

BACKGROUND

Coxsackievirus A16 (CVA16) is one of the major etiological agents of hand, foot and mouth disease (HFMD). This study aimed to investigate the molecular epidemiology and evolutionary characteristics of CVA16.

METHODS

Throat swabs were collected from children with HFMD and suspected HFMD during 2010-2019. Enteroviruses (EVs) were detected and typed by real-time reverse transcription-polymerase chain reaction (RT-PCR) and RT-PCR. The genotype, evolutionary rate, the most recent common ancestor, population dynamics and selection pressure of CVA16 were analyzed based on viral protein gene (VP1) by bioinformatics software.

RESULTS

A total of 4709 throat swabs were screened. EVs were detected in 3180 samples and 814 were CVA16 positive. More than 81% of CVA16-positive children were under 5 years old. The prevalence of CVA16 showed obvious periodic fluctuations with a high level during 2010-2012 followed by an apparent decline during 2013-2017. However, the activities of CVA16 increased gradually during 2018-2019. All the Beijing CVA16 strains belonged to sub-genotype B1, and B1b was the dominant strain. One B1c strain was detected in Beijing for the first time in 2016. The estimated mean evolutionary rate of VP1 gene was 4.49 × 10^-3 substitution/site/year. Methionine gradually fixed at site-23 of VP1 since 2012. Two sites were detected under episodic positive selection, one of which (site-223) located in neutralizing linear epitope PEP71.

CONCLUSIONS

The dominant strains of CVA16 belonged to clade B1b and evolved in a fast evolutionary rate during 2010-2019 in Beijing. To provide more favorable data for HFMD prevention and control, it is necessary to keep attention on molecular epidemiological and evolutionary characteristics of CVA16.

Two distant helicases in one mycovirus: evidence of horizontal gene transfer between mycoviruses, coronaviruses and other nidoviruses

Nidovirales, which accommodates viruses with the largest RNA genomes, includes the notorious coronaviruses; however, the evolutionary route for nidoviruses is not well understood. We have characterized a positive-sense (+) single-stranded (ss) RNA mycovirus, Rhizoctonia solani hypovirus 2 (RsHV2), from the phytopathogenic fungus Rhizoctonia solani. RsHV2 has the largest RNA genome size of 22,219 nucleotides, excluding the poly(A) tail, in all known mycoviruses, and contains two open reading frames (ORF1 and ORF2). ORF1 encodes a protein of 2,009 amino acid (aa) that includes a conserved helicase domain belonging to helicase superfamily I (SFI). In contrast, ORF2 encodes a polyprotein of 4459 aa containing the hallmark genes of hypoviruses. The latter includes a helicase belonging to SFII. Following phylogenetic analysis, the ORF1-encoded helicase (Hel1) unexpectedly clustered in an independent evolutionary branch together with nidovirus helicases, including coronaviruses, and bacteria helicases. Thus, Hel1 presence indicates the occurrence of horizontal gene transfer between viruses and bacteria. These findings also suggest that RsHV2 is most likely a recombinant arising between hypoviruses and nidoviruses.

Deformed wing virus: using reverse genetics to tackle unanswered questions about the most important viral pathogen of honey bees

Deformed wing virus (DWV) is the most important viral pathogen of honey bees. It usually causes asymptomatic infections but, when vectored by the ectoparasitic mite Varroa destructor, it is responsible for the majority of overwintering colony losses globally. Although DWV was discovered four decades ago, research has been hampered by the absence of an in vitro cell culture system or the ability to culture pure stocks of the virus. The recent developments of reverse genetic systems for DWV go some way to addressing these limitations. They will allow the investigation of specific questions about strain variation, host tropism and pathogenesis to be answered, and are already being exploited to study tissue tropism and replication in Varroa and non-Apis pollinators. Three areas neatly illustrate the advances possible with reverse genetic approaches: (i) strain variation and recombination, in which reverse genetics has highlighted similarities rather than differences between virus strains; (ii) analysis of replication kinetics in both honey bees and Varroa, in studies that likely explain the near clonality of virus populations often reported; and (iii) pathogen spillover to non-Apis pollinators, using genetically tagged viruses to accurately monitor replication and infection.

Molecular characterization of porcine epidemic diarrhea virus associated with outbreaks in southwest China during 2014-2018

Porcine epidemic diarrhea virus (PEDV), which re-emerged in China since 2010, has swept across the whole country leading to tremendous economic losses. In this study, a total of 645 diarrhea samples collected from 156 pig farms in Sichuan and Guizhou province during 2014-2018 were tested for PEDV. We found that samples from 47.66% (84/156) of the farms were positive for PEDV with an overall detection rate of 35.81% (231/645). Fifty-two strains were selected for full-length S gene analyses, and these strains were classified into three subgroups, an S-INDEL subgroup (G1c), and two non-S-INDEL subgroups (G2b, AJ1102-like and G2c), accounting for 15.38% (8/52), 23.08% (12/52) and 59.62% (31/52) of the total analysed strains, respectively. We found these three subgroups of PEDV coexisted in Sichuan province, and the S-INDEL strain was detected in Guizhou. Further antigenic variation analysis of the neutralizing epitopes (S10, COE, SS2, SS6 and 2C10) on the spike protein revealed that the S-INDEL and non-S-INDEL strains shared similar variation features in COE and SS6, but exhibited distinct variation patterns in the S10 domain. Unique variation patterns on N-glycosylation sites in the S protein were also observed for the S-INDEL and non-S-INDEL strains. Moreover, nine strains (three from each subgroup) were subjected to full-genome characterization. Complete genome phylogeny showed an inconsistent tree topology for genotyping, with two G2c strains grouped into the GII-b (AH2012-like) genogroup and the remaining seven strains including three S-INDEL strains grouped into the GII-c genogroup. Further recombination analyses indicated that six of the GII-c strains probably originated from intra-genogroup recombinations. Notably, three newly emerged S-INDEL strains with novel recombination patterns were first identified. Together, our data revealed a new status of PEDV in southwest China, which can increase understanding of the prevalence, genetic characteristics and evolutionary profiles of circulating PEDV strains in China.

Genomic recombination events may reveal the evolution of coronavirus and the origin of SARS-CoV-2

To trace the evolution of coronaviruses and reveal the possible origin of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes the coronavirus disease 2019 (COVID-19), we collected and thoroughly analyzed 29,452 publicly available coronavirus genomes, including 26,312 genomes of SARS-CoV-2 strains. We observed coronavirus recombination events among different hosts including 3 independent recombination events with statistical significance between some isolates from humans, bats and pangolins. Consistent with previous records, we also detected putative recombination between strains similar or related to Bat-CoV-RaTG13 and Pangolin-CoV-2019. The putative recombination region is located inside the receptor-binding domain (RBD) of the spike glycoprotein (S protein), which may represent the origin of SARS-CoV-2. Population genetic analyses provide estimates suggesting that the putative introduced DNA within the RBD is undergoing directional evolution. This may result in the adaptation of the virus to hosts. Unsurprisingly, we found that the putative recombination region in S protein was highly diverse among strains from bats. Bats harbor numerous coronavirus subclades that frequently participate in recombination events with human coronavirus. Therefore, bats may provide a pool of genetic diversity for the origin of SARS-CoV-2.

Recombination and Positive Selection Differentially Shaped the Diversity of Betacoronavirus Subgenera

The Betacoronavirus genus of mammal-infecting viruses includes three subgenera (Sarbecovirus, Embecovirus, and Merbecovirus), in which most known human coronaviruses, including SARS-CoV-2, cluster. Coronaviruses are prone to host shifts, with recombination and positive selection possibly contributing to their high zoonotic potential. We analyzed the role of these two forces in the evolution of viruses belonging to the Betacoronavirus genus. The results showed that recombination has been pervasive during sarbecovirus evolution, and it is more widespread in this subgenus compared to the other two. In both sarbecoviruses and merbecoviruses, recombination hotspots are clearly observed. Conversely, positive selection was a less prominent force in sarbecoviruses compared to embecoviruses and merbecoviruses and targeted distinct genomic regions in the three subgenera, with S being the major target in sarbecoviruses alone. Overall, the results herein indicate that Betacoronavirus subgenera evolved along different trajectories, which might recapitulate their host preferences or reflect the origins of the presently available coronavirus sequences.

Polycistronic Genome Segment Evolution and Gain and Loss of FAST Protein Function during Fusogenic Orthoreovirus Speciation

The Reoviridae family is the only non-enveloped virus family with members that use syncytium formation to promote cell–cell virus transmission. Syncytiogenesis is mediated by a fusion-associated small transmembrane (FAST) protein, a novel family of viral membrane fusion proteins. Previous evidence suggested the fusogenic reoviruses arose from an ancestral non-fusogenic virus, with the preponderance of fusogenic species suggesting positive evolutionary pressure to acquire and maintain the fusion phenotype. New phylogenetic analyses that included the atypical waterfowl subgroup of avian reoviruses and recently identified new orthoreovirus species indicate a more complex relationship between reovirus speciation and fusogenic capacity, with numerous predicted internal indels and 5’-terminal extensions driving the evolution of the orthoreovirus’ polycistronic genome segments and their encoded FAST and fiber proteins. These inferred recombination events generated bi- and tricistronic genome segments with diverse gene constellations, they occurred pre- and post-orthoreovirus speciation, and they directly contributed to the evolution of the four extant orthoreovirus FAST proteins by driving both the gain and loss of fusion capability. We further show that two distinct post-speciation genetic events led to the loss of fusion in the waterfowl isolates of avian reovirus, a recombination event that replaced the p10 FAST protein with a heterologous, non-fusogenic protein and point substitutions in a conserved motif that destroyed the p10 assembly into multimeric fusion platforms.

Characterizing transcriptional regulatory sequences in coronaviruses and their role in recombination

Novel coronaviruses, including SARS-CoV-2, SARS, and MERS, often originate from recombination events. The mechanism of recombination in RNA viruses is template switching. Coronavirus transcription also involves template switching at specific regions, called transcriptional regulatory sequences (TRS). It is hypothesized but not yet verified that TRS sites are prone to recombination events. Here, we developed a tool called SuPER to systematically identify TRS in coronavirus genomes and then investigated whether recombination is more common at TRS. We ran SuPER on 506 coronavirus genomes and identified 465 TRS-L and 3509 TRS-B. We found that the TRS-L core sequence (CS) and the secondary structure of the leader sequence are generally conserved within coronavirus genera but different between genera. By examining the location of recombination breakpoints with respect to TRS-B CS, we observed that recombination hotspots are more frequently co-located with TRS-B sites than expected.

Coinfections and their molecular consequences in the porcine respiratory tract

Understudied, coinfections are more frequent in pig farms than single infections. In pigs, the term “Porcine Respiratory Disease Complex” (PRDC) is often used to describe coinfections involving viruses such as swine Influenza A Virus (swIAV), Porcine Reproductive and Respiratory Syndrome Virus (PRRSV), and Porcine CircoVirus type 2 (PCV2) as well as bacteria like Actinobacillus pleuropneumoniae, Mycoplasma hyopneumoniae and Bordetella bronchiseptica. The clinical outcome of the various coinfection or superinfection situations is usually assessed in the studies while in most of cases there is no clear elucidation of the fine mechanisms shaping the complex interactions occurring between microorganisms. In this comprehensive review, we aimed at identifying the studies dealing with coinfections or superinfections in the pig respiratory tract and at presenting the interactions between pathogens and, when possible, the mechanisms controlling them. Coinfections and superinfections involving viruses and bacteria were considered while research articles including protozoan and fungi were excluded. We discuss the main limitations complicating the interpretation of coinfection/superinfection studies, and the high potential perspectives in this fascinating research field, which is expecting to gain more and more interest in the next years for the obvious benefit of animal health.

Dissemination of Internal Ribosomal Entry Sites (IRES) Between Viruses by Horizontal Gene Transfer

Members of Picornaviridae and of the Hepacivirus, Pegivirus and Pestivirus genera of Flaviviridae all contain an internal ribosomal entry site (IRES) in the 5'-untranslated region (5'UTR) of their genomes. Each class of IRES has a conserved structure and promotes 5'-end-independent initiation of translation by a different mechanism. Picornavirus 5'UTRs, including the IRES, evolve independently of other parts of the genome and can move between genomes, most commonly by intratypic recombination. We review accumulating evidence that IRESs are genetic entities that can also move between members of different genera and even between families. Type IV IRESs, first identified in the Hepacivirus genus, have subsequently been identified in over 25 genera of Picornaviridae, juxtaposed against diverse coding sequences. In several genera, members have either type IV IRES or an IRES of type I, II or III. Similarly, in the genus Pegivirus, members contain either a type IV IRES or an unrelated type; both classes of IRES also occur in members of the genus Hepacivirus. IRESs utilize different mechanisms, have different factor requirements and contain determinants of viral growth, pathogenesis and cell type specificity. Their dissemination between viruses by horizontal gene transfer has unexpectedly emerged as an important facet of viral evolution.

Clinical sequelae of the novel coronavirus: does COVID-19 infection predispose patients to cancer?

As cancer patients are clinically known to be predisposed to COVID-19 infection, a corollary question of whether COVID-19 infection predisposes to cancer is explored. This article seeks to establish an association between novel coronavirus sequelae and cancer. A literature review on COVID-19 mechanisms of action, molecular responses it elicits upon infection and tumorigenesis pathways is conducted to establish this association. Major signaling pathways implicated in aberrant cellular growth are activated, the ensuing cytokine storm weakens the immune system response to tumors, and patients may develop cancer as a result of superimposed mutagenic and/or carcinogenic events. Future work needs to be performed to support this hypothesis, both in in vitro models and preclinical studies. COVID-19 patients may need to be monitored post-infection for developing cancer.

Genomic Analyses of Human Sapoviruses Detected over a 40-Year Period Reveal Disparate Patterns of Evolution among Genotypes and Genome Regions

Human sapovirus is a causative agent of acute gastroenteritis in all age groups. The use of full-length viral genomes has proven beneficial to investigate evolutionary dynamics and transmission chains. In this study, we developed a full-length genome sequencing platform for human sapovirus and sequenced the oldest available strains (collected in the 1970s) to analyse diversification of sapoviruses. Sequence analyses from five major genotypes (GI.1, GI.2, GII.1, GII.3, and GIV.1) showed limited intra-genotypic diversification for over 20–40 years. The accumulation of amino acid mutations in VP1 was detected for GI.2 and GIV.1 viruses, while having a similar rate of nucleotide evolution to the other genotypes. Differences in the phylogenetic clustering were detected between RdRp and VP1 sequences of our archival strains as well as other reported putative recombinants. However, the lack of the parental strains and differences in diversification among genomic regions suggest that discrepancies in the phylogenetic clustering of sapoviruses could be explained, not only by recombination, but also by disparate nucleotide substitution patterns between RdRp and VP1 sequences. Together, this study shows that, contrary to noroviruses, sapoviruses present limited diversification by means of intra-genotype variation and recombination.

Coding potential and sequence conservation of SARS-CoV-2 and related animal viruses

In December 2019, a novel human-infecting coronavirus (SARS-CoV-2) was recognized in China. In a few months, SARS-CoV-2 has caused thousands of disease cases and deaths in several countries. Phylogenetic analyses indicated that SARS-CoV-2 clusters with SARS-CoV in the Sarbecovirus subgenus and viruses related to SARS-CoV-2 were identified from bats and pangolins. Coronaviruses have long and complex genomes with high plasticity in terms of gene content. To date, the coding potential of SARS-CoV-2 remains partially unknown. We thus used available sequences of bat and pangolin viruses to determine the selective events that shaped the genome structure of SARS-CoV-2 and to assess its coding potential. By searching for signals of significantly reduced variability at synonymous sites (dS), we identified six genomic regions, one of these corresponding to the programmed −1 ribosomal frameshift. The most prominent signal of dS reduction was observed within the E gene. A genome-wide analysis of conserved RNA structures indicated that this region harbors a putative functional RNA element that is shared with the SARS-CoV lineage. Additional signals of reduced dS indicated the presence of internal ORFs. Whereas the presence ORF9a (internal to N) was previously proposed by homology with a well characterized protein of SARS-CoV, ORF3h (for hypothetical, within ORF3a) was not previously described. The predicted product of ORF3h has 90% identity with the corresponding predicted product of SARS-CoV and displays features suggestive of a viroporin. Finally, analysis of the putative ORF10 revealed high dN/dS (3.82) in SARS-CoV-2 and related coronaviruses. In the SARS-CoV lineage, the ORF is predicted to encode a truncated protein and is neutrally evolving. These data suggest that ORF10 encodes a functional protein in SARS-CoV-2 and that positive selection is driving its evolution. Experimental analyses will be necessary to validate and characterize the coding and non-coding functional elements we identified.

•

We analyzed the coding region of SARS-CoV-2 and related bat/pangolin viruses.

•

We identified six regions of significantly low variability at sysnonymous sites.

•

One of these corresponds to a conserved RNA structure shared with the SARS-CoV lineage.

•

The dS reduction within ORF3a corresponds to a potential ORF encoding a viroporin.

•

In SARS-CoV-2 and related viruses, the putative 3′ terminal ORF10 has high dN/dS.