Novel Picornavirus in Lambs with Severe Encephalomyelitis

Identification and Genomic Characterization of Bovine Boosepivirus A in the United States and Canada

Boosepivirus is a new genus in the Picornaviridae family. Boosepiviruses (BooVs) are genetically classified into three species: A, B, and C. Initially, Boosepivirus A and B were identified in cattle, whereas Boosepivirus C was detected in sheep. Recent evidence showed that Boosepivirus B was detected in sheep and Boosepivirus C was identified in goats, suggesting that Boosepvirus might cross the species barrier to infect different hosts. Different from BooV B, BooV A is less studied. In the present study, we reported identification of two North American BooV A strains from cattle. Genomic characterization revealed that US IL33712 (GenBank accession #PP035161) and Canada 1087562 (GenBank accession #PP035162) BooV A strains are distantly related to each other, and US IL33712 is more closely correlated to two Asian BooV A strains. US-strain-specific insertions, NorthAmerican-strain-specific insertions, and species A-specific insertions are observed and could contribute to viral pathogenicity and host adaptation. Our findings highlight the importance of continued surveillance of BooV A in animals.

Identification and Genomic Characterization of Bovine Boosepivirus A in the United States and Canada

Boosepivirus is a new genus in the Picornaviridae family. Boosepiviruses (BooVs) are genetically classified into three species: A, B, and C. Initially, Boosepivirus A and B were identified in cattle, whereas Boosepivirus C was detected in sheep. Recent evidence showed that Boosepivirus B was detected in sheep and Boosepivirus C was identified in goats, suggesting that Boosepvirus might cross the species barrier to infect different hosts. Different from BooV B, BooV A is less studied. In the present study, we reported identification of two North American BooV A strains from cattle. Genomic characterization revealed that US IL33712 (GenBank accession #PP035161) and Canada 1087562 (GenBank accession #PP035162) BooV A strains are distantly related to each other, and US IL33712 is more closely correlated to two Asian BooV A strains. US-strain-specific insertions, NorthAmerican-strain-specific insertions, and species A-specific insertions are observed and could contribute to viral pathogenicity and host adaptation. Our findings highlight the importance of continued surveillance of BooV A in animals.

The genomic and epidemiological investigations of enteric viruses of domestic caprine (Capra hircus) revealed the presence of multiple novel viruses related to known strains of humans and ruminant livestock species

IMPORTANCE

Compared with other domestic animals, the virome and viral diversity of small ruminants especially in caprine are less studied even of its zoonotic potential. In this study, the enteric virome of caprine was investigated in detail using next-generation sequencing and reverse transcription PCR techniques. The complete or nearly complete genomes of seven novel viruses were determined which show a close phylogenetic relationship to known human and ruminant viruses. The high similarity between the identified caprine tusavirus (family Parvoviridae) and an unassigned CRESS DNA virus with closely related human strains could indicate the (reverse) zoonotic potential of these viruses. Others, like astroviruses (family Astroviridae), enteroviruses, or novel caripiviruses (named after the term caprine picornavirus) of family Picornaviridae found mostly in multiple co-infections in caprine and ovine, could indicate the cross-species transmission capabilities of these viruses between small ruminants.

Detection and characterization of bopiviruses in domestic and wild ruminants

Highly divergent picornaviruses (PVs) classified in the genus Bopivirus have been recently discovered on faecal samples from sheep and goats in Hungary and from fallow and red deer in Australia. In this study, we investigated the epidemiology of these novel viruses in domestic and wild ruminants from Northwestern Italian Alps by testing archival faecal samples collected from 128 sheep, 167 goats, 61 red deer (Cervus elaphus), 77 roe deer (Capreolus capreolus), 43 chamois (Rupicapra rupicapra) and 32 Alpine ibex (Capra ibex). Bopivirus RNA was detected in a total of 19 animals, including 14 sheep (10.9%), 2 red deer (3.3%), 1 roe deer (1.3%), 1 chamois (2.3 %) and 1 Alpine ibex (3.3 %), but not in goats. Upon sequence analysis of the 3D^RdRp region, the sequences generated from chamois, roe deer, Alpine ibex and ovine faecal samples showed the highest nucleotide identity (96.8-100%) to bopiviruses detected in goats and sheep from Hungarian farms, whereas strains found in red deer displayed the closest relatedness (90.8%-91.2%) to bopiviruses identified in fallow and red deer in Australia. The nearly complete genome sequence of strains 12/2020/ITA (ON497046) and 14-73/2020/ITA (ON497047) detected in an Alpine ibex and in a sheep, respectively, was determined by combining a modified 3'-RACE protocol with Oxford Nanopore Technologies sequencing platform. On phylogenetic analysis based on the complete polyprotein, both strains segregated into the candidate species Bopivirus B along with ovine and caprine strains detected in Hungary (90.0-94.6% nucleotide and 94.6-98.0% amino acid identities). The findings of this study expand the host range of these novel viruses and hint to a possible virus circulation between domestic ruminants and wild animals.

A novel enterovirus in lambs with poliomyelitis and brain stem encephalitis

An Austrian organic dairy sheep farm experienced cases of recumbency and sudden deaths in 3- to 4-week-old lambs. Two animals were subjected to thorough clinical and pathological investigations. Pathohistological analysis identified severe nonsuppurative myelitis and mild nonsuppurative encephalitis. A reverse-transcription quantitative PCR (RT-qPCR) assay for the recently discovered ovine picornavirus causing comparable lesions scored negative. By next-generation sequencing-based metagenomics, a nearly complete genome of a novel enterovirus could be detected and assembled. In situ hybridization using a specifically designed probe revealed robust signals in affected motoneurons of the spinal cord suggesting a causative role of the novel virus.

In action-an early warning system for the detection of unexpected or novel pathogens

Proactive approaches in preventing future epidemics include pathogen discovery prior to their emergence in human and/or animal populations. Playing an important role in pathogen discovery, high-throughput sequencing (HTS) enables the characterization of microbial and viral genetic diversity within a given sample. In particular, metagenomic HTS allows the unbiased taxonomic profiling of sequences; hence, it can identify novel and highly divergent pathogens such as viruses. Newly discovered viral sequences must be further investigated using genomic characterization, molecular and serological screening, and/or invitro and invivo characterization. Several outbreak and surveillance studies apply unbiased generic HTS to characterize the whole genome sequences of suspected pathogens. In contrast, this study aimed to screen for novel and unexpected pathogens in previously generated HTS datasets and use this information as a starting point for the establishment of an early warning system (EWS). As a proof of concept, the EWS was applied to HTS datasets and archived samples from the 2018–9 West Nile virus (WNV) epidemic in Germany. A metagenomics read classifier detected sequences related to genome sequences of various members of Riboviria. We focused the further EWS investigation on viruses belonging to the families Peribunyaviridae and Reoviridae, under suspicion of causing co-infections in WNV-infected birds. Phylogenetic analyses revealed that the reovirus genome sequences clustered with sequences assigned to the species Umatilla virus (UMAV), whereas a new peribunyavirid, tentatively named ‘Hedwig virus’ (HEDV), belonged to a putative novel genus of the family Peribunyaviridae. In follow-up studies, newly developed molecular diagnostic assays detected fourteen UMAV-positive wild birds from different German cities and eight HEDV-positive captive birds from two zoological gardens. UMAV was successfully cultivated in mosquito C6/36 cells inoculated with a blackbird liver. In conclusion, this study demonstrates the power of the applied EWS for the discovery and characterization of unexpected viruses in repurposed sequence datasets, followed by virus screening and cultivation using archived sample material. The EWS enhances the strategies for pathogen recognition before causing sporadic cases and massive outbreaks and proves to be a reliable tool for modern outbreak preparedness.

Next-generation diagnostics: virus capture facilitates a sensitive viral diagnosis for epizootic and zoonotic pathogens including SARS-CoV-2

BACKGROUND

The detection of pathogens in clinical and environmental samples using high-throughput sequencing (HTS) is often hampered by large amounts of background information, which is especially true for viruses with small genomes. Enormous sequencing depth can be necessary to compile sufficient information for identification of a certain pathogen. Generic HTS combining with in-solution capture enrichment can markedly increase the sensitivity for virus detection in complex diagnostic samples.

METHODS

A virus panel based on the principle of biotinylated RNA baits was developed for specific capture enrichment of epizootic and zoonotic viruses (VirBaits). The VirBaits set was supplemented by a SARS-CoV-2 predesigned bait set for testing recent SARS-CoV-2-positive samples. Libraries generated from complex samples were sequenced via generic HTS (without enrichment) and afterwards enriched with the VirBaits set. For validation, an internal proficiency test for emerging epizootic and zoonotic viruses (African swine fever virus, Ebolavirus, Marburgvirus, Nipah henipavirus, Rift Valley fever virus) was conducted.

RESULTS

The VirBaits set consists of 177,471 RNA baits (80-mer) based on about 18,800 complete viral genomes targeting 35 epizootic and zoonotic viruses. In all tested samples, viruses with both DNA and RNA genomes were clearly enriched ranging from about 10-fold to 10,000-fold for viruses including distantly related viruses with at least 72% overall identity to viruses represented in the bait set. Viruses showing a lower overall identity (38% and 46%) to them were not enriched but could nonetheless be detected based on capturing conserved genome regions. The internal proficiency test supports the improved virus detection using the combination of HTS plus targeted enrichment but also points to the risk of cross-contamination between samples.

CONCLUSIONS

The VirBaits approach showed a high diagnostic performance, also for distantly related viruses. The bait set is modular and expandable according to the favored diagnostics, health sector, or research question. The risk of cross-contamination needs to be taken into consideration. The application of the RNA-baits principle turned out to be user friendly, and even non-experts can easily use the VirBaits workflow. The rapid extension of the established VirBaits set adapted to actual outbreak events is possible as shown for SARS-CoV-2. Video abstract.

Next-generation diagnostics: virus capture facilitates a sensitive viral diagnosis for epizootic and zoonotic pathogens including SARS-CoV-2

BACKGROUND

The detection of pathogens in clinical and environmental samples using high-throughput sequencing (HTS) is often hampered by large amounts of background information, which is especially true for viruses with small genomes. Enormous sequencing depth can be necessary to compile sufficient information for identification of a certain pathogen. Generic HTS combining with in-solution capture enrichment can markedly increase the sensitivity for virus detection in complex diagnostic samples.

METHODS

A virus panel based on the principle of biotinylated RNA baits was developed for specific capture enrichment of epizootic and zoonotic viruses (VirBaits). The VirBaits set was supplemented by a SARS-CoV-2 predesigned bait set for testing recent SARS-CoV-2-positive samples. Libraries generated from complex samples were sequenced via generic HTS (without enrichment) and afterwards enriched with the VirBaits set. For validation, an internal proficiency test for emerging epizootic and zoonotic viruses (African swine fever virus, Ebolavirus, Marburgvirus, Nipah henipavirus, Rift Valley fever virus) was conducted.

RESULTS

The VirBaits set consists of 177,471 RNA baits (80-mer) based on about 18,800 complete viral genomes targeting 35 epizootic and zoonotic viruses. In all tested samples, viruses with both DNA and RNA genomes were clearly enriched ranging from about 10-fold to 10,000-fold for viruses including distantly related viruses with at least 72% overall identity to viruses represented in the bait set. Viruses showing a lower overall identity (38% and 46%) to them were not enriched but could nonetheless be detected based on capturing conserved genome regions. The internal proficiency test supports the improved virus detection using the combination of HTS plus targeted enrichment but also points to the risk of cross-contamination between samples.

CONCLUSIONS

The VirBaits approach showed a high diagnostic performance, also for distantly related viruses. The bait set is modular and expandable according to the favored diagnostics, health sector, or research question. The risk of cross-contamination needs to be taken into consideration. The application of the RNA-baits principle turned out to be user friendly, and even non-experts can easily use the VirBaits workflow. The rapid extension of the established VirBaits set adapted to actual outbreak events is possible as shown for SARS-CoV-2. Video abstract.

Multiple Types of Novel Enteric Bopiviruses (Picornaviridae) with the Possibility of Interspecies Transmission Identified from Cloven-Hoofed Domestic Livestock (Ovine, Caprine and Bovine) in Hungary

Most picornaviruses of the family Picornaviridae are relatively well known, but there are certain “neglected” genera like Bopivirus, containing a single uncharacterised sequence (bopivirus A1, KM589358) with very limited background information. In this study, three novel picornaviruses provisionally called ovipi-, gopi- and bopivirus/Hun (MW298057-MW298059) from enteric samples of asymptomatic ovine, caprine and bovine respectively, were determined using RT-PCR and dye-terminator sequencing techniques. These monophyletic viruses share the same type II-like IRES, NPGP-type 2A, similar genome layout (4-3-4) and cre-localisations. Culture attempts of the study viruses, using six different cell lines, yielded no evidence of viral growth in vitro. Genomic and phylogenetic analyses show that bopivirus/Hun of bovine belongs to the species Bopivirus A, while the closely related ovine-origin ovipi- and caprine-origin gopivirus could belong to a novel species “Bopivirus B” in the genus Bopivirus. Epidemiological investigation of N = 269 faecal samples of livestock (ovine, caprine, bovine, swine and rabbit) from different farms in Hungary showed that bopiviruses were most prevalent among <12-month-old ovine, caprine and bovine, but undetectable in swine and rabbit. VP1 capsid-based phylogenetic analyses revealed the presence of multiple lineages/genotypes, including closely related ovine/caprine strains, suggesting the possibility of ovine–caprine interspecies transmission of certain bopiviruses.

A theoretical and generalized approach for the assessment of the sample-specific limit of detection for clinical metagenomics

Metagenomics is a powerful tool to identify novel or unexpected pathogens, since it is generic and relatively unbiased. The limit of detection (LOD) is a critical parameter for the routine application of methods in the clinical diagnostic context. Although attempts for the determination of LODs for metagenomics next-generation sequencing (mNGS) have been made previously, these were only applicable for specific target species in defined samples matrices. Therefore, we developed and validated a generalized probability-based model to assess the sample-specific LOD of mNGS experiments (LOD_mNGS). Initial rarefaction analyses with datasets of Borna disease virus 1 human encephalitis cases revealed a stochastic behavior of virus read detection. Based on this, we transformed the Bernoulli formula to predict the minimal necessary dataset size to detect one virus read with a probability of 99%. We validated the formula with 30 datasets from diseased individuals, resulting in an accuracy of 99.1% and an average of 4.5 ± 0.4 viral reads found in the calculated minimal dataset size. We demonstrated by modeling the virus genome size, virus-, and total RNA-concentration that the main determinant of mNGS sensitivity is the virus-sample background ratio. The predicted LOD_mNGS for the respective pathogenic virus in the datasets were congruent with the virus-concentration determined by RT-qPCR. Theoretical assumptions were further confirmed by correlation analysis of mNGS and RT-qPCR data from the samples of the analyzed datasets. This approach should guide standardization of mNGS application, due to the generalized concept of LOD_mNGS.

Highly efficient library preparation for Ion Torrent sequencing using Y-adapters

Library preparation is a crucial step in next-generation sequencing workflows. Key determinants of successful library preparation are the available amount of input DNA and the efficiency of the conversion of this DNA into functional library molecules. While the standard blunt-end ligation protocol for Ion Torrent libraries has a theoretical maximum efficiency of 25%, Y-adapters enable highly efficient library preparation by (i) sticky-end ligation and (ii) rendering both DNA strands functional for sequencing, hence resulting in a theoretical efficiency of up to 100%. Moreover, the generation of adapter dimers is reduced. Therefore, we designed, optimized and validated Y-adapters compatible with Ion Torrent sequencing. These facilitate higher library yields combined with overall high sequencing performance regarding the key characteristics read-length, base quality, and library complexity.