Hidden Viral Sequences in Public Sequencing Data and Warning for Future Emerging Diseases

Highly diverse RNA viruses and phage sequences concealed within birds

The diversity of birds in most parts of the world is very high, and thus, they may carry different types of highly differentiated and unknown viruses. Thanks to advanced sequencing technologies, studies on the diversity of bird-associated viruses have increased over the past few years. In this study, a large-scale viral metagenomics survey was performed on cloacal swabs of 2,990 birds from nine provinces of the Chinese mainland. To detect undescribed RNA viruses in birds, more than 1,800 sequences sharing relatively low (<60%) amino acid sequence identity with the best match in the GenBank database were screened. Potentially novel viruses related to vertebrates have been identified, and several potential recombination signals were found. Additionally, hundreds of RNA viral sequences related to plants, fungi, and insects were detected, including previously unknown viruses. Furthermore, we investigated the novelty, functionality, and classification of the phages examined in this study. These viruses occupied topological positions on the evolutionary trees to a certain extent and might form novel putative families, genera, or species, thus providing information to fill the phylogenetic gaps of related viruses. These findings provided new insights into bird-associated viruses, but the interactions among these viruses remain unknown and require further investigation.IMPORTANCEStudying the diversity of RNA viruses in birds and mammals is crucial due to their potential impact on human health and the global ecosystem. Many RNA viruses, such as influenza and coronaviruses, have been shown to cross the species barrier and cause zoonotic diseases. In this metagenomics study involving 2,990 birds from at least 82 species, we identified over 1,800 RNA sequences with distant relationships to known viruses, some of which are rare in birds. The study highlights the scope and diversity of RNA viruses in birds, providing data to predict disease risks and monitor potential viral threats to wildlife, livestock, and human health. This information can aid in the development of strategies for disease prevention and control.

Isolation and Genomic Characterization of Kadipiro Virus from Mosquitoes in Yunnan, China

Background: Kadipiro virus (KDV) is a species of the new 12 segmented RNA virus grouped under the genus Seadornavirus within the Reoviridae family. It has previously been isolated or detected from mosquito, Odonata, and bat feces in Indonesia, China, and Denmark, respectively. Here, we describe the isolation and characterization of a viral strain from mosquitoes in Yunnan Province, China. Methods: Mosquitoes were collected overnight using light traps in Shizong county, on July 17, 2023. Virus was isolated from the mosquito homogenate and grown using baby hamster kidney and Aedes albopictus (C6/36) cells. Preliminary identification of the virus was performed by agarose gel electrophoresis (AGE). The full-genome sequences of the strain were determined by full-length amplification of cDNAs and sequenced using next-generation sequencing. Results: We isolated a viral strain (SZ_M48) from mosquitoes (Culex tritaeniorhynchus Giles) that caused cytopathogenic effects in C6/36 cells. AGE analysis indicated a genome consisting of 12 segments of double-stranded RNA that demonstrated a "6-5-1" pattern, similar to the migrating bands of KDV. Phylogenetic analysis based on the full-genome sequence revealed that SZ_M48 is more clustered with KDV isolates from Hubei and Shangdong in China than with Indonesian and Danish strains. The identity between SZ_M48 and SDKL1625 (Shandong, China) is slightly lower than that of QTM27331 (Hubei, China), and the identity with JKT-7075 (Indonesia) and 21164-6/M.dau/DK (Denmark) is the lowest. Conclusion: The full-genome sequence of the new KDV strain described in this study may be useful for surveillance of the evolutionary characteristics of KDVs. Moreover, these findings extend the knowledge about the genomic diversity, potential vectors, and the distribution of KDVs in China.

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.

Divergent Genotype of Hepatitis A Virus in Alpacas, Bolivia, 2019

Hepatitis A virus (HAV) is a common human pathogen found exclusively in primates. In a molecular and serologic study of 64 alpacas in Bolivia, we detected RNA of distinct HAV in ≈9% of animals and HAV antibodies in ≈64%. Complete-genome analysis suggests a long association of HAV with alpacas.

Metagenomic analysis of herbivorous mammalian viral communities in the Northwest Plateau

BACKGROUND

Mammals are potential hosts for many infectious diseases. However, studies on the viral communities of herbivorous mammals in the Northwest Plateau are limited. Here, we studied the viral communities of herbivorous mammals in the Northwest Plateau using virus metagenomic analysis to analyze and compare the viral community composition of seven animal species.

RESULTS

By library construction and next-generation sequencing, contigs and singlets reads with similar viral sequences were classified into 24 viral families. Analyzed from the perspective of sampling areas, the virus community composition was relatively similar in two areas of Wuwei and Jinchang, Gansu Province. Analyzed from the perspective of seven animal species, the viral reads of seven animal species were mostly ssDNA and dominated by CRESS-DNA viruses. Phylogenetic analysis based on viral marker genes indicated that CRESS-DNA viruses and microviruses have high genetic diversity. In addition to DNA viruses, nodaviruses, pepper mild mottle viruses and picornaviruses were RNA viruses that we performed by phylogenetic analysis. The CRESS-DNA viruses and nodaviruses are believed to infect plants and insects, and microviruses can infect bacteria, identifying that they were likely from the diet of herbivorous mammals. Notably, two picornaviruses were identified from red deer and wild horse, showing that the picornavirus found in red deer had the relatively high similarity with human hepatitis A virus, and the picornavirus carried by wild horse could potentially form a new species within the Picornaviridae family.

CONCLUSIONS

This study explored the herbivorous mammalian virus community in the Northwest Plateau and the genetic characteristics of viruses that potentially threaten human health. It reveals the diversity and stability of herbivorous mammalian virus communities in the Northwest Plateau and helps to expand our knowledge of various herbivorous mammalian potentially pathogenic viruses.

Benchmarking State-of-the-Art Approaches for Norovirus Genome Assembly in Metagenome Sample

A recently published article in BMCGenomics by Fuentes-Trillo et al. contains a comparison of assembly approaches of several noroviral samples via different tools and preprocessing strategies. It turned out that the study used outdated versions of tools as well as tools that were not designed for the viral assembly task. In order to improve the suboptimal assemblies, authors suggested different sophisticated preprocessing strategies that seem to make only minor contributions to the results. We have reproduced the analysis using state-of-the-art tools designed for viral assembly, and we demonstrate that tools from the SPAdes toolkit (rnaviralSPAdes and coronaSPAdes) allow one to assemble the samples from the original study into a single contig without any additional preprocessing.

Attenuation hotspots in neurotropic human astroviruses

During the last decade, the detection of neurotropic astroviruses has increased dramatically. The MLB genogroup of astroviruses represents a genetically distinct group of zoonotic astroviruses associated with gastroenteritis and severe neurological complications in young children, the immunocompromised, and the elderly. Using different virus evolution approaches, we identified dispensable regions in the 3' end of the capsid-coding region responsible for attenuation of MLB astroviruses in susceptible cell lines. To create recombinant viruses with identified deletions, MLB reverse genetics (RG) and replicon systems were developed. Recombinant truncated MLB viruses resulted in imbalanced RNA synthesis and strong attenuation in iPSC-derived neuronal cultures confirming the location of neurotropism determinants. This approach can be used for the development of vaccine candidates using attenuated astroviruses that infect humans, livestock animals, and poultry.

Prevalence and genetic diversity in bovine parechovirus infecting Japanese cattle

The first bovine parechovirus (Bo_ParV) was reported in 2021, and currently, only two nearly complete genome sequences of Bo_ParV are available. In this study, we detected Bo_ParVs in 10 out of 158 bovine fecal samples tested using real-time RT-PCR, and Bo_ParVs were isolated from three of these samples using MA104 cells. Analysis of the P1 region revealed that Bo_ParVs shared high pairwise amino acid sequence similarity (≥ 95.7% identity), suggesting antigenic similarity among Bo_ParVs, whereas nucleotide sequence identity values (≥ 84.8%) indicated more variability. A recombination breakpoint was identified in the 2B region, which may influence the evolution of this virus.

The First Identification of a Narnavirus in Bigyra, a Marine Protist

Current information on the diversity and evolution of eukaryotic RNA viruses is biased towards host lineages, such as animals, plants, and fungi. Although protists represent the majority of eukaryotic diversity, our understanding of the protist RNA virosphere is still limited. To reveal untapped RNA viral diversity, we screened RNA viruses from 30 marine protist isolates and identified a novel RNA virus named Haloplacidia narnavirus 1 (HpNV1). A phylogenetic ana-lysis revealed that HpNV1 is a new member of the family Narnaviridae. The present study filled a gap in the distribution of narnaviruses and implies their wide distribution in Stramenopiles.

Discovery of novel papillomaviruses in the critically endangered Malayan and Chinese pangolins

Pangolins are scaly and toothless mammals which are distributed across Africa and Asia. Currently, the Malayan, Chinese and Philippine pangolins are designated as critically endangered species. Although few pangolin viruses have been described, their viromes have received more attention following the discovery that they harbour sarbecoviruses related to SARS-CoV-2. Using large-scale genome mining, we discovered novel lineages of papillomaviruses infecting the Malayan and Chinese pangolins. We were able to assemble three complete circular papillomavirus genomes with an intact coding capacity and five additional L1 genes encoding the major capsid protein. Phylogenetic analysis revealed that seven out of eight L1 sequences formed a monophyletic group which is the sister lineage to the Tupaia belangeri papillomavirus 1, isolated from Yunnan province in China. Additionally, a single L1 sequence assembled from a Chinese pangolin was placed in a clade closer to Alphapapillomavirus and Omegapapillomavirus. Examination of the SRA data from 95 re-sequenced genomes revealed that 49.3% of Malayan pangolins and 50% of Chinese pangolins were positive for papillomavirus reads. Our results indicate that pangolins in South-East Asia are the hosts of diverse and highly prevalent papillomaviruses, and highlight the value of in silico mining of host sequencing data for the discovery of novel viruses.

Large-scale investigation of zoonotic viruses in the era of high-throughput sequencing

Zoonotic diseases considerably impact public health and socioeconomics. RNA viruses reportedly caused approximately 94% of zoonotic diseases documented from 1990 to 2010, emphasizing the importance of investigating RNA viruses in animals. Furthermore, it has been estimated that hundreds of thousands of animal viruses capable of infecting humans are yet to be discovered, warning against the inadequacy of our understanding of viral diversity. High-throughput sequencing (HTS) has enabled the identification of viral infections with relatively little bias. Viral searches using both symptomatic and asymptomatic animal samples by HTS have revealed hidden viral infections. This review introduces the history of viral searches using HTS, current analytical limitations, and future potentials. We primarily summarize recent research on large-scale investigations on viral infections reusing HTS data from public databases. Furthermore, considering the accumulation of uncultivated viruses, we discuss current studies and challenges for connecting viral sequences to their phenotypes using various approaches: performing data analysis, developing predictive modeling, or implementing high-throughput platforms of virological experiments. We believe that this article provides a future direction in large-scale investigations of potential zoonotic viruses using the HTS technology.

Is the Intergenic Region of Aedes aegypti Totivirus a Recombination Hotspot?

The genus totivirus in the family Totiviridae contains double-stranded RNA viruses. Their genome has two open reading frames (ORFs) that encode capsid protein (CP) and RNA-dependent RNA polymerase (RdRp). The toti-like viruses recently identified in Anopheles sp. and Aedes aegypti mosquitoes (AaTV) share the same genome organization as other totiviruses. The AaTVs that have been described in distinct geographical regions are monophyletic. In this study, we show that AaTV sequences can be grouped into at least three phylogenetic clades (named A, B, and C). Clades A and B are composed of AaTV sequences from mosquitoes collected in the Caribbean region (Guadeloupe), and clade C contains sequences from the USA. These clades may represent AaTV lineages that are locally adapted to their host populations. We also identified three recombinant AaTV strains circulating in mosquitoes in Guadeloupe. Although these strains have different chimeric patterns, the position of the recombination breakpoint was identical in all strains. Interestingly, this breakpoint is located in a hairpin-like structure in the intergenic region of the AaTV genome. This RNA structure may stall RNA polymerase processivity and consequently induce template switching. In vitro studies should be conducted to further investigate the biological significance of AaTV's intergenic region as a recombination hotspot.

First isolation and genomic characterization of bovine parechovirus from faecal samples of cattle in Japan

A novel picornavirus was isolated from the faeces of a diarrhoeic cow using MA-104 cells at the third blind passage. This virus, named Den1/2021/JPN, was completely sequenced using total RNA from the cell culture supernatant by deep sequencing. The genome of Den1/2021/JPN had a standard picornavirus genome organisation with conserved picornaviral motifs. The 5' untranslated region harboured a type-II internal ribosomal entry site. Den1/2021/JPN was most closely related to a bovine parechovirus (Bo_ParV) named cow/2018/4, which has been recently identified in publicly available databases. Phylogenetic analyses and pairwise sequence comparison revealed that Den1/2021/JPN and Bo_ParV cow/2018/4 clustered with parechoviruses and were most closely related to Parechovirus E identified in birds of prey, exhibiting nucleotide sequence similarity of 64.2-64.5 %, 58.6-59.7 % and 66.3-66.4 % in the polyprotein, P1 and 2C+3 CD coding regions, respectively. This study presents the first report on the isolation of Bo_ParV. Den1/2021/JPN and Bo_ParV cow/2018/4, which are candidates for a novel species in the genus Parechovirus.

Association of Common Zoonotic Pathogens With Concentrated Animal Feeding Operations

Animal farming has intensified significantly in recent decades, with the emergence of concentrated animal feeding operations (CAFOs) in industrialized nations. The congregation of susceptible animals in CAFOs can lead to heavy environmental contamination with pathogens, promoting the emergence of hyper-transmissible, and virulent pathogens. As a result, CAFOs have been associated with emergence of highly pathogenic avian influenza viruses, hepatitis E virus, Escherichia coli O157:H7, Streptococcus suis, livestock-associated methicillin-resistant Staphylococcus aureus, and Cryptosporidium parvum in farm animals. This has led to increased transmission of zoonotic pathogens in humans and changes in disease patterns in general communities. They are exemplified by the common occurrence of outbreaks of illnesses through direct and indirect contact with farm animals, and wide occurrence of similar serotypes or subtypes in both humans and farm animals in industrialized nations. Therefore, control measures should be developed to slow down the dispersal of zoonotic pathogens associated with CAFOs and prevent the emergence of new pathogens of epidemic and pandemic potential.

Novel hepatoviruses in synanthropic bats in the upper Midwestern United States

A nearly complete genome sequence of hepatovirus G was isolated from an Eptesicus fuscus bat submitted for rabies virus testing due to human exposure in South Dakota. The predicted polyprotein sequence was 78.2% and 74.4% identical to genotypes G1 and G2, respectively, recovered from bats in Ghana. Quantitative PCR on 90 E. fuscus bats showed that eight (8.9%) were positive for hepatovirus G. Targeted sequencing of the VP2 region of the genome for five positive samples showed >99% identity to hepatovirus G strain Ef15893, demonstrating that hepatovirus G commonly circulates in E. fuscus bats in the upper Midwest.

[Yezo virus and emerging orthonairovirus diseases]

A new etiological agent of an acute febrile illness following tick bite has been found in Hokkaido, Japan, in 2019 and designated as Yezo virus. Seven cases of Yezo virus infection were identified from 2014 to 2020 by passive and retrospective surveillance. Yezo virus is classified into the genus Orthonairovirus, family Nairoviridae and forms Sulina genogroup together with Sulina virus, which was identified in ticks in Romania. The Sulina genogroup viruses are closely related to the Tamdy genogroup viruses recently reported as causative agents of febrile illness in China and distant from known orthonairovirus pathogens, such as Crimean-Congo hemorrhagic fever virus. Since only limited information is available for the emerging orthonairovirus diseases, including Yezo virus infection, their occurrence should be carefully monitored.