Deep mining of the Sequence Read Archive reveals major genetic innovations in coronaviruses and other nidoviruses of aquatic vertebrates.

Virus discovery by genomics and metagenomics empowered studies of viromes, facilitated characterization of pathogen epidemiology, and redefined our understanding of the natural genetic diversity of viruses with profound functional and structural implications. Here we employed a data-driven virus discovery approach that directly queries unprocessed sequencing data in a highly parallelized way and involves a targeted viral genome assembly strategy in a wide range of sequence similarity. By screening more than 269,000 datasets of numerous authors from the Sequence Read Archive and using two metrics that quantitatively assess assembly quality, we discovered 40 nidoviruses from six virus families whose members infect vertebrate hosts. They form 13 and 32 putative viral subfamilies and genera, respectively, and include 11 coronaviruses with bisegmented genomes from fishes and amphibians, a giant 36.1 kilobase coronavirus genome with a duplicated spike glycoprotein (S) gene, 11 tobaniviruses and 17 additional corona-, arteri-, cremega-, nanhypo- and nangoshaviruses. Genome segmentation emerged in a single evolutionary event in the monophyletic lineage encompassing the subfamily Pitovirinae. We recovered the bisegmented genome sequences of two coronaviruses from RNA samples of 69 infected fishes and validated the presence of poly(A) tails at both segments using 3'RACE PCR and subsequent Sanger sequencing. We report a genetic linkage between accessory and structural proteins whose phylogenetic relationships and evolutionary distances are incongruent with the phylogeny of replicase proteins. We rationalize these observations in a model of inter-family S recombination involving at least five ancestral corona- and tobaniviruses of aquatic hosts. In support of this model, we describe an individual fish co-infected with members from the families Coronaviridae and Tobaniviridae. Our results expand the scale of the known extraordinary evolutionary plasticity in nidoviral genome architecture and call for revisiting fundamentals of genome expression, virus particle biology, host range and ecology of vertebrate nidoviruses.

Delving into the Aftermath of a Disease-Associated Near-Extinction Event: A Five-Year Study of a Serpentovirus (Nidovirus) in a Critically Endangered Turtle Population.

Bellinger River virus (BRV) is a serpentovirus (nidovirus) that was likely responsible for the catastrophic mortality of the Australian freshwater turtle Myuchelys georgesi in February 2015. From November 2015 to November 2020, swabs were collected from turtles during repeated river surveys to estimate the prevalence of BRV RNA, identify risk factors associated with BRV infection, and refine sample collection. BRV RNA prevalence at first capture was significantly higher in M. georgesi (10.8%) than in a coexisting turtle, Emydura macquarii (1.0%). For M. georgesi, various risk factors were identified depending on the analysis method, but a positive BRV result was consistently associated with a larger body size. All turtles were asymptomatic when sampled and conjunctival swabs were inferred to be optimal for ongoing monitoring. Although the absence of disease and recent BRV detections suggests a reduced ongoing threat, the potential for the virus to persist in an endemic focus or resurge in cyclical epidemics cannot be excluded. Therefore, BRV is an ongoing potential threat to the conservation of M. georgesi, and strict adherence to biosecurity principles is essential to minimise the risk of reintroduction or spread of BRV or other pathogens.

Co-infecting viruses of species Bovine rhinitis B virus (Picornaviridae) and Bovine nidovirus 1 (Tobaniviridae) identified for the first time from a post-mortem respiratory sample of a sheep (Ovis aries) in Hungary.

In this study, a picornavirus and a nidovirus were identified from a single available nasopharyngeal swab (NPS) sample of a freshly deceased sheep, as the only vertebrate viruses found with viral metagenomics and next-generation sequencing methods. The sample was originated from a mixed feedlot farm in Hungary where sheep and cattle were held together but in separate stalls. Most of the sheep had respiratory signs (coughing and increased respiratory effort) at the time of sampling. Other NPS were not, but additional enteric samples were collected from sheep (n = 27) and cattle (n = 11) of the same farm at that time. The complete/nearly complete genomes of the identified viruses were determined using RT-PCR and Nanopore (MinION-Flonge) / Dye-terminator sequencing techniques. The results of detailed genomic and phylogenetic analyses indicate that the identified picornavirus most likely belongs to a type 4 genotype of species Bovine rhinitis B virus (BRBV-4, OR885914) of genus Aphthovirus, family Picornaviridae while the ovine nidovirus (OvNV, OR885915) - as a novel variant - could belong to the recently created Bovine nidovirus 1 (BoNV) species of genus Bostovirus, family Tobaniviridae. None of the identified viruses were detectable in the enteric samples using RT-PCR and generic screening primer pairs. Both viruses are well-known respiratory pathogens of cattle, but their presence was not demonstrated before in other animals, like sheep. Furthermore, neither BRBV-4 nor BoNVs were investigated in European cattle and/or sheep flocks, therefore it cannot be determined whether the presence of these viruses in sheep was a result of a single host species switch/spillover event or these viruses are circulating in not just cattle but sheep populations as well. Further studies required to investigate the spread of these viruses in Hungarian and European sheep and cattle populations and to identify their pathogenic potential in sheep.

Serpentoviruses Exhibit Diverse Organization and ORF Composition with Evidence of Recombination.

Serpentoviruses are a subfamily of positive sense RNA viruses in the order Nidovirales, family Tobaniviridae, associated with respiratory disease in multiple clades of reptiles. While the broadest viral diversity is reported from captive pythons, other reptiles, including colubrid snakes, turtles, and lizards of captive and free-ranging origin are also known hosts. To better define serpentoviral diversity, eleven novel serpentovirus genomes were sequenced with an Illumina MiSeq and, when necessary, completed with other Sanger sequencing methods. The novel serpentoviral genomes, along with 57 other previously published serpentovirus genomes, were analyzed alongside four outgroup genomes. Genomic analyses included identifying unique genome templates for each serpentovirus clade, as well as analysis of coded protein composition, potential protein function, protein glycosylation sites, differences in phylogenetic history between open-reading frames, and recombination. Serpentoviral genomes contained diverse protein compositions. In addition to the fundamental structural spike, matrix, and nucleoprotein proteins required for virion formation, serpentovirus genomes also included 20 previously uncharacterized proteins. The uncharacterized proteins were homologous to a number of previously characterized proteins, including enzymes, transcription factors, scaffolding, viral resistance, and apoptosis-related proteins. Evidence for recombination was detected in multiple instances in genomes from both captive and free-ranging snakes. These results show serpentovirus as a diverse clade of viruses with genomes that code for a wide diversity of proteins potentially enhanced by recombination events.

Expression dynamics of the aplysia abyssovirus.

Two recent studies documented the genome of a novel, extremely large (35.9 kb), nidovirus in RNA sequence databases from the marine neural model Aplysia californica. The goal of the present study was to document the distribution and transcriptional dynamics of this virus, Aplysia abyssovirus 1 (AAbV), in maricultured and wild animals. We confirmed previous findings that AAbV RNA is widespread and reaches extraordinary levels in apparently healthy animals. Transmission electron microscopy identified viral replication factories in ciliated gill epithelial cells but not in neurons where viral RNA is most highly expressed. Viral transcripts do not exhibit evidence of discontinuous RNA synthesis as in coronaviruses but are consistent with production of a single leaderless subgenomic RNA, as in the Gill-associated virus of Penaeus monodon. Splicing patterns in chronically infected adults suggested high levels of defective genomes, possibly explaining the lack of obvious disease signs in high viral load animals.

Detection of nidoviruses in samples collected from captive snakes in Europe between 2016 and 2021.

BACKGROUND: Nidoviruses are increasingly detected in various snake species worldwide, but much remains to be learned about their distribution and the factors influencing their epidemiology. METHODS: This retrospective study evaluated the results of routine nidovirus testing, by PCR, of 5210 swab samples from pet snakes from various European countries that were submitted to a commercial veterinary laboratory in Germany between 2016 and 2021. RESULTS: The overall detection rate was 19.96%. However, the detection rate varied significantly depending on the snake species (p < 0.0001), with the highest rate in Indian pythons (Python molurus) (42.24%). Rates also varied depending on the season of sample collection (p < 0.0001), with the highest rate in winter (24.46%), and the country of sample origin (p < 0.0001), with the highest rate in Austria (36.69%). The detection rate also decreased significantly (p = 0.0003) over the 6-year observation period, from 26.43% to 17.64%. LIMITATION: No information on clinical signs was available for most of the sampled snakes. CONCLUSION: The present study supplies new information on the distribution of python nidoviruses (subgenus Roypretovirus) in pet snakes in Europe and indicates a dynamic situation with possible changes in prevalence over time.

The Isolation and Full-Length Transcriptome Sequencing of a Novel Nidovirus and Response of Its Infection in Japanese Flounder (Paralichthys olivaceus).

A novel nidovirus, CSBV Bces-Po19, was isolated from the marine fish, Japanese flounder (Paralichthys olivaceus). The viral genome was 26,597 nucleotides long and shared 98.62% nucleotide identity with CSBV WHQSR4345. PacBio Sequel and Illumina sequencing were used to perform full-length transcriptome sequencing on CSBV Bces-Po19-sensitive (S) and -resistant (R) Japanese flounder. The results of negative staining revealed bacilliform and spherical virions. There were in total 1444 different genes between CSBV Bces-Po19 S and R groups, with 935 being up-regulated and 513 being down-regulated. Metabolism-, immune-, and RNA-related pathways were significantly enriched. Furthermore, CSBV Bces-Po19 infection induced alternative splicing (AS) events in Japanese flounder; the S group had a higher numbers of AS events (12,352) than the R group (11,452). The number of long non-coding RNA (lncRNA) in the S group, on the other hand, was significantly lower than in the R group. In addition to providing valuable information that sheds more light on CSBV Bces-Po19 infection, these research findings provide further clues for CSBV Bces-Po19 prevention and treatment.

Genomic evolution of the Coronaviridae family.

The current outbreak of coronavirus disease-2019 (COVID-19) caused by SARS-CoV-2 poses unparalleled challenges to global public health. SARS-CoV-2 is a Betacoronavirus, one of four genera belonging to the Coronaviridae subfamily Orthocoronavirinae. Coronaviridae, in turn, are members of the order Nidovirales, a group of enveloped, positive-stranded RNA viruses. Here we present a systematic phylogenetic and evolutionary study based on protein domain architecture, encompassing the entire proteomes of all Orthocoronavirinae, as well as other Nidovirales. This analysis has revealed that the genomic evolution of Nidovirales is associated with extensive gains and losses of protein domains. In Orthocoronavirinae, the sections of the genomes that show the largest divergence in protein domains are found in the proteins encoded in the amino-terminal end of the polyprotein (PP1ab), the spike protein (S), and many of the accessory proteins. The diversity among the accessory proteins is particularly striking, as each subgenus possesses a set of accessory proteins that is almost entirely specific to that subgenus. The only notable exception to this is ORF3b, which is present and orthologous over all Alphacoronaviruses. In contrast, the membrane protein (M), envelope small membrane protein (E), nucleoprotein (N), as well as proteins encoded in the central and carboxy-terminal end of PP1ab (such as the 3C-like protease, RNA-dependent RNA polymerase, and Helicase) show stable domain architectures across all Orthocoronavirinae. This comprehensive analysis of the Coronaviridae domain architecture has important implication for efforts to develop broadly cross-protective coronavirus vaccines.

Readdressing the genetic diversity and taxonomy of the Mesoniviridae family, as well as its relationships with other nidoviruses and putative mesonivirus-like viral sequences.

Research on the recently established Mesoniviridae family (Order Nidovirales), RNA genome insect-specific viruses, has been steadily growing in the last decade. However, after the last detailed phylogenetic characterization of mesoniviruses in 2014, numerous new sequences, even in organisms other than mosquitos, have been identified and characterized. In this study, we analyzed nucleotide and protein sequences of mesoniviruses with a wide range of molecular tools including genetic distance, Shannon entropy, selective pressure analysis, polymorphism identification, principal coordinate analysis, likelihood mapping and phylodynamic reconstruction. We also sought to revaluate new mesoniviruses sequence positions within the family, proposing a taxonomic revision. The different sub-lineages of mosquito mesoniviruses sequences presented low sequence diversity and entropy, with incongruences to the existing taxonomy being found after an extensive phylogenetic characterization. High sequence discrepancy and differences in genome organization were found between mosquito mesoniviruses and other mesoniviruses, so their future classification, as other meso-like viruses that are found in other organisms, should be approached with caution. No evidence of frequent recombination was found, and mesonivirus genomes seem to evolve under strong purifying selection. Insufficient data by root-to-tip analysis did not yet allow for an adequate phylogeographic reconstruction.

Protein Nucleotidylylation in +ssRNA Viruses.

Nucleotidylylation is a post-transcriptional modification important for replication in the picornavirus supergroup of RNA viruses, including members of the Caliciviridae, Coronaviridae, Picornaviridae and Potyviridae virus families. This modification occurs when the RNA-dependent RNA polymerase (RdRp) attaches one or more nucleotides to a target protein through a nucleotidyl-transferase reaction. The most characterized nucleotidylylation target is VPg (viral protein genome-linked), a protein linked to the 5' end of the genome in Caliciviridae, Picornaviridae and Potyviridae. The nucleotidylylation of VPg by RdRp is a critical step for the VPg protein to act as a primer for genome replication and, in Caliciviridae and Potyviridae, for the initiation of translation. In contrast, Coronaviridae do not express a VPg protein, but the nucleotidylylation of proteins involved in replication initiation is critical for genome replication. Furthermore, the RdRp proteins of the viruses that perform nucleotidylylation are themselves nucleotidylylated, and in the case of coronavirus, this has been shown to be essential for viral replication. This review focuses on nucleotidylylation within the picornavirus supergroup of viruses, including the proteins that are modified, what is known about the nucleotidylylation process and the roles that these modifications have in the viral life cycle.

Structural analysis of viral ExoN domains reveals polyphyletic hijacking events.

Nidoviruses and arenaviruses are the only known RNA viruses encoding a 3'-5' exonuclease domain (ExoN). The proofreading activity of the ExoN domain has played a key role in the growth of nidoviral genomes, while in arenaviruses this domain partakes in the suppression of the host innate immune signaling. Sequence and structural homology analyses suggest that these proteins have been hijacked from cellular hosts many times. Analysis of the available nidoviral ExoN sequences reveals a high conservation level comparable to that of the viral RNA-dependent RNA polymerases (RdRp), which are the most conserved viral proteins. Two highly preserved zinc fingers are present in all nidoviral exonucleases, while in the arenaviral protein only one zinc finger can be identified. This is in sharp contrast with the reported lack of zinc fingers in cellular ExoNs, and opens the possibility of therapeutic strategies in the struggle against COVID-19.

A nidovirus perspective on SARS-CoV-2.

Two pandemics of respiratory distress diseases associated with zoonotic introductions of the species Severe acute respiratory syndrome-related coronavirus in the human population during 21st century raised unprecedented interest in coronavirus research and assigned it unseen urgency. The two viruses responsible for the outbreaks, SARS-CoV and SARS-CoV-2, respectively, are in the spotlight, and SARS-CoV-2 is the focus of the current fast-paced research. Its foundation was laid down by studies of many corona- and related viruses that collectively form the vast order Nidovirales. Comparative genomics of nidoviruses played a key role in this advancement over more than 30 years. It facilitated the transfer of knowledge from characterized to newly identified viruses, including SARS-CoV and SARS-CoV-2, as well as contributed to the dissection of the nidovirus proteome and identification of patterns of variations between different taxonomic groups, from species to families. This review revisits selected cases of protein conservation and variation that define nidoviruses, illustrates the remarkable plasticity of the proteome during nidovirus adaptation, and asks questions at the interface of the proteome and processes that are vital for nidovirus reproduction and could inform the ongoing research of SARS-CoV-2.

Serpentovirus (Nidovirus) and Orthoreovirus Coinfection in Captive Veiled Chameleons (Chamaeleo calyptratus) with Respiratory Disease.

Serpentoviruses are an emerging group of nidoviruses known to cause respiratory disease in snakes and have been associated with disease in other non-avian reptile species (lizards and turtles). This study describes multiple episodes of respiratory disease-associated mortalities in a collection of juvenile veiled chameleons (Chamaeleo calyptratus). Histopathologic lesions included rhinitis and interstitial pneumonia with epithelial proliferation and abundant mucus. Metagenomic sequencing detected coinfection with two novel serpentoviruses and a novel orthoreovirus. Veiled chameleon serpentoviruses are most closely related to serpentoviruses identified in snakes, lizards, and turtles (approximately 40-50% nucleotide and amino acid identity of ORF1b). Veiled chameleon orthoreovirus is most closely related to reptilian orthoreoviruses identified in snakes (approximately 80-90% nucleotide and amino acid identity of the RNA-dependent RNA polymerase). A high prevalence of serpentovirus infection (>80%) was found in clinically healthy subadult and adult veiled chameleons, suggesting the potential for chronic subclinical carriers. Juvenile veiled chameleons typically exhibited a more rapid progression compared to subadults and adults, indicating a possible age association with morbidity and mortality. This is the first description of a serpentovirus infection in any chameleon species. A causal relationship between serpentovirus infection and respiratory disease in chameleons is suspected. The significance of orthoreovirus coinfection remains unknown.

Genetic, Morphological and Antigenic Relationships between Mesonivirus Isolates from Australian Mosquitoes and Evidence for Their Horizontal Transmission.

The Mesoniviridae are a newly assigned family of viruses in the order Nidovirales. Unlike other nidoviruses, which include the Coronaviridae, mesoniviruses are restricted to mosquito hosts and do not infect vertebrate cells. To date there is little information on the morphological and antigenic characteristics of this new group of viruses and a dearth of mesonivirus-specific research tools. In this study we determined the genetic relationships of recent Australian isolates of Alphamesonivirus 4 (Casuarina virus-CASV) and Alphamesonivirus 1 (Nam Dinh virus-NDiV), obtained from multiple mosquito species. Australian isolates of NDiV showed high-level similarity to the prototype NDiV isolate from Vietnam (99% nucleotide (nt) and amino acid (aa) identity). Isolates of CASV from Central Queensland were genetically very similar to the prototype virus from Darwin (95-96% nt and 91-92% aa identity). Electron microscopy studies demonstrated that virion diameter (≈80 nm) and spike length (≈10 nm) were similar for both viruses. Monoclonal antibodies specific to CASV and NDiV revealed a close antigenic relationship between the two viruses with 13/34 mAbs recognising both viruses. We also detected NDiV RNA on honey-soaked nucleic acid preservation cards fed on by wild mosquitoes supporting a possible mechanism of horizontal transmission between insects in nature.

High-Quality Resolution of the Outbreak-Related Zika Virus Genome and Discovery of New Viruses Using Ion Torrent-Based Metatranscriptomics.

Arboviruses, including the Zika virus, have recently emerged as one of the most important threats to human health. The use of metagenomics-based approaches has already proven valuable to aid surveillance of arboviral infections, and the ability to reconstruct complete viral genomes from metatranscriptomics data is key to the development of new control strategies for these diseases. Herein, we used RNA-based metatranscriptomics associated with Ion Torrent deep sequencing to allow for the high-quality reconstitution of an outbreak-related Zika virus (ZIKV) genome (10,739 nt), with extended 5'-UTR and 3'-UTR regions, using a newly-implemented bioinformatics approach. Besides allowing for the assembly of one of the largest complete ZIKV genomes to date, our strategy also yielded high-quality complete genomes of two arthropod-infecting viruses co-infecting C6/36 cell lines, namely: Alphamesonivirus 1 strain Salvador (20,194 nt) and Aedes albopictus totivirus-like (4618 nt); the latter likely represents a new viral species. Altogether, our results demonstrate that our bioinformatics approach associated with Ion Torrent sequencing allows for the high-quality reconstruction of known and unknown viral genomes, overcoming the main limitation of RNA deep sequencing for virus identification.

A case report of reptile-associated nidovirus (serpentovirus) in a ball python (Python regius) in Taiwan.

Reptile-associated nidoviruses (serpentoviruses) have been reported to cause proliferative interstitial pneumonia in pythons and other reptile species. A captive, younger than 2 years old, intact female ball python (Python regius) showed increased oral mucus, wheezing, and audible breathing with weight loss. Gross and microscopic examination revealed large amounts of mucus in the esophagus and proliferative interstitial pneumonia. Serpentovirus genes were detected from the lung tissues by polymerase chain reaction. The current serpentoviruses was phylogenetically grouped with the serpentovirus previously identified in the US. No case of serpentovirus infection has been reported in Asia. The present report provides information of complete genome sequence and global distribution of serpentovirus.

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2.

The present outbreak of a coronavirus-associated acute respiratory disease called coronavirus disease 19 (COVID-19) is the third documented spillover of an animal coronavirus to humans in only two decades that has resulted in a major epidemic. The Coronaviridae Study Group (CSG) of the International Committee on Taxonomy of Viruses, which is responsible for developing the classification of viruses and taxon nomenclature of the family Coronaviridae, has assessed the placement of the human pathogen, tentatively named 2019-nCoV, within the Coronaviridae. Based on phylogeny, taxonomy and established practice, the CSG recognizes this virus as forming a sister clade to the prototype human and bat severe acute respiratory syndrome coronaviruses (SARS-CoVs) of the species Severe acute respiratory syndrome-related coronavirus, and designates it as SARS-CoV-2. In order to facilitate communication, the CSG proposes to use the following naming convention for individual isolates: SARS-CoV-2/host/location/isolate/date. While the full spectrum of clinical manifestations associated with SARS-CoV-2 infections in humans remains to be determined, the independent zoonotic transmission of SARS-CoV and SARS-CoV-2 highlights the need for studying viruses at the species level to complement research focused on individual pathogenic viruses of immediate significance. This will improve our understanding of virus­host interactions in an ever-changing environment and enhance our preparedness for future outbreaks.

Investigations into the presence of nidoviruses in pythons.

BACKGROUND: Pneumonia and stomatitis represent severe and often fatal diseases in different captive snakes. Apart from bacterial infections, paramyxo-, adeno-, reo- and arenaviruses cause these diseases. In 2014, new viruses emerged as the cause of pneumonia in pythons. In a few publications, nidoviruses have been reported in association with pneumonia in ball pythons and a tiger python. The viruses were found using new sequencing methods from the organ tissue of dead animals. METHODS: Severe pneumonia and stomatitis resulted in a high mortality rate in a captive breeding collection of green tree pythons. Unbiased deep sequencing lead to the detection of nidoviral sequences. A developed RT-qPCR was used to confirm the metagenome results and to determine the importance of this virus. A total of 1554 different boid snakes, including animals suffering from respiratory diseases as well as healthy controls, were screened for nidoviruses. Furthermore, in addition to two full-length sequences, partial sequences were generated from different snake species. RESULTS: The assembled full-length snake nidovirus genomes share only an overall genome sequence identity of less than 66.9% to other published snake nidoviruses and new partial sequences vary between 99.89 and 79.4%. Highest viral loads were detected in lung samples. The snake nidovirus was not only present in diseased animals, but also in snakes showing no typical clinical signs. CONCLUSIONS: Our findings further highlight the possible importance of snake nidoviruses in respiratory diseases and proof multiple circulating strains with varying disease potential. Nidovirus detection in clinical healthy individuals might represent testing during the incubation period or reconvalescence. Our investigations show new aspects of nidovirus infections in pythons. Nidoviruses should be included in routine diagnostic workup of diseased reptiles.

Endangered wild salmon infected by newly discovered viruses.

The collapse of iconic, keystone populations of sockeye (Oncorhynchus nerka) and Chinook (Oncorhynchus tshawytscha) salmon in the Northeast Pacific is of great concern. It is thought that infectious disease may contribute to declines, but little is known about viruses endemic to Pacific salmon. Metatranscriptomic sequencing and surveillance of dead and moribund cultured Chinook salmon revealed a novel arenavirus, reovirus and nidovirus. Sequencing revealed two different arenavirus variants which each infect wild Chinook and sockeye salmon. In situ hybridisation localised arenavirus mostly to blood cells. Population surveys of >6000 wild juvenile Chinook and sockeye salmon showed divergent distributions of viruses, implying different epidemiological processes. The discovery in dead and dying farmed salmon of previously unrecognised viruses that are also widely distributed in wild salmon, emphasizes the potential role that viral disease may play in the population dynamics of wild fish stocks, and the threat that these viruses may pose to aquaculture.