Novel viruses discovered in metatranscriptomic analysis of farmed barramundi in Asia and Australia.

Barramundi aquaculture is at risk of severe disease outbreaks and massive production losses. Here we used bioinformatics to screen 84 farmed barramundi transcriptomes to identify novel viruses that could threaten barramundi aquaculture and to establish a barramundi aquaculture virome. We discovered five novel viruses: latid herpesvirus 1 (LatHV-1) from the Alloherpesviridae family, barramundi parvovirus 1 (BParV1) from the Parvoviridae family, barramundi calicivirus 1 (BCaV1) from the Caliciviridae family, and barramundi associated picorna-like virus 1 and 2 (BPicV1 and BPicV2) from the Picornaviridae family. LatHV-1, BCaV1, and BParV1 are closely related to pathogenic viruses found in other fish species that can cause mass mortality in farms. To aid in future viral surveillance, we also designed and successfully tested an RT-PCR assay for the detection of BCaV1. Overall, we discovered a range of pathogenic viruses in barramundi aquaculture, paving the way for developing effective detection methods to assist early outbreak management.

A Statovirus-like virus from respiratory tracts of patients, China.

The emerging evidence of human infections with emerging viruses suggests their potential public health importance. A novel taxon of viruses named Statoviruses (for stool-associated Tombus-like viruses) was recently identified in the gastrointestinal tracts of multiple mammals. Here we report the discovery of respiratory Statovirus-like viruses (provisionally named Restviruses) from the respiratory tracts of five patients experiencing acute respiratory disease with Human coronavirus OC43 infection through the retrospective analysis of meta-transcriptomic data. Restviruses shared 53.1%-98.8% identities of genomic sequences with each other and 39.9%-44.3% identities with Statoviruses. The phylogenetic analysis revealed that Restviruses together with a Stato-like virus from nasal-throat swabs of Vietnamese patients with acute respiratory disease, formed a well-supported clade distinct from the taxon of Statoviruses. However, the consistent genome characteristics of Restviruses and Statoviruses suggested that they might share similar evolutionary trajectories. These findings warrant further studies to elucidate the etiological and epidemiological significance of the emerging Restviruses.

Benchmarking informatics approaches for virus discovery: caution is needed when combining in silico identification methods.

Understanding the ecological impacts of viruses on natural and engineered ecosystems relies on the accurate identification of viral sequences from community sequencing data. To maximize viral recovery from metagenomes, researchers frequently combine viral identification tools. However, the effectiveness of this strategy is unknown. Here, we benchmarked combinations of six widely used informatics tools for viral identification and analysis (VirSorter, VirSorter2, VIBRANT, DeepVirFinder, CheckV, and Kaiju), called "rulesets." Rulesets were tested against mock metagenomes composed of taxonomically diverse sequence types and diverse aquatic metagenomes to assess the effects of the degree of viral enrichment and habitat on tool performance. We found that six rulesets achieved equivalent accuracy [Matthews Correlation Coefficient (MCC) = 0.77, Padj ≥ 0.05]. Each contained VirSorter2, and five used our "tuning removal" rule designed to remove non-viral contamination. While DeepVirFinder, VIBRANT, and VirSorter were each found once in these high-accuracy rulesets, they were not found in combination with each other: combining tools does not lead to optimal performance. Our validation suggests that the MCC plateau at 0.77 is partly caused by inaccurate labeling within reference sequence databases. In aquatic metagenomes, our highest MCC ruleset identified more viral sequences in virus-enriched (44%-46%) than in cellular metagenomes (7%-19%). While improved algorithms may lead to more accurate viral identification tools, this should be done in tandem with careful curation of sequence databases. We recommend using the VirSorter2 ruleset and our empirically derived tuning removal rule. Our analysis provides insight into methods for in silico viral identification and will enable more robust viral identification from metagenomic data sets. IMPORTANCE: The identification of viruses from environmental metagenomes using informatics tools has offered critical insights in microbial ecology. However, it remains difficult for researchers to know which tools optimize viral recovery for their specific study. In an attempt to recover more viruses, studies are increasingly combining the outputs from multiple tools without validating this approach. After benchmarking combinations of six viral identification tools against mock metagenomes and environmental samples, we found that these tools should only be combined cautiously. Two to four tool combinations maximized viral recovery and minimized non-viral contamination compared with either the single-tool or the five- to six-tool ones. By providing a rigorous overview of the behavior of in silico viral identification strategies and a pipeline to replicate our process, our findings guide the use of existing viral identification tools and offer a blueprint for feature engineering of new tools that will lead to higher-confidence viral discovery in microbiome studies.

Untangling an insect's virome from its endogenous viral elements.

BACKGROUND: Insects are an important reservoir of viral biodiversity, but the vast majority of viruses associated with insects have not been discovered. Recent studies have employed high-throughput RNA sequencing, which has led to rapid advances in our understanding of insect viral diversity. However, insect genomes frequently contain transcribed endogenous viral elements (EVEs) with significant homology to exogenous viruses, complicating the use of RNAseq for viral discovery. METHODS: In this study, we used a multi-pronged sequencing approach to study the virome of an important agricultural pest and prolific vector of plant pathogens, the potato aphid Macrosiphum euphorbiae. We first used rRNA-depleted RNAseq to characterize the microbes found in individual insects. We then used PCR screening to measure the frequency of two heritable viruses in a local aphid population. Lastly, we generated a quality draft genome assembly for M. euphorbiae using Illumina-corrected Nanopore sequencing to identify transcriptionally active EVEs in the host genome. RESULTS: We found reads from two insect-specific viruses (a Flavivirus and an Ambidensovirus) in our RNAseq data, as well as a parasitoid virus (Bracovirus), a plant pathogenic virus (Tombusvirus), and two phages (Acinetobacter and APSE). However, our genome assembly showed that part of the 'virome' of this insect can be attributed to EVEs in the host genome. CONCLUSION: Our work shows that EVEs have led to the misidentification of aphid viruses from RNAseq data, and we argue that this is a widespread challenge for the study of viral diversity in insects.

The RNA virome of echinoderms.

Echinoderms are a phylum of marine invertebrates that include model organisms, keystone species, and animals commercially harvested for seafood. Despite their scientific, ecological, and economic importance, there is little known about the diversity of RNA viruses that infect echinoderms compared to other invertebrates. We screened over 900 transcriptomes and viral metagenomes to characterize the RNA virome of 38 echinoderm species from all five classes (Crinoidea, Holothuroidea, Asteroidea, Ophiuroidea and Echinoidea). We identified 347 viral genome fragments that were classified to genera and families within nine viral orders - Picornavirales, Durnavirales, Martellivirales, Nodamuvirales, Reovirales, Amarillovirales, Ghabrivirales, Mononegavirales, and Hepelivirales. We compared the relative viral representation across three life stages (embryo, larvae, adult) and characterized the gene content of contigs which encoded complete or near-complete genomes. The proportion of viral reads in a given transcriptome was not found to significantly differ between life stages though the majority of viral contigs were discovered from transcriptomes of adult tissue. This study illuminates the biodiversity of RNA viruses from echinoderms, revealing the occurrence of viral groups in natural populations.

Adenovirus surveillance in wild carnivores from Brazil.

Landscape transformation favors the spread of new pathogens that can be shared between domestic and wild animals. Certain adenoviruses (e.g., canine adenovirus 1 and 2, family Adenoviridae) can infect domestic and wild carnivores. In domestic canids, these viruses are associated with hepatic and respiratory diseases (among others). Nevertheless, information regarding adenovirus pathogenicity and molecular features in wild carnivores is still limited. Herein we surveyed adenovirus in free-ranging carnivores from Brazil. Total DNA was extracted from and subsequently tested by a nested panPCR in spleen and/or lung of 52 carnivores, representing species of the following families: Canidae (n = 4), Felidae (n = 3), Mustelidae (n = 2) and Procyonidae (n = 2). The obtained sequences were compared to others available at GenBank. Available tissue samples from the positive cases were evaluated histopathologically. One out of 52 (1.9%, CI 95%, 0.0-5.7%) carnivores was positive; a roadkilled ocelot (Leopardus pardalis). The obtained sequence presented a low deduced amino acid (78.1%) similarity with the closest adenovirus, identified in a pinniped from the United States of America. This fact and its detection in a novel host suggest it may be representative of a novel species and denominated ocelot adenovirus 1. None of the gross and microscopic findings of the positive case were associated with adenovirus. To the authors' knowledge, this is the first report of adenovirus in wild felids of South America and the second worldwide. Further studies are necessary to assess the epidemiology and potential pathogenicity of this agent in wild carnivores.

Discovery of Novel Viruses Associated With the Invasive Cane Toad (Rhinella marina) in Its Native and Introduced Ranges.

Cane toads (Rhinella marina) are notoriously successful invaders: from 101 individuals brought to Australia in 1935, poisonous toads now cover an area >1.2 million km2 with adverse effects on native fauna. Despite extensive research on the role of macroparasites in cane toad invasion, viral research is lagging. We compared viral prevalence and diversity between toads in their native range (French Guiana, n=25) and two introduced ranges: Australia (n=151) and Hawai'i (n=10) with a metatranscriptomic and metagenomic approach combined with PCR screening. Australian toads almost exclusively harbor one of seven viruses detected globally. Rhimavirus-A (Picornaviridae) exhibited low genetic diversity and likely actively infected 9% of sampled Australian toads extending across ~2,000km of Northern Australia and up to the current invasion front. In native range cane toads, we identified multiple phylogenetically distinct viruses (Iridoviridae, Picornaviridae, Papillomaviridae, and Nackedna-like virus). None of the same viruses was detected in both ranges, suggesting that Australian cane toads have largely escaped the viral infection experienced by their native range counterparts. The novel native range viruses described here are potential biocontrol agents, as Australian toads likely lack prior immunological exposure to these viruses. Overall, our evidence suggests that there may be differences between viruses infecting cane toads in their native vs. introduced ranges, which lays the groundwork for further studies on how these viruses have influenced the toads' invasion history.

Single mosquito metatranscriptomics identifies vectors, emerging pathogens and reservoirs in one assay.

Mosquitoes are major infectious disease-carrying vectors. Assessment of current and future risks associated with the mosquito population requires knowledge of the full repertoire of pathogens they carry, including novel viruses, as well as their blood meal sources. Unbiased metatranscriptomic sequencing of individual mosquitoes offers a straightforward, rapid, and quantitative means to acquire this information. Here, we profile 148 diverse wild-caught mosquitoes collected in California and detect sequences from eukaryotes, prokaryotes, 24 known and 46 novel viral species. Importantly, sequencing individuals greatly enhanced the value of the biological information obtained. It allowed us to (a) speciate host mosquito, (b) compute the prevalence of each microbe and recognize a high frequency of viral co-infections, (c) associate animal pathogens with specific blood meal sources, and (d) apply simple co-occurrence methods to recover previously undetected components of highly prevalent segmented viruses. In the context of emerging diseases, where knowledge about vectors, pathogens, and reservoirs is lacking, the approaches described here can provide actionable information for public health surveillance and intervention decisions.

ContigExtender: a new approach to improving de novo sequence assembly for viral metagenomics data.

BACKGROUND: Metagenomics is the study of microbial genomes for pathogen detection and discovery in human clinical, animal, and environmental samples via Next-Generation Sequencing (NGS). Metagenome de novo sequence assembly is a crucial analytical step in which longer contigs, ideally whole chromosomes/genomes, are formed from shorter NGS reads. However, the contigs generated from the de novo assembly are often very fragmented and rarely longer than a few kilo base pairs (kb). Therefore, a time-consuming extension process is routinely performed on the de novo assembled contigs. RESULTS: To facilitate this process, we propose a new tool for metagenome contig extension after de novo assembly. ContigExtender employs a novel recursive extending strategy that explores multiple extending paths to achieve highly accurate longer contigs. We demonstrate that ContigExtender outperforms existing tools in synthetic, animal, and human metagenomics datasets. CONCLUSIONS: A novel software tool ContigExtender has been developed to assist and enhance the performance of metagenome de novo assembly. ContigExtender effectively extends contigs from a variety of sources and can be incorporated in most viral metagenomics analysis pipelines for a wide variety of applications, including pathogen detection and viral discovery.

Diversity of Sea Star-Associated Densoviruses and Transcribed Endogenous Viral Elements of Densovirus Origin.

A viral etiology of sea star wasting syndrome (SSWS) was originally explored with virus-sized material challenge experiments, field surveys, and metagenomics, leading to the conclusion that a densovirus is the predominant DNA virus associated with this syndrome and, thus, the most promising viral candidate pathogen. Single-stranded DNA viruses are, however, highly diverse and pervasive among eukaryotic organisms, which we hypothesize may confound the association between densoviruses and SSWS. To test this hypothesis and assess the association of densoviruses with SSWS, we compiled past metagenomic data with new metagenomic-derived viral genomes from sea stars collected from Antarctica, California, Washington, and Alaska. We used 179 publicly available sea star transcriptomes to complement our approaches for densovirus discovery. Lastly, we focus the study on sea star-associated densovirus (SSaDV), the first sea star densovirus discovered, by documenting its biogeography and putative tissue tropism. Transcriptomes contained only endogenized densovirus elements similar to the NS1 gene, while numerous extant densoviral genomes were recovered from viral metagenomes. SSaDV was associated with nearly all tested species from southern California to Alaska, and in contrast to previous work, we show that SSaDV is one genotype among a high diversity of densoviruses present in sea stars across the West Coast of the United States and globally that are commonly associated with grossly normal (i.e., healthy or asymptomatic) animals. The diversity and ubiquity of these viruses in sea stars confound the original hypothesis that one densovirus is the etiological agent of SSWS.IMPORTANCE The primary interest in sea star densoviruses, specifically SSaDV, has been their association with sea star wasting syndrome (SSWS), a disease that has decimated sea star populations across the West Coast of the United States since 2013. The association of SSaDV with SSWS was originally drawn from metagenomic analysis, which was further studied through field surveys using quantitative PCR (qPCR), with the conclusion that it was the most likely viral candidate in the metagenomic data based on its representation in symptomatic sea stars compared to asymptomatic sea stars. We reexamined the original metagenomic data with additional genomic data sets and found that SSaDV was 1 of 10 densoviruses present in the original data set and was no more represented in symptomatic sea stars than in asymptomatic sea stars. Instead, SSaDV appears to be a widespread, generalist virus that exists among a large diversity of densoviruses present in sea star populations.

Characterization of a novel picornavirus isolated from moribund aquacultured clownfish.

Over the last decade, a number of USA aquaculture facilities have experienced periodic mortality events of unknown aetiology in their clownfish (Amphiprion ocellaris). Clinical signs of affected individuals included lethargy, altered body coloration, reduced body condition, tachypnea, and abnormal positioning in the water column. Samples from outbreaks were processed for routine parasitological, bacteriological, and virological diagnostic testing, but no consistent parasitic or bacterial infections were observed. Histopathological evaluation revealed individual cell necrosis and mononuclear cell inflammation in the branchial cavity, pharynx, oesophagus and/or stomach of four examined clownfish, and large basophilic inclusions within the pharyngeal mucosal epithelium of one fish. Homogenates from pooled external and internal tissues from these outbreaks were inoculated onto striped snakehead (SSN-1) cells for virus isolation and cytopathic effects were observed, resulting in monolayer lysis in the initial inoculation and upon repassage. Transmission electron microscopy of infected SSN-1 cells revealed small round particles (mean diameter=20.0-21.7 nm) within the cytoplasm, consistent with the ultrastructure of a picornavirus. Full-genome sequencing of the purified virus revealed a novel picornavirus most closely related to the bluegill picornavirus and other members of the genus Limnipivirus. Additionally, pairwise protein alignments between the clownfish picornavirus (CFPV) and other known members of the genus Limnipivirus yielded results in accordance with the current International Committee on Taxonomy of Viruses criteria for members of the same genus. Thus, CFPV represents a proposed new limnipivirus species. Future experimental challenge studies are needed to determine the role of CFPV in disease.

A Highly Prevalent and Pervasive Densovirus Discovered among Sea Stars from the North American Atlantic Coast.

The etiology of sea star wasting syndrome is hypothesized to be caused by a densovirus, sea star-associated densovirus (SSaDV), that has previously been reported on the Pacific and Atlantic Coasts of the United States. In this study, we reevaluated the presence of SSaDV among sea stars from the North American Atlantic Coast and in doing so discovered a novel densovirus that we have named Asterias forbesi-associated densovirus (AfaDV), which shares 78% nucleotide pairwise identity with SSaDV. In contrast to previous studies, SSaDV was not detected in sea stars from the North American Atlantic Coast. Using a variety of PCR-based techniques, we investigated the tissue tropism, host specificity, and prevalence of AfaDV among populations of sea stars at five locations along the Atlantic Coast. AfaDV was detected in three sea star species (Asterias forbesi, Asterias rubens, and Henricia sp.) found in this region and was highly prevalent (>80% of individuals tested; n = 134), among sampled populations. AfaDV was detected in the body wall, gonads, and pyloric caeca (digestive gland) of specimens but was not detected in their coelomic fluid. A significant difference in viral load (copies mg-1) was found between tissue types, with the pyloric caeca having the highest viral loads. Further investigation of Asterias forbesi gonad tissue found germ line cells (oocytes) to be virus positive, suggesting a potential route of vertical transmission. Taken together, these observations show that the presence of AfaDV is not an indicator of sea star wasting syndrome because AfaDV is a common constituent of these animals' microbiome, regardless of health.IMPORTANCE Sea star wasting syndrome is a disease primarily observed on the Pacific and Atlantic Coasts of North America that has significantly impacted sea star populations. The etiology of this disease is unknown, although it is hypothesized to be caused by a densovirus, SSaDV. However, previous studies have not found a correlation between SSaDV and sea star wasting syndrome on the North American Atlantic Coast. This study suggests that this observation may be explained by the presence of a genetically similar densovirus, AfaDV, that may have confounded previous studies. SSaDV was not present in sea stars screened in this study, and instead, AfaDV was commonly found in sea star populations across the New England region, with no apparent signs of disease. These results suggest that sea star densoviruses may be common constituents of the animals' microbiome, and the diversity and extent of these viruses among wild populations may be greater than previously recognized.

Mash Screen: high-throughput sequence containment estimation for genome discovery.

The MinHash algorithm has proven effective for rapidly estimating the resemblance of two genomes or metagenomes. However, this method cannot reliably estimate the containment of a genome within a metagenome. Here, we describe an online algorithm capable of measuring the containment of genomes and proteomes within either assembled or unassembled sequencing read sets. We describe several use cases, including contamination screening and retrospective analysis of metagenomes for novel genome discovery. Using this tool, we provide containment estimates for every NCBI RefSeq genome within every SRA metagenome and demonstrate the identification of a novel polyomavirus species from a public metagenome.

High Diversity and Novel Enteric Viruses in Fecal Viromes of Healthy Wild and Captive Thai Cynomolgus Macaques (Macaca fascicularis).

Cynomolgus macaques are common across South East Asian countries including Thailand. The National Primate Research Center of Thailand, Chulalongkorn University (NPRCT-CU) captures wild-borne cynomolgus macaque for research use. Limited information is available on the enteric viruses and possible zoonotic infections into or from cynomolgus macaques. We characterized and compare the fecal virome of two populations; healthy wild-originated captive cynomolgus macaques (n = 43) reared in NPRCT-CU and healthy wild cynomolgus macaques (n = 35). Over 90% of recognized viral sequence reads amplified from feces were from bacterial viruses. Viruses from seven families of mammalian viruses were also detected (Parvoviridae, Anelloviridae, Picornaviridae, Adenoviridae, Papillomaviridae, Herpesviridae, and Caliciviridae). The genomes of a member of a new picornavirus genus we named Mafapivirus, a primate chapparvovirus, and a circular Rep-encoding single-strand (CRESS) DNA virus were also characterized. Higher abundance of CRESS DNA viruses of unknown tropism and invertebrate-tropic ambidensovirus were detected in wild versus captive macaques likely reflecting dietary differences. Short term rearing in captivity did not have a pronounced effect on the diversity of mammalian viruses of wild cynomolgus macaques. This study is the first report of the fecal virome of cynomolgus macaques, non-human primates frequently used in biomedical research and vaccination studies.

A decade of RNA virus metagenomics is (not) enough.

It is hard to overemphasize the role that metagenomics has had on our recent understanding of RNA virus diversity. Metagenomics in the 21st century has brought with it an explosion in the number of RNA virus species, genera, and families far exceeding that following the discovery of the microscope in the 18th century for eukaryotic life or culture media in the 19th century for bacteriology or the 20th century for virology. When the definition of success in organism discovery is measured by sequence diversity and evolutionary distance, RNA viruses win. This review explores the history of RNA virus metagenomics, reasons for the successes so far in RNA virus metagenomics, and methodological concerns. In addition, the review briefly covers clinical metagenomics and environmental metagenomics and highlights some of the critical accomplishments that have defined the fast pace of RNA virus discoveries in recent years. Slightly more than a decade in, the field is exhausted from its discoveries but knows that there is yet even more out there to be found.

Genomes of viral isolates derived from different mosquitos species.

Eleven viral isolates derived mostly in albopictus C6/36 cells from mosquito pools collected in Southeast Asia and the Americas between 1966 and 2014 contained particles with electron microscopy morphology typical of reoviruses. Metagenomics analysis yielded the near complete genomes of three novel reoviruses, Big Cypress orbivirus, Ninarumi virus, and High Island virus and a new tetravirus, Sarawak virus. Strains of previously characterized Sathuvarachi, Yunnan, Banna and Parry's Lagoon viruses (Reoviridae), Bontang virus (Mesoniviridae), and Culex theileri flavivirus (Flaviviridae) were also characterized. The availability of these mosquito virus genomes will facilitate their detection by metagenomics or PCR to better determine their geographic range, extent of host tropism, and possible association with arthropod or vertebrate disease.

The potential of viral metagenomics in blood transfusion safety.

Thanks to the significant advent of high throughput sequencing in the last ten years, it is now possible via metagenomics to define the spectrum of the microbial sequences present in human blood samples. Therefore, metagenomics sequencing appears as a promising approach for the identification and global surveillance of new, emerging and/or unexpected viruses that could impair blood transfusion safety. However, despite considerable advantages compared to the traditional methods of pathogen identification, this non-targeted approach presents several drawbacks including a lack of sensitivity and sequence contaminant issues. With further improvements, especially to increase sensitivity, metagenomics sequencing should become in a near future an additional diagnostic tool in infectious disease field and especially in blood transfusion safety.

Statoviruses, A novel taxon of RNA viruses present in the gastrointestinal tracts of diverse mammals.

Next-generation sequencing has expanded our understanding of the viral populations that constitute the mammalian virome. We describe a novel taxon of viruses named Statoviruses, for Stool associated Tombus-like viruses, present in multiple metagenomic datasets. These viruses define a novel clade that is phylogenetically related to the RNA virus families Tombusviridae and Flaviviridae. Five distinct statovirus types were identified in human, macaque, mouse, and cow gastrointestinal tract samples. The prototype genome, statovirus A, was frequently identified in macaque stool samples from multiple geographically distinct cohorts. Another genome, statovirus C1, was discovered in a stool sample from a human child with fever, cough, and rash. Further experimental data will clarify whether these viruses are infectious to mammals or if they originate from another source present in the mammalian gastrointestinal tract.

Development of a high-throughput multiplexed real time RT-PCR assay for detection of human pegivirus 1 and 2.

Human Pegivirus 2 (HPgV-2) was recently identified in the bloodstream of HCV-infected and multiply transfused individuals. Initial reports show HPgV-2 circulates at a low prevalence in HCV co-infected individuals, necessitating testing of large cohorts of samples to identify infected persons. The identification of additional HPgV-2 cases was facilitated by the development of a high throughput and reliable molecular reverse transcription polymerase chain reaction (RT-PCR) assay intended for use on the automated Abbott m2000 system with a capability of extracting and testing 96 samples at once. A dual target approach was taken to reduce the risk of a false-negative result, amplifying sequences within the 5' UTR and NS2/3 coding regions of HPgV-2. The assay was expanded to multiplex detection of the other human Pegivirus, HPgV-1 (formerly GBV-C), to allow simultaneous prevalence comparison. The limit of detection (LOD; 95% detection) for HPgV-2 was experimentally determined to be 126 copies/mL. Through use of the newly developed multiplex assay, 21 strains of HPgV-2 circulating in HCV past or present infections were identified, with all strains confirmed by next generation sequencing. The multiplexed assay has high specificity and showed no cross-reactivity of HPgV-2 with HPgV-1 or other Flaviviruses. This automated assay will be instrumental in future studies addressing HPgV-2 pathogenicity, prevalence, and sequence diversity.