A highly divergent picornavirus in a marine mammal

Structure and function of type IV IRES in picornaviruses: a systematic review

The Picornaviridae is a family of icosahedral viruses with single-stranded, highly diverse positive-sense RNA genomes. Virions consist of a capsid, without envelope, surrounding a core of RNA genome. A typical genome of picornavirus harbors a well-conserved and highly structured RNA element known as the internal ribosome entry site (IRES), functionally essential for viral replication and protein translation. Based on differences in their structures and mechanisms of action, picornaviral IRESs have been categorized into five types: type I, II, III, IV, and V. Compared with the type IV IRES, the others not only are structurally complicated, but also involve multiple initiation factors for triggering protein translation. The type IV IRES, often referred to as hepatitis C virus (HCV)-like IRES due to its structural resemblance to the HCV IRES, exhibits a simpler and more compact structure than those of the other four. The increasing identification of picornaviruses with the type IV IRES suggests that this IRES type seems to reveal strong retention and adaptation in terms of viral evolution. Here, we systematically reviewed structural features and biological functions of the type IV IRES in picornaviruses. A comprehensive understanding of the roles of type IV IRESs will contribute to elucidating the replication mechanism and pathogenesis of picornaviruses.

The 28S rRNA RT-qPCR assay for host depletion evaluation to enhance avian virus detection in Illumina and Nanopore sequencing

Abundant host and bacterial sequences can obscure the detection of less prevalent viruses in untargeted next-generation sequencing (NGS). Efficient removal of these non-targeted sequences is vital for accurate viral detection. This study presents a novel 28S ribosomal RNA (rRNA) RT-qPCR assay designed to assess the efficiency of avian rRNA depletion before conducting costly NGS for the detection of avian RNA viruses. The comprehensive evaluation of this 28S-test focuses on substituting DNase I with alternative DNases in our established depletion protocols and finetuning essential parameters for reliable host rRNA depletion. To validate the effectiveness of the 28S-test, we compared its performance with NGS results obtained from both Illumina and Nanopore sequencing platforms. This evaluation utilized swab samples from chickens infected with highly pathogenic avian influenza virus, subjected to established and modified depletion protocols. Both methods significantly reduced host rRNA levels, but using the alternative DNase had superior performance. Additionally, utilizing the 28S-test, we explored cost- and time-effective strategies, such as reduced probe concentrations and other alternative DNase usage, assessed the impact of filtration pre-treatment, and evaluated various experimental parameters to further optimize the depletion protocol. Our findings underscore the value of the 28S-test in optimizing depletion methods for advancing improvements in avian disease research through NGS.

A widespread picornavirus affects the hemocytes of the noble pen shell (Pinna nobilis), leading to its immunosuppression

Introduction

The widespread mass mortality of the noble pen shell (Pinna nobilis) has occurred in several Mediterranean countries in the past 7 years. Single-stranded RNA viruses affecting immune cells and leading to immune dysfunction have been widely reported in human and animal species. Here, we present data linking P. nobilis mass mortality events (MMEs) to hemocyte picornavirus (PV) infection. This study was performed on specimens from wild and captive populations.

Methods

We sampled P. nobilis from two regions of Spain [Catalonia (24 animals) and Murcia (four animals)] and one region in Italy [Venice (6 animals)]. Each of them were analyzed using transmission electron microscopy (TEM) to describe the morphology and self-assembly of virions. Illumina sequencing coupled to qPCR was performed to describe the identified virus and part of its genome.

Results and discussion

In 100% of our samples, ultrastructure revealed the presence of a virus (20 nm diameter) capable of replicating within granulocytes and hyalinocytes, leading to the accumulation of complex vesicles of different dimensions within the cytoplasm. As the PV infection progressed, dead hemocytes, infectious exosomes, and budding of extracellular vesicles were visible, along with endocytic vesicles entering other cells. The THC (total hemocyte count) values observed in both captive (eight animals) (3.5 × 104-1.60 × 105 ml-1 cells) and wild animals (14 samples) (1.90-2.42 × 105 ml-1 cells) were lower than those reported before MMEs. Sequencing of P. nobilis (six animals) hemocyte cDNA libraries revealed the presence of two main sequences of Picornavirales, family Marnaviridae. The highest number of reads belonged to animals that exhibited active replication phases and abundant viral particles from transmission electron microscopy (TEM) observations. These sequences correspond to the genus Sogarnavirus-a picornavirus identified in the marine diatom Chaetoceros tenuissimus (named C. tenuissimus RNA virus type II). Real-time PCR performed on the two most abundant RNA viruses previously identified by in silico analysis revealed positive results only for sequences similar to the C. tenuissimus RNA virus. These results may not conclusively identify picornavirus in noble pen shell hemocytes; therefore, further study is required. Our findings suggest that picornavirus infection likely causes immunosuppression, making individuals prone to opportunistic infections, which is a potential cause for the MMEs observed in the Mediterranean.

Rapid whole genome sequencing methods for RNA viruses

RNA viruses are the etiological agents of many infectious diseases. Since RNA viruses are error-prone during genome replication, rapid, accurate and economical whole RNA viral genome sequence determination is highly demanded. Next-generation sequencing (NGS) techniques perform whole viral genome sequencing due to their high-throughput sequencing capacity. However, the NGS techniques involve a significant burden for sample preparation. Since to generate complete viral genome coverage, genomic nucleic acid enrichment is required by reverse transcription PCR using virus-specific primers or by viral particle concentration. Furthermore, conventional NGS techniques cannot determine the 5' and 3' terminal sequences of the RNA viral genome. Therefore, the terminal sequences are determined one by one using rapid amplification of cDNA ends (RACE). However, since some RNA viruses have segmented genomes, the burden of the determination using RACE is proportional to the number of segments. To date, there is only one study attempting whole genome sequencing of multiple RNA viruses without using above mentioned methods, but the generated sequences' accuracy compared to the reference sequences was up to 97% and did not reach 100% due to the low read depth. Hence, we established novel methods, named PCR-NGS and RCA-NGS, that were optimized for an NGS machine, MinION. These methods do not require nucleic acid amplification with virus-specific PCR primers, physical viral particle enrichment, and RACE. These methods enable whole RNA viral genome sequencing by combining the following techniques: (1) removal of unwanted DNA and RNA other than the RNA viral genome by nuclease treatment; (2) the terminal of viral genome sequence determination by barcoded linkers ligation; (3) amplification of the viral genomic cDNA using ligated linker sequences-specific PCR or an isothermal DNA amplification technique, such as rolling circle amplification (RCA). The established method was evaluated using isolated RNA viruses with single-stranded, double-stranded, positive-stranded, negative-stranded, non-segmented or multi-segmented genomes. As a result, all the viral genome sequences could be determined with 100% accuracy, and these mean read depths were greater than 2,500×, at least using either of the methods. This method should allow for easy and economical determination of accurate RNA viral genomes.

Expansion of viral genomes with viral protein genome linked copies

Virus protein-linked genome (VPg) proteins are required for replication. VPgs are duplicated in a subset of RNA viruses however their roles are not fully understood and the extent of viral genomes containing VPg copies has not been investigated in detail. Here, we generated a novel bioinformatics approach to identify VPg sequences in viral genomes using hidden Markov models (HMM) based on alignments of dicistrovirus VPg sequences. From metagenomic datasets of dicistrovirus genomes, we identified 717 dicistrovirus genomes containing VPgs ranging from a single copy to 8 tandem copies. The VPgs are classified into nine distinct types based on their sequence and length. The VPg types but not VPg numbers per viral genome followed specific virus clades, thus suggesting VPgs co-evolved with viral genomes. We also identified VPg duplications in aquamavirus and mosavirus genomes. This study greatly expands the number of viral genomes that contain VPg copies and indicates that duplicated viral sequences are more widespread than anticipated.

An Emerging Duck Egg-Reducing Syndrome Caused by a Novel Picornavirus Containing Seven Putative 2A Peptides

Since 2016, frequent outbreaks of egg-reducing syndromes caused by an unknown virus in duck farms have resulted in huge economic losses in China. The causative virus was isolated and identified as a novel species in Avihepatovirus of the picornavirus family according to the current guidelines of the International Committee on Taxonomy of Viruses (ICVT), and was named the duck egg-reducing syndrome virus (DERSV). The DERSV was most closely related to wild duck avihepatovirus-like virus (WDALV) with 64.0%, 76.8%, 77.5%, and 70.7% of amino acid identities of P1, 2C, 3C, and 3D proteins, respectively. The DERSV had a typical picornavirus-like genomic structure, but with the longest 2A region in the reported picornaviruses so far. Importantly, the clinical symptoms were successfully observed by artificially infecting ducks with DERSV, even in the contact exposed ducks, which suggested that DERSV transmitted among ducks by direct contact. The antibody levels of DERSV were correlated with the emergence of the egg-reducing syndromes in ducks in field. These results indicate that DERSV is a novel emerging picornavirus causing egg-reducing syndrome in ducks.

Virome in Fecal Samples From Wild Giant Pandas (Ailuropoda Melanoleuca)

The giant panda (Ailuropoda melanoleuca) is one of the most endangered mammals in the world; anthropogenic habitat loss and poaching still threaten the survival of wild pandas. Viral infection has become one of the potential threats to the health of these animals, but the available information related to these infections is still limited. In order to detect possible vertebrate viruses, the virome in the fecal samples of seven wild giant pandas from Qinling Mountains was investigated by using the method of viral metagenomics. From the fecal virome of wild giant pandas, we determined six nearly complete genomes belonging to the order Picornavirales, two of which may be qualified as a novel virus family or genus. In addition, four complete genomes belonging to the Genomoviridae family were also fully characterized. This virological investigation has increased our understanding of the gut viral community in giant pandas. Whether these viruses detected in fecal samples can really infect giant panda needs further research.

Virus-derived sequences from the transcriptomes of two snail vectors of schistosomiasis, Biomphalaria pfeifferi and Bulinus globosus from Kenya

Schistosomiasis, which infects more than 230 million people, is vectored by freshwater snails. We identified viral sequences in the transcriptomes of Biomphalaria pfeifferi (BP) and Bulinus globosus (BuG), two of the world's most important schistosomiasis vectors in Africa. Sequences from 26 snails generated using Illumina Hi-Seq or 454 sequencing were assembled using Trinity and CAP3 and putative virus sequences were identified using a bioinformatics pipeline. Phylogenetic analyses were performed using viral RNA-dependent RNA polymerase and coat protein sequences to establish relatedness between virus sequences identified and those of known viruses. Viral sequences were identified from the entire snail holobiont, including symbionts, ingested material and organisms passively associated with the snails. Sequences derived from more than 17 different viruses were found including five near full-length genomes, most of which were small RNA viruses with positive sense RNA genomes (i.e., picorna-like viruses) and some of which are likely derived from adherent or ingested diatoms. Based on phylogenetic analysis, five of these viruses (including BPV2 and BuGV2) along with four Biomphalaria glabrata viruses reported previously, cluster with known invertebrate viruses and are putative viruses of snails. The presence of RNA sequences derived from four of these novel viruses in samples was confirmed. Identification of the genome sequences of candidate snail viruses provides a first step toward characterization of additional gastropod viruses, including from species of biomedical significance.

Identification and genome analysis of a novel picornavirus from captive belugas (Delphinapterus leucas) in China

The discovery of new viruses is important for predicting their potential threats to the health of humans and other animals. A novel picornavirus was identified from oral, throat, and anal swab samples collected from belugas (Delphinapterus leucas), from Dalian Sun Asia Tourism Holding Co., China, between January and December 2018, using a metagenomics approach. The genome of this novel PicoV-HMU-1 strain was 8197 nucleotides (nt) in length, with a open reading frame (from 1091 to 8074 nt) that encoded a polyprotein precursor of 2328 amino acids. Moreover, the genomic length and GC content of PicoV-HMU-1 were within the ranges found in other picornaviruses, and the genome organization was also similar. Nevertheless, PicoV-HMU-1 had a lower amino acid identity and distinct host species compared with other members of the Picornaviridae family. Phylogenetic trees were constructed based on the P1 and 3D amino acid sequences of PicoV-HMU-1 along with representative members of the Picornaviridae family, which showed that PicoV-HMU-1 was related to unclassified bat picornaviruses groups. These findings suggest that the PicoV-HMU-1 strain represents a potentially novel genus of picornavirus. These data can enhance our understanding of the picornavirus genetic diversity and evolution.

Epidemiology and genetic characteristics of murine kobuvirus from faecal samples of Rattus losea, Rattus tanezumi and Rattus norvegicus in southern China

Recently, murine kobuvirus (MuKV), a novel member of the family Picornaviridae, was identified in faecal samples of Rattus norvegicus in China. The limited information on the circulation of MuKV in other murine rodent species prompted us to investigate its prevalence and conduct a genetic characterization of MuKV in Rattus losea, Rattus tanezumi and Rattus norvegicus in China. Between 2015 and 2017, 243 faecal samples of these three murine rodent species from three regions in southern China were screened for the presence of MuKV. The overall prevalence was 23.0% (56/243). Three complete MuKV polyprotein sequences were acquired, and the genome organization was determined. Phylogenetic analyses suggested that our sequences were closely related to Chinese strains and belong to the species Aichivirus A in the genus Kobuvirus. Additional studies are required to understand the true prevalence of MuKV in murine rodent populations in China.

Picornaviruses: A View from 3A

Picornaviruses are comprised of a positive-sense RNA genome surrounded by a protein shell (or capsid). They are ubiquitous in vertebrates and cause a wide range of important human and animal diseases. The genome encodes a single large polyprotein that is processed to structural (capsid) and non-structural proteins. The non-structural proteins have key functions within the viral replication complex. Some, such as 3D^pol (the RNA dependent RNA polymerase) have conserved functions and participate directly in replicating the viral genome, whereas others, such as 3A, have accessory roles. The 3A proteins are highly divergent across the Picornaviridae and have specific roles both within and outside of the replication complex, which differ between the different genera. These roles include subverting host proteins to generate replication organelles and inhibition of cellular functions (such as protein secretion) to influence virus replication efficiency and the host response to infection. In addition, 3A proteins are associated with the determination of host range. However, recent observations have challenged some of the roles assigned to 3A and suggest that other viral proteins may carry them out. In this review, we revisit the roles of 3A in the picornavirus life cycle. The 3AB precursor and mature 3A have distinct functions during viral replication and, therefore, we have also included discussion of some of the roles assigned to 3AB.

Persistent Detection of Cosavirus and Saffold Cardiovirus in Riachuelo River, Argentina

Cosaviruses (CoSV) and Saffold cardiovirus (SAFV) are novel members of the Picornaviridae family. The Matanza-Riachuelo river basin covers a total area of 2200 km² with approximately 60 km long. Its last section is called Riachuelo River. The aim of this study was to describe the circulation of both picornaviruses and their relationship with the environmental situation of the Riachuelo River using 274 samples collected from 2005 to 2015. CoSV and SAFV were investigated in samples available by two periods: 2005-2006 and 2014-2015 (103 and 101, respectively). Physicochemical and bacteriological parameters confirmed very high levels of human fecal contamination during the 11 years evaluated. CoSV was detected in 85.7% (66/77) and 65.4% (17/26) of the samples collected in 2005-2006 and 2014-2015 periods, respectively. Species A and D were identified, the first one being widely predominant: 74.1% (20/27) and 75.0% (3/4) in both periods. SAFV virus was detected in 47.1% (32/68) and 52.6% (10/19) in periods 2005-2006 and 2014-2015, respectively. SAFV-6 was the most identified genotype in the entire study, while SAFV-3 was predominant in 2005-2006. The contribution of genotypes 1, 2, 4 and 8 was minor. The high prevalence of CoSV and SAFV suggests that both viruses have been circulating in Argentina at least since 2005. Our results show that a watercourse with high rates of human fecal contamination can become a persistent source of new viruses which capacity to produce human diseases is unknown.

Genomic Characterization of Picornaviruses Isolated From Ribbon (Histriophoca fasciata) and Harbor (Phoca vitulina) Seals

The seal picornavirus 1, species Aquamavirus A, is currently the only recognized member of the genus Aquamavirus within the family Picornaviridae. The bear picornavirus 1 was recently proposed as the second species in the genus under the name aquamavirus B. Herein, we determined the complete genomes of two novel pinniped picornaviruses, the harbor seal picornavirus (HsPV) and the ribbon seal picornavirus (RsPV). The HsPV and the RsPV were isolated in Vero.DogSLAMtag cells from samples collected from stranded harbor (Phoca vitulina) and ribbon (Histriophoca fasciata) seals. RsPV-infected Vero.DogSLAMtag cells displaying extensive cytopathic effects were processed for transmission electron microscopy and revealed non-enveloped viral particles aggregated into paracrystalline arrays in the cytoplasm. A next-generation sequencing approach was used to recover the complete genomes of the HsPV and the RsPV (6,709 and 6,683 bp, respectively). Phylogenetic and genetic analyses supported the HsPV and the RsPV as members of the Aquamavirus genus. Based on these results, RsPV represents a novel strain of Aquamavirus A, while the HsPV is a novel strain of the proposed species aquamavirus B. These discoveries provide information on the evolutionary relationships and ultrastructure of aquamaviruses and expands the known host range of those viruses. Our results underscore the importance of the application of classical virology and pathology techniques coupled with high-throughput sequencing technologies for the discovery and characterization of pathogens in wild marine mammals.

Next-generation sequencing in the diagnosis of viral encephalitis: sensitivity and clinical limitations

Identification of pathogens causing viral encephalitis remains challenging, and in over 50% of cases the etiologic factor remains undetermined. Next-generation sequencing (NGS) based metagenomics has been successfully used to detect novel and rare infections, but its value for routine diagnosis of encephalitis remains unclear. The aim of the present study was to determine the sensitivity of shotgun metagenomic sequencing protocols, which include preamplification, and testing it against cerebrospinal fluid (CSF) samples from encephalitis patients. For sensitivity testing HIV and HBV positive sera were serially diluted in CSF from an uninfected patient. NGS repeatedly detected HIV and HBV sequences present at concentrations from 10⁵ to 10² and from 10⁵ to 10 viral copies/reaction, respectively. However, when the same protocols were applied to RT-PCR/PCR positive CSF samples from 6 patients with enteroviral encephalitis (median viral load 47 copies/ml) and 15 patients with HSV, CMV or VZV encephalitis (median viral load 148 copies/ml), only 7 (28.6%) were identified as positive. In conclusions, while NGS has the advantage of being able to identify a wide range of potential pathogens it seems to be less sensitive compared to the standard amplification-based assays in the diagnosis of encephalitis, where low viral loads are common.

The Husavirus Posa-Like Viruses in China, and a New Group of Picornavirales

Novel posa-like viral genomes were first identified in swine fecal samples using metagenomics and were designated as unclassified viruses in the order Picornavirales. In the present study, nine husavirus strains were identified in China. Their genomes share 94.1–99.9% similarity, and alignment of these nine husavirus strains identified 697 nucleotide polymorphism sites across their full-length genomes. These nine strains were directly clustered with the Husavirus 1 lineage, and their genomic arrangement showed similar characteristics. These posa-like viruses have undergone a complex evolutionary process, and have a wide geographic distribution, complex host spectrum, deep phylogenetic divergence, and diverse genomic organizations. The clade of posa-like viruses forms a single group, which is evolutionarily distinct from other known families and could represent a distinct family within the Picornavirales. The genomic arrangement of Picornavirales and the new posa-like viruses are different, whereas the posa-like viruses have genomic modules similar to the families Dicistroviridae and Marnaviridae. The present study provides valuable genetic evidence of husaviruses in China, and clarifies the phylogenetic dynamics and the evolutionary characteristics of Picornavirales.

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.

Asian black bear (Ursus thibetanus) picornavirus related to seal aquamavirus A

The complete genome of a bear picornavirus 1 (BePV-1) in the viscera of an Asian black bear (Ursus thibetanus) from China was characterized using viral metagenomics and RT-PCR/Sanger sequencing. The genome of BePV1 is 6703 nt long, contains a type-IV IRES 5'UTR with the '8-like' motif, encodes a 2053-aa-long polyprotein showing a 3-4-4 organization pattern and two 2A genes. BePV-1 showed the highest overall genome nucleotide sequence identity of 71.7% to a picornavirus genome from an Arctic ringed seal (Phoca hispida) from Canada, classified as a member of the species Aquamavirus A, currently the only one in the genus Aquamavirus. Phylogenetic and genetic distance analyses of P1 and 3D indicated that Asian bear picornavirus (aquamavirus B) represents the second sequenced member of the genus Aquamavirus.

Limited Correlation of Shotgun Metagenomics Following Host Depletion and Routine Diagnostics for Viruses and Bacteria in Low Concentrated Surrogate and Clinical Samples

The etiologic cause of encephalitis, meningitis or meningo-encephalitis is unknown in up to 70% of cases. Clinical shotgun metagenomics combined with host depletion is a promising technique to identify infectious etiologies of central nervous system (CNS) infections. We developed a straightforward eukaryotic host nucleic acid depletion method that preserves intact viruses and bacteria for subsequent shotgun metagenomics screening of clinical samples, focusing on cerebrospinal fluid (CSF). A surrogate CSF sample for a CNS infection paradigm was used to evaluate the proposed depletion method consisting of selective host cell lysis, followed by enzymatic degradation of the liberated genomic DNA for final depletion with paramagnetic beads. Extractives were subjected to reverse transcription, followed by whole genome amplification and next generation sequencing. The effectiveness of the host depletion method was demonstrated in surrogate CSF samples spiked with three 1:100 dilutions of Influenza A H3N2 virus (qPCR Ct-values 20.7, 28.8, >42/negative). Compared to the native samples, host depletion increased the amount of the virus subtype reads by factor 7127 and 132, respectively, while in the qPCR negative sample zero vs. 31 (1.4E-4 %) virus subtype reads were detected (native vs. depleted). The workflow was applied to thirteen CSF samples of patients with meningo-/encephalitis (two bacterial, eleven viral etiologies), a serum of an Andes virus infection and a nose swab of a common cold patient. Unlike surrogate samples, host depletion of the thirteen human CSF samples and the nose swab did not result in more reads indicating presence of damaged pathogens due to, e.g., host immune response. Nevertheless, previously diagnosed pathogens in the human CSF samples (six viruses, two bacteria), the serum, and the nose swab (Human rhinovirus A31) were detected in the depleted and/or the native samples. Unbiased evaluation of the taxonomic profiles supported the diagnosed pathogen in two native CSF samples and the native and depleted serum and nose swab, while detecting various contaminations that interfered with pathogen identification at low concentration levels. In summary, damaged pathogens and contaminations complicated analysis and interpretation of clinical shotgun metagenomics data. Still, proper consideration of these issues may enable future application of metagenomics for clinical diagnostics.

Detection by Direct Next Generation Sequencing Analysis of Emerging Enterovirus D68 and C109 Strains in an Environmental Sample From Scotland

Background: Human enteroviruses (EVs) have been linked with severe disease and syndromes as varied as acute respiratory illness, myocarditis, and flaccid paralysis. With global polio eradication on sight the focus of clinical investigations has expanded to the identification of other EV serotypes associated with severe neurological conditions such as EV-D68, responsible for large outbreaks in 2014 and 2016 that spread worldwide and were related with severe respiratory disease leading to acute myelitis in some cases. New EV serotypes with epidemic potential continue to emerge such as EV-C104, EV-C105, EV-C109, and EV-C117 identified in respiratory samples in recent years.

Methods: We used a next generation sequencing (NGS) approach to detect multiple EV serotypes directly in a sewage concentrate from Glasgow (Scotland, United Kingdom) generating whole-capsid nucleotide sequences that were compared to sequences of cell culture isolates from this sewage sample and clinical EV isolates from GenBank.

Results: Thirteen different serotypes belonging to all four A, B, C, and D EV species were identified in the sewage concentrate. EV strains closely related to EV-D68 epidemic isolates of B3 lineage reported in the United States and Europe in 2016 and to EV-C109 respiratory isolates found in Denmark and Netherlands in 2015 were identified.

Conclusion: Environmental surveillance (ES) can effectively detect EV circulation in human populations. The use of NGS for ES can help overcoming the limitations of traditional cell culture and sequencing methods, which are selective and biased, and can contribute to the early detection and assessment of spread of emerging EV pathogens.

Identification and whole-genome characterization of a recombinant Enterovirus B69 isolated from a patient with Acute Flaccid Paralysis in Niger, 2015

Enterovirus B69 (EV-B69) is a rarely reported type and till date, only the full-length genome sequence of the prototype strain is available. Besides the prototype strain, only limited VP1 sequences of this virus from Africa and India are available in GenBank. In this study, we analyzed the full-length genome sequence of an EV-B69 strain recovered from a patient with acute flaccid paralysis in Niger. Compared with the EV-B69 prototype strain, it had 79.6% and 76.3% nucleotide identity in the complete genome and VP1 coding region, respectively. VP1 sequence analyses revealed also high variation in nucleotide similarity (68.9%-82.8%) with previously isolated EV-B69 strains in India and Africa. The great genetic divergence among EV-B69 strains indicates that this type is not a newly emergent virus, but has circulated for many years at low epidemic strength. Phylogenetic incongruity between structural and non-structural regions and similarity plot analyses revealed that multiple recombination events occurred during its evolution. This study expands the number of EV-B69 whole genome sequences which would help genomic comparison for future studies to understand the biological and pathogenic properties of this virus, assess its potential public health impact and comprehend the role of recombination in the evolution of enteroviruses.

Virome Sequencing of Stool Samples

Next-generation sequencing has opened avenues to studying complex populations such as the bacteriome (all bacteria), mycobiome (all fungi), and virome (all viruses in a given sample). Viromes are less often investigated as compared to bacteriomes. The reasons are mostly methodological: because no common pan-viral sequence signature exists, metagenomic sequencing remains the only option. This brings about the need of laborious virus enrichment, multiple signal amplification steps with virtually no possibility of interim quality control, and complicated bioinformatic analysis of the ensuing sequence data. Nevertheless, over the past decade virome sequencing has been enormously successful in identifying new agents in human and animal diseases, and in characterizing viruses in various ecological niches. Recently, virome sequencing has been also employed in studies of non-infectious diseases, which has brought about new challenges of sensitivity, costs, and reproducibility in testing of large sets of samples. Here, we present a detailed protocol that has been utilized in virome studies where hundreds of samples had to be reliably tested in order to assess the association of the stool virome with susceptibility to type 1 diabetes, a non-infectious autoimmune disease.

The Complete Genome Sequence of a Human Parechovirus from a Child with Diarrhea in China Revealed Intertypic Recombination

A human parechovirus (HPeV), CH-ZXY1, was detected in feces from a child with diarrhea. Phylogenetic trees over three different genomic regions revealed discordant topological structures. Recombination analysis indicates that CH-ZXY1 is a recombinant resulting from recombination between HPeV5 and HPeV1, which was confirmed by PCR covering the recombination breakpoint.

Identification of Novel Rosavirus Species That Infects Diverse Rodent Species and Causes Multisystemic Dissemination in Mouse Model

While novel picornaviruses are being discovered in rodents, their host range and pathogenicity are largely unknown. We identified two novel picornaviruses, rosavirus B from the street rat, Norway rat, and rosavirus C from five different wild rat species (chestnut spiny rat, greater bandicoot rat, Indochinese forest rat, roof rat and Coxing's white-bellied rat) in China. Analysis of 13 complete genome sequences showed that “Rosavirus B” and “Rosavirus C” represent two potentially novel picornavirus species infecting different rodents. Though being most closely related to rosavirus A, rosavirus B and C possessed distinct protease cleavage sites and variations in Yn-Xm-AUG sequence in 5’UTR and myristylation site in VP4. Anti-rosavirus B VP1 antibodies were detected in Norway rats, whereas anti-rosavirus C VP1 and neutralizing antibodies were detected in Indochinese forest rats and Coxing's white-bellied rats. While the highest prevalence was observed in Coxing's white-bellied rats by RT-PCR, the detection of rosavirus C from different rat species suggests potential interspecies transmission. Rosavirus C isolated from 3T3 cells causes multisystemic diseases in a mouse model, with high viral loads and positive viral antigen expression in organs of infected mice after oral or intracerebral inoculation. Histological examination revealed alveolar fluid exudation, interstitial infiltration, alveolar fluid exudate and wall thickening in lungs, and hepatocyte degeneration and lymphocytic/monocytic inflammatory infiltrates with giant cell formation in liver sections of sacrificed mice. Since rosavirus A2 has been detected in fecal samples of children, further studies should elucidate the pathogenicity and emergence potential of different rosaviruses.

We identified two novel picornaviruses, rosavirus B and C, infecting street and wild rats respectively in China. While rosavirus B was detected from Norway rats, rosavirus C was detected from five different wild rat species (chestnut spiny rat, greater bandicoot rat, Indochinese forest rat, roof rat and Coxing's white-bellied rat) by RT-PCR. Anti-rosavirus B antibodies were detected in Norway rats, whereas anti-rosavirus C antibodies were detected in Indochinese forest rats and Coxing's white-bellied rats, supporting potential interspecies transmission of rosavirus C. Genome analysis supported the classification of rosavirus B and C as two novel picornavirus species, with genome features distinct from rosavirus A. Rosavirus C isolated from 3T3 cells causes multisystemic diseases in a mouse model, with viruses and pathologies detected in various organs of infected mice after oral or intracerebral inoculation. Our results extend our knowledge on the host range and pathogenicity of rodent picornaviruses.

[Natural reservoirs of viruses of the genus Hepacivirus, Flaviviridae]

HCV is a cause of acute and chronic liver diseases, including chronic hepatitis, cirrhosis, and hepatocellular carcinoma. Under natural conditions, HCV is able to infect only humans, and only chimpanzees are sensitive to experimental infection. In recent years, viruses genetically related to HCV were discovered in wild mammals (rodents, bats, rabbits), as well as in domestic animals living in close contact with humans (dogs, horses, cows). The hepacivirus genus of the family Flaviviridae, previously represented only by HCV and, presumably, by GBV-B, now includes new related viruses of animals. The results of the study of molecular-genetic and biological properties of the hepaciviruses provide an opportunity to understand the history, evolution, and the origin of HCV. It also opens up the prospect of using HCV homologues of non-primates as a laboratory model for preclinical medical and prophylactic drugs against hepatitis c. It was found that the hepacivirus of horses is the most closely related to HCV among currently known HCV homologues.

Molecular epidemiology of canine picornavirus in Hong Kong and Dubai and proposal of a novel genus in Picornaviridae

Previously, we reported the discovery of a novel canine picornavirus (CanPV) in the fecal sample of a dog. In this molecular epidemiology study, CanPV was detected in 15 (1.11%) of 1347 canine fecal samples from Hong Kong and one (0.76%) of 131 canine fecal samples from Dubai, with viral loads 1.06×10(3) to 6.64×10(6) copies/ml. Complete genome sequencing and phylogenetic analysis showed that CanPV was clustered with feline picornavirus (FePV), bat picornavirus (BatPV) 1 to 3, Ia io picornavirus 1 (IaioPV1) and bovine picornavirus (BoPV), and this cluster was most closely related to the genera Enterovirus and Sapelovirus. The Ka/Ks ratios of all the coding regions were <0.1. According to the definition of the Picornavirus Study Group of ICTV, CanPV, FePV, BatPV 1 to 3, IaioPV1 and BoPV should constitute a novel genus in Picornaviridae. BEAST analysis showed that this genus diverged from its most closely related genus, Sapelovirus, about 49 years ago.

Metagenomic Survey of Viral Diversity Obtained from Feces of Subantarctic and South American Fur Seals

The Brazilian South coast seasonally hosts numerous marine species, observed particularly during winter months. Some animals, including fur seals, are found dead or debilitated along the shore and may harbor potential pathogens within their microbiota. In the present study, a metagenomic approach was performed to evaluate the viral diversity in feces of fur seals found deceased along the coast of the state of Rio Grande do Sul. The fecal virome of two fur seal species was characterized: the South American fur seal (Arctocephalus australis) and the Subantarctic fur seal (Arctocephalus tropicalis). Fecal samples from 10 specimens (A. australis, n = 5; A. tropicalis, n = 5) were collected and viral particles were purified, extracted and amplified with a random PCR. The products were sequenced through Ion Torrent and Illumina platforms and assembled reads were submitted to BLASTx searches. Both viromes were dominated by bacteriophages and included a number of potentially novel virus genomes. Sequences of picobirnaviruses, picornaviruses and a hepevirus-like were identified in A. australis. A rotavirus related to group C, a novel member of the Sakobuvirus and a sapovirus very similar to California sea lion sapovirus 1 were found in A. tropicalis. Additionally, sequences of members of the Anelloviridae and Parvoviridae families were detected in both fur seal species. This is the first metagenomic study to screen the fecal virome of fur seals, contributing to a better understanding of the complexity of the viral community present in the intestinal microbiota of these animals.

A novel gemycircularvirus in an unexplained case of child encephalitis

Background

Recently, a diverse group of viruses with circular, replication initiator protein(Rep) encoding, single stranded DNA (CRESS-DNA) genomes, were discovered from wide range of eukaryotic organisms ranging from mammals to fungi. Gemycircularvirus belongs to a distinct group of CRESS-DNA genomes and is classified under the genus name of Gemycircularvirus.

Findings

Here, a novel gemycircularvirus named GeTz1 from cerebrospinal fluid sample of a child with unexplainable encephalitis was characterized. The novel gemycircularvirus encodes two major proteins, including a capsid protein (Cap) and a replication-associated protein (Rep). Phylogenetic analysis based on the amino acid sequence of Rep indicated that GeTz1 clusters with one gemycircularvirus discovered from bird (KF371633), sharing 46.6 % amino acid sequence identity with each other.

Conclusion

A novel gemycircularvirus was discovered from cerebrospinal fluid sample of a child with unexplainable encephalitis. Further studies, such as testing human sera for specific antibodies, should be performed to investigate whether gemycircularvirus infects human and is associated with encephalitis.

HERPESVIRUSES INCLUDING NOVEL GAMMAHERPESVIRUSES ARE WIDESPREAD AMONG PHOCID SEAL SPECIES IN CANADA

Little is known about herpesviruses in Canadian pinnipeds. We measured prevalence of antibodies to herpesviruses in the sera from Canadian phocid seals by an indirect enzyme-linked immunosorbent assay. Wild harbor seals (Phoca vitulina) and captive harbor seals were positive for antibodies to Phocid herpesvirus 1 (PhoHV-1) at prevalences of 91% and 100%, respectively. Sera from wild hooded seals (Cystophora cristata), harp seals (Pagophilus groenlandica), and grey seals (Halichoerus grypus) were positive for antibodies to PhoHV-1 antigenically related herpesvirus antigens at 73%, 79%, and 96%, respectively. We isolated new herpesviruses in cell culture from two hunter-harvested ringed seals (Pusa hispida) in poor body condition from Ulukhaktok, Northwest Territories, Canada; one lethargic hooded seal from the St. Lawrence Estuary, Québec, Canada; and one captive, asymptomatic harp seal from the Magdalen Islands, Québec. Partial sequencing of the herpesvirus DNA polymerase gene revealed that all four virus isolates were closely related to PhoHV-2, a member of the Gammaherpesvirinae subfamily, with nucleotide similarity ranging between 92.8% and 95.3%. The new seal herpesviruses were genetically related to other known pinniped herpesviruses, such as PhoHV-1, Otariid herpesvirus 3, Hawaiian monk (Monachus schauinslandi) seal herpesvirus, and Phocid herpesvirus 5 with 47-48%, 55%, 77%, and 70-77% nucleotide similarities, respectively. The harp seal herpesvirus and both ringed seal herpesviruses were almost identical to each other, whereas the hooded seal herpesvirus was genetically different from the three others (92.8% nucleotide similarity), indicating detection of at least two novel seal herpesviruses. These findings are the first isolation, partial genome sequencing, and identification of seal gammaherpesviruses in three species of Canadian phocid seals; two species of which were suspected of exposure to one or more antigenically related herpesviruses based on serologic analyses.

Virome Capture Sequencing Enables Sensitive Viral Diagnosis and Comprehensive Virome Analysis

Insensitivity and technical complexity have impeded the implementation of high-throughput nucleic acid sequencing in differential diagnosis of viral infections in clinical laboratories. Here, we describe the development of a virome capture sequencing platform for vertebrate viruses (VirCapSeq-VERT) that increases the sensitivity of sequence-based virus detection and characterization. The system uses ~2 million probes that cover the genomes of members of the 207 viral taxa known to infect vertebrates, including humans. A biotinylated oligonucleotide library was synthesized on the NimbleGen cleavable array platform and used for solution-based capture of viral nucleic acids present in complex samples containing variable proportions of viral and host nucleic acids. The use of VirCapSeq-VERT resulted in a 100- to 10,000-fold increase in viral reads from blood and tissue homogenates compared to conventional Illumina sequencing using established virus enrichment procedures, including filtration, nuclease treatments, and RiboZero rRNA subtraction. VirCapSeq-VERT had a limit of detection comparable to that of agent-specific real-time PCR in serum, blood, and tissue extracts. Furthermore, the method identified novel viruses whose genomes were approximately 40% different from the known virus genomes used for designing the probe library. The VirCapSeq-VERT platform is ideally suited for analyses of virome composition and dynamics.

Importance  VirCapSeq-VERT enables detection of viral sequences in complex sample backgrounds, including those found in clinical specimens, such as serum, blood, and tissue. The highly multiplexed nature of the system allows both the simultaneous identification and the comprehensive genetic characterization of all known vertebrate viruses, their genetic variants, and novel viruses. The operational simplicity and efficiency of the VirCapSeq-VERT platform may facilitate transition of high-throughput sequencing to clinical diagnostic as well as research applications.

VirCapSeq-VERT enables detection of viral sequences in complex sample backgrounds, including those found in clinical specimens, such as serum, blood, and tissue. The highly multiplexed nature of the system allows both the simultaneous identification and the comprehensive genetic characterization of all known vertebrate viruses, their genetic variants, and novel viruses. The operational simplicity and efficiency of the VirCapSeq-VERT platform may facilitate transition of high-throughput sequencing to clinical diagnostic as well as research applications.

A highly divergent picornavirus in an amphibian, the smooth newt (Lissotriton vulgaris)

Genetically highly divergent picornavirus (Newt/2013/HUN, KP770140) was detected using viral metagenomics in faecal samples of free-living smooth newts (Lissotriton vulgaris). Newt picornavirus was identified by reverse transcription-polymerase chain reaction (RT-PCR) in six (25 %) of the 24 samples originating from individuals caught in two out of the six investigated natural ponds in Hungary. The first picornavirus in amphibians expands the host range of members of the Picornaviridae, and opens a new research field in picornavirus evolution in lower vertebrates. Newt picornavirus represents a novel species in a novel genus within the family Picornaviridae, provisionally named genus Ampivirus (amphibian picornavirus).

A novel dromedary camel enterovirus in the family Picornaviridae from dromedaries in the Middle East

The recent emergence of Middle East respiratory syndrome coronavirus from the Middle East and the discovery of the virus from dromedary camels have boosted interest in the search for novel viruses in dromedaries. Whilst picornaviruses are known to infect various animals, their existence in dromedaries was unknown. We describe the discovery of a novel picornavirus, dromedary camel enterovirus (DcEV), from dromedaries in Dubai. Among 215 dromedaries, DcEV was detected in faecal samples of four (1.9 %) dromedaries [one (0.5 %) adult dromedary and three (25 %) dromedary calves] by reverse transcription PCR. Analysis of two DcEV genomes showed that DcEV was clustered with other species of the genus Enterovirus and was most closely related to and possessed highest amino acid identities to the species Enterovirus E and Enterovirus F found in cattle. The G+C content of DcEV was 45 mol%, which differed from that of Enterovirus E and Enterovirus F (49-50 mol%) by 4-5 %. Similar to other members of the genus Enterovirus, the 5' UTR of DcEV possessed a putative type I internal ribosome entry site. The low ratios of the number of nonsynonymous substitutions per non-synonymous site to the number of synonymous substitutions per synonymous site (Ka/Ks) of various coding regions suggested that dromedaries are the natural reservoir in which DcEV has been stably evolving. These results suggest that DcEV is a novel species of the genus Enterovirus in the family Picornaviridae. Western blot analysis using recombinant DcEV VP1 polypeptide showed a high seroprevalence of 52 % among serum samples from 172 dromedaries for IgG, concurring with its much higher infection rates in dromedary calves than in adults. Further studies are important to understand the pathogenicity, epidemiology and genetic evolution of DcEV in this unique group of animals.

Widespread distribution and structural diversity of Type IV IRESs in members of Picornaviridae

Picornavirus genomes contain internal ribosomal entry sites (IRESs) that promote end-independent translation initiation. Five structural classes of picornavirus IRES have been identified, but numerous IRESs remain unclassified. Here, previously unrecognized Type IV IRESs were identified in members of three proposed picornavirus genera (Limnipivirus, Pasivirus, Rafivirus) and four recognized genera (Kobuvirus, Megrivirus, Sapelovirus, Parechovirus). These IRESs are ~230-420 nucleotides long, reflecting heterogeneity outside a common structural core. Closer analysis yielded insights into evolutionary processes that have shaped contemporary IRESs. The presence of related IRESs in diverse genera supports the hypothesis that they are heritable genetic elements that spread by horizontal gene transfer. Recombination likely also accounts for the exchange of some peripheral subdomains, suggesting that IRES evolution involves incremental addition of elements to a pre-existing core. Nucleotide conservation is concentrated in ribosome-binding sites, and at the junction of helical domains, likely to ensure orientation of subdomains in an active conformation.

Picornavirus RNA polyadenylation by 3D(pol), the viral RNA-dependent RNA polymerase

•

Picornaviral RdRPs are responsible for the polyadenylation of viral RNA.

•

Reiterative transcription mechanisms occur during replication of poly(A) tails.

•

Conserved RdRP structures influence the size of poly(A) tails.

•

Common features of picornavirus RdRPs and telomerase reverse transcriptase.

•

Poly(A) tails are a telomere of picornavirus RNA genomes.

Poly(A) tails are functionally important features of all picornavirus RNA genomes. Some viruses have genomes with relatively short poly(A) tails (encephalomyocarditis virus) whereas others have genomes with longer poly(A) tails (polioviruses and rhinoviruses). Here we review the polyadenylation of picornavirus RNA as it relates to the structure and function of 3D^pol. Poliovirus 3D^pol uses template-dependent reiterative transcription mechanisms as it replicates the poly(A) tails of viral RNA (Steil et al., 2010). These mechanisms are analogous to those involved in the polyadenylation of vesicular stomatitis virus and influenza virus mRNAs. 3D^pol residues intimately associated with viral RNA templates and products regulate the size of poly(A) tails in viral RNA (Kempf et al., 2013). Consistent with their ancient evolutionary origins, picornavirus 3D^pol and telomerase reverse transcriptase (TERT) share structural and functional features. Structurally, both 3D^pol and TERT assume a “right-hand” conformation with thumb, palm and fingers domains encircling templates and products. Functionally, both 3D^pol and TERT use template-dependent reiterative transcription mechanisms to synthesize repetitive sequences: poly(A) tails in the case of picornavirus RNA genomes and DNA telomeres in the case of eukaryotic chromosomes. Thus, picornaviruses and their eukaryotic hosts (humans and animals) maintain the 3′ ends of their respective genomes via evolutionarily related mechanisms.

Metagenomic analysis of the shrew enteric virome reveals novel viruses related to human stool-associated viruses

Shrews are small insectivorous mammals that are distributed worldwide. Similar to rodents, shrews live on the ground and are commonly found near human residences. In this study, we investigated the enteric virome of wild shrews in the genus Crocidura using a sequence-independent viral metagenomics approach. A large portion of the shrew enteric virome was composed of insect viruses, whilst novel viruses including cyclovirus, picornavirus and picorna-like virus were also identified. Several cycloviruses, including variants of human cycloviruses detected in cerebrospinal fluid and stools, were detected in wild shrews at a high prevalence rate. The identified picornavirus was distantly related to human parechovirus, inferring the presence of a new genus in this family. The identified picorna-like viruses were characterized as different species of calhevirus 1, which was discovered previously in human stools. Complete or nearly complete genome sequences of these novel viruses were determined in this study and then were subjected to further genetic characterization. Our study provides an initial view of the diversity and distinctiveness of the shrew enteric virome and highlights unique novel viruses related to human stool-associated viruses.

DNA microarray for detection of gastrointestinal viruses

Gastroenteritis is a clinical illness of humans and other animals that is characterized by vomiting and diarrhea and caused by a variety of pathogens, including viruses. An increasing number of viral species have been associated with gastroenteritis or have been found in stool samples as new molecular tools have been developed. In this work, a DNA microarray capable in theory of parallel detection of more than 100 viral species was developed and tested. Initial validation was done with 10 different virus species, and an additional 5 species were validated using clinical samples. Detection limits of 1 × 10(3) virus particles of Human adenovirus C (HAdV), Human astrovirus (HAstV), and group A Rotavirus (RV-A) were established. Furthermore, when exogenous RNA was added, the limit for RV-A detection decreased by one log. In a small group of clinical samples from children with gastroenteritis (n = 76), the microarray detected at least one viral species in 92% of the samples. Single infection was identified in 63 samples (83%), and coinfection with more than one virus was identified in 7 samples (9%). The most abundant virus species were RV-A (58%), followed by Anellovirus (15.8%), HAstV (6.6%), HAdV (5.3%), Norwalk virus (6.6%), Human enterovirus (HEV) (9.2%), Human parechovirus (1.3%), Sapporo virus (1.3%), and Human bocavirus (1.3%). To further test the specificity and sensitivity of the microarray, the results were verified by reverse transcription-PCR (RT-PCR) detection of 5 gastrointestinal viruses. The RT-PCR assay detected a virus in 59 samples (78%). The microarray showed good performance for detection of RV-A, HAstV, and calicivirus, while the sensitivity for HAdV and HEV was low. Furthermore, some discrepancies in detection of mixed infections were observed and were addressed by reverse transcription-quantitative PCR (RT-qPCR) of the viruses involved. It was observed that differences in the amount of genetic material favored the detection of the most abundant virus. The microarray described in this work should help in understanding the etiology of gastroenteritis in humans and animals.

Epidemiology of human parechovirus, Aichi virus and salivirus in fecal samples from hospitalized children with gastroenteritis in Hong Kong

BACKGROUND

Emerging human picornaviruses, including human parechovirus (HPeV), Aichi virus (AiV) and salivirus (SalV) were found to be associated with gastroenteritis, but their roles in enteric infections are not fully understood. In addition, no report on the circulation of these viruses in Hong Kong is available. The objective of this study was to investigate the prevalence and genetic diversity of HPeV, AiV and SalV in fecal samples from hospitalized children with gastroenteritis in Hong Kong.

METHODS

Fecal samples from hospitalized children with gastroenteritis were subject to detection of HPeV, AiV and SalV by RT-PCR using consensus primers targeted to their 5'UTRs. Positive samples were subject to capsid and/or 3CD region analysis for genotype determination. The epidemiology of HPeV, AiV and SalV infections was analyzed.

RESULTS

Among 1,708 fecal samples subjected to RT-PCR using primers targeted to 5'UTR of HPeV, AiV and SalV, viruses were detected in 55 samples, with 50 positive for HPeV only, 3 positive for AiV only, 1 positive for both HPeV and AiV, and 1 positive for both HPeV and SalV. Phylogenetic analysis of the partial VP1 gene of the 33 HPeV strains revealed the presence of genotypes of HPeV- 1, 3, 4, 5, 7, 10, among which HPeV-1 was the predominant genotype circulating in our population. The peak activity of HPeV infection was in fall. Of the 3 children with AiV infection, the 3 AiV strains were found to belong to genotype A based on the phylogenetic analysis of their partial VP1 and 3CD regions. The genotype of a SalV strain detected in this study could not be determined. Co-detection of different pathogens was observed in 24 samples (43.6%) of 55 fecal samples positive for HPeV, AiV and SalV.

CONCLUSIONS

HPeV, AiV and SalV were detected in fecal samples of hospitalized children with gastroenteritis in Hong Kong, with the former having the highest prevalence. HPeV-1 was the predominant genotype among HPeVs, while genotype A was the predominant genotype among AiVs in this study.

Avian picornaviruses: molecular evolution, genome diversity and unusual genome features of a rapidly expanding group of viruses in birds

Picornaviridae is one of the most diverse families of viruses infecting vertebrate species. In contrast to the relative small number of mammal species compared to other vertebrates, the abundance of mammal-infecting picornaviruses was significantly overrepresented among the presently known picornaviruses. Therefore most of the current knowledge about the genome diversity/organization patterns and common genome features were based on the analysis of mammal-infecting picornaviruses. Beside the well known reservoir role of birds in case of several emerging viral pathogens, little is known about the diversity of picornaviruses circulating among birds, although in the last decade the number of known avian picornavirus species with complete genome was increased from one to at least 15. However, little is known about the geographic distribution, host spectrum or pathogenic potential of the recently described picornaviruses of birds. Despite the low number of known avian picornaviruses, the phylogenetic and genome organization diversity of these viruses were remarkable. Beside the common L-4-3-4 and 4-3-4 genome layouts unusual genome patterns (3-4-4; 3-5-4, 3-6-4; 3-8-4) with variable, multicistronic 2A genome regions were found among avian picornaviruses. The phylogenetic and genomic analysis revealed the presence of several conserved structures at the untranslated regions among phylogenetically distant avian and non-avian picornaviruses as well as at least five different avian picornavirus phylogenetic clusters located in every main picornavirus lineage with characteristic genome layouts which suggests the complex evolution history of these viruses. Based on the remarkable genetic diversity of the few known avian picornaviruses, the emergence of further divergent picornaviruses causing challenges in the current taxonomy and also in the understanding of the evolution and genome organization of picornaviruses will be strongly expected. In this review we would like to summarize the current knowledge about the taxonomy, pathogenic potential, phylogenetic/genomic diversity and evolutional relationship of avian picornaviruses.

Chickens host diverse picornaviruses originated from potential interspecies transmission with recombination

While chickens are an important reservoir for emerging pathogens such as avian influenza viruses, little is known about the diversity of picornaviruses in poultry. We discovered a previously unknown diversity of picornaviruses in chickens in Hong Kong. Picornaviruses were detected in 87 cloacal and 7 tracheal samples from 93 of 900 chickens by reverse transcription-PCR, with their partial 3D(pol) gene sequences forming five distinct clades (I to V) among known picornaviruses. Analysis of eight genomes from different clades revealed seven different picornaviruses, including six novel picornavirus species (ChPV1 from clade I, ChPV2 and ChPV3 from clade II, ChPV4 and ChPV5 from clade III, ChGV1 from clade IV) and one existing species (Avian encephalomyelitis virus from clade V). The six novel chicken picornavirus genomes exhibited distinct phylogenetic positions and genome features different from related picornaviruses, supporting their classification as separate species. Moreover, ChPV1 may potentially belong to a novel genus, with low sequence homologies to related picornaviruses, especially in the P1 and P2 regions, including the predicted L and 2A proteins. Nevertheless, these novel picornaviruses were most closely related to picornaviruses of other avian species (ChPV1 related to Passerivirus A, ChPV2 and ChPV3 to Avisivirus A and Duck hepatitis A virus, ChPV4 and ChPV5 to Melegrivirus A, ChGV1 to Gallivirus A). Furthermore, ChPV5 represented a potential recombinant picornavirus, with its P2 and P3 regions possibly originating from Melegrivirus A. Chickens are an important reservoir for diverse picornaviruses that may cross avian species barriers through mutation or recombination.

Complete genome sequence of a novel calicivirus from a goose

A novel goose calicivirus (GoCV) was sequenced. The 8013-nt-long genome was organized into two open reading frames that were in the same frame and separated by 3 nucleotides. This feature is similar to what has been observed in turkey calicivirus (TuCV). Comparison of GoCV with other caliciviruses showed that it shared the highest amino acid sequence identities of 62, 38, and 52% in the nonstructural protein, VP1, and VP2, respectively, with TuCV. Phylogenetic analysis based on the amino acid sequences of nonstructural protein and VP1 demonstrated that GoCV was most closely related to but distinct from TuCV. Thus, GoCV was identified as a novel member in the proposed genus Nacovirus.

Genetic characterization of a novel duck-origin picornavirus with six 2A proteins

A novel virus was detected from diseased ducks and completely determined. The virus was shown to have a picornavirus-like genome layout. Interestingly, the genome contained a total of up to six 2As, including four 2As (2A1-2A4) each having an NPGP motif, an AIG1-like 2A5, and a parechovirus-like 2A6. The 5'UTR was predicted to possess a hepacivirus/pestivirus-like internal ribosome entry site (IRES). However, the subdomain IIIe consisted of a 3 nt stem and five unpaired bases, distinct from those found in all other HP-like IRESs. The virus was most closely related to duck hepatitis A virus, with amino acid identities of 37.7 %, 39 % and 43.7 % in the P1, P2 and P3 regions, respectively. Based on these investigations, together with phylogenetic analyses, the virus could be considered as the founding member of a novel picornavirus genus that we tentatively named 'Aalivirus', with 'Aalivirus A' as the type species.

Isolation and molecular characterization of a novel picornavirus from baitfish in the USA

During both regulatory and routine surveillance sampling of baitfish from the states of Illinois, Minnesota, Montana, and Wisconsin, USA, isolates (n = 20) of a previously unknown picornavirus were obtained from kidney/spleen or entire viscera of fathead minnows (Pimephales promelas) and brassy minnows (Hybognathus hankinsoni). Following the appearance of a diffuse cytopathic effect, examination of cell culture supernatant by negative contrast electron microscopy revealed the presence of small, round virus particles (∼ 30-32 nm), with picornavirus-like morphology. Amplification and sequence analysis of viral RNA identified the agent as a novel member of the Picornaviridae family, tentatively named fathead minnow picornavirus (FHMPV). The full FHMPV genome consisted of 7834 nucleotides. Phylogenetic analysis based on 491 amino acid residues of the 3D gene showed 98.6% to 100% identity among the 20 isolates of FHMPV compared in this study while only 49.5% identity with its nearest neighbor, the bluegill picornavirus (BGPV) isolated from bluegill (Lepomis macrochirus). Based on complete polyprotein analysis, the FHMPV shared 58% (P1), 33% (P2) and 43% (P3) amino acid identities with BGPV and shared less than 40% amino acid identity with all other picornaviruses. Hence, we propose the creation of a new genus (Piscevirus) within the Picornaviridae family. The impact of FHMPV on the health of fish populations is unknown at present.

Characterization of a new picornavirus isolated from the freshwater fish Lepomis macrochirus

The freshwater fish Lepomis macrochirus (bluegill) is common to North American waters, and important both ecologically and as a sport fish. In 2001 an unknown virus was isolated from bluegills following a bluegill fish kill. This virus was identified as a picornavirus [termed bluegill picornavirus (BGPV)] and a diagnostic reverse transcriptase PCR was developed. A survey of bluegills in Wisconsin waters showed the presence of BGPV in 5 of 17 waters sampled, suggesting the virus is widespread in bluegill populations. Experimental infections of bluegills confirmed that BGPV can cause morbidity and mortality in bluegills. Molecular characterization of BGPV revealed several distinct genome characteristics, the most unusual of which is the presence of a short poly(C) tract in the 3' UTR. Additionally, the genome encodes a polyprotein lacking a leader peptide and a VP0 maturation cleavage site, and is predicted to encode two distinct 2A proteins. Sequence comparison showed that the virus is most closely related to a phylogenetic cluster of picornaviruses that includes the genera Aquamavirus, Avihepatovirus and Parechovirus. However, it is distinct enough, for example sharing only about 38% sequence identity to the parechoviruses in the 3D region, that it may represent a new genus in the family Picornaviridae.

Novel B19-like parvovirus in the brain of a harbor seal

Using random PCR in combination with next-generation sequencing, a novel parvovirus was detected in the brain of a young harbor seal (Phoca vitulina) with chronic non-suppurative meningo-encephalitis that was rehabilitated at the Seal Rehabilitation and Research Centre (SRRC) in the Netherlands. In addition, two novel viruses belonging to the family Anelloviridae were detected in the lungs of this animal. Phylogenetic analysis of the coding sequence of the novel parvovirus, tentatively called Seal parvovirus, indicated that this virus belonged to the genus Erythrovirus, to which human parvovirus B19 also belongs. Although no other seals with similar signs were rehabilitated in SRRC in recent years, a prevalence study of tissues of seals from the same area collected in the period 2008-2012 indicated that the Seal parvovirus has circulated in the harbor seal population at least since 2008. The presence of the Seal parvovirus in the brain was confirmed by real-time PCR and in vitro replication. Using in situ hybridization, we showed for the first time that a parvovirus of the genus Erythrovirus was present in the Virchow-Robin space and in cerebral parenchyma adjacent to the meninges. These findings showed that a parvovirus of the genus Erythrovirus can be involved in central nervous system infection and inflammation, as has also been suspected but not proven for human parvovirus B19 infection.

Characterization of a novel picornavirus isolate from a diseased European eel (Anguilla anguilla)

A novel picornavirus was isolated from specimens of a diseased European eel (Anguilla anguilla). This virus induced a cytopathic effect in eel embryonic kidney cells and high mortality in a controlled transmission study using elvers. Eel picornavirus has a genome of 7,496 nucleotides that encodes a polyprotein of 2,259 amino acids. It has a typical picornavirus genome layout, but its low similarity to known viral proteins suggests a novel species in the family Picornaviridae.

The origin of biased sequence depth in sequence-independent nucleic acid amplification and optimization for efficient massive parallel sequencing

Sequence Independent Single Primer Amplification is one of the most widely used random amplification approaches in virology for sequencing template preparation. This technique relies on oligonucleotides consisting of a 3' random part used to prime complementary DNA synthesis and a 5' defined tag sequence for subsequent amplification. Recently, this amplification method was combined with next generation sequencing to obtain viral sequences. However, these studies showed a biased distribution of the resulting sequence reads over the analyzed genomes. The aim of this study was to elucidate the mechanisms that lead to biased sequence depth when using random amplification. Avian paramyxovirus type 8 was used as a model RNA virus to investigate these mechanisms. We showed, based on in silico analysis of the sequence depth in relation to GC-content, predicted RNA secondary structure and sequence complementarity to the 3' part of the tag sequence, that the tag sequence has the main contribution to the observed bias in sequence depth. We confirmed this finding experimentally using both fragmented and non-fragmented viral RNAs as well as primers differing in random oligomer length (6 or 12 nucleotides) and in the sequence of the amplification tag. The observed oligonucleotide annealing bias can be reduced by extending the random oligomer sequence and by in silico combining sequence data from SISPA experiments using different 5' defined tag sequences. These findings contribute to the optimization of random nucleic acid amplification protocols that are currently required for downstream applications such as viral metagenomics and microarray analysis.

Evaluation of rapid and simple techniques for the enrichment of viruses prior to metagenomic virus discovery

•

The effect of simple virus enrichment methods were tested on a metagenomics dataset.

•

Centrifugation, filtration or nuclease-treatment was evaluated.

•

A multi-step enrichment method increased the proportion of virus sequences.

•

This evaluation guides researchers in their choice of enrichment methodology.

The discovery of new or divergent viruses using metagenomics and high-throughput sequencing has become more commonplace. The preparation of a sample is known to have an effect on the representation of virus sequences within the metagenomic dataset yet comparatively little attention has been given to this. Physical enrichment techniques are often applied to samples to increase the number of viral sequences and therefore enhance the probability of detection. With the exception of virus ecology studies, there is a paucity of information available to researchers on the type of sample preparation required for a viral metagenomic study that seeks to identify an aetiological virus in an animal or human diagnostic sample. A review of published virus discovery studies revealed the most commonly used enrichment methods, that were usually quick and simple to implement, namely low-speed centrifugation, filtration, nuclease-treatment (or combinations of these) which have been routinely used but often without justification. These were applied to a simple and well-characterised artificial sample composed of bacterial and human cells, as well as DNA (adenovirus) and RNA viruses (influenza A and human enterovirus), being either non-enveloped capsid or enveloped viruses. The effect of the enrichment method was assessed by both quantitative real-time PCR and metagenomic analysis that incorporated an amplification step. Reductions in the absolute quantities of bacteria and human cells were observed for each method as determined by qPCR, but the relative abundance of viral sequences in the metagenomic dataset remained largely unchanged. A 3-step method of centrifugation, filtration and nuclease-treatment showed the greatest increase in the proportion of viral sequences. This study provides a starting point for the selection of a purification method in future virus discovery studies, and highlights the need for more data to validate the effect of enrichment methods on different sample types, amplification, bioinformatics approaches and sequencing platforms. This study also highlights the potential risks that may attend selection of a virus enrichment method without any consideration for the sample type being investigated.