Comparing viral metagenomics methods using a highly multiplexed human viral pathogens reagent

Analysis of viral diversity in dogs with acute gastroenteritis from Brazilian Amazon

Viral gastroenteritis is commonly reported in dogs and involves a great diversity of enteric viruses. In this research, viral diversity was investigated in dogs with diarrhea in Northern Brazil using shotgun metagenomics. Furthermore, the presence of norovirus (NoV) was investigated in 282 stool/rectal swabs of young/adult dogs with or without diarrhea from two public kennels, based on one-step reverse transcription polymerase chain reaction (RT-PCR) for genogroup VI and VII (GVI and GVII) and real-time RT-PCR for GI, GII, and GIV. Thirty-one viral families were identified, including bacteriophages. Phylogenetic analyses showed twelve complete or nearly complete genomes belonging to the species of Protoparvovirus carnivoran1, Mamastrovirus 5, Aichivirus A2, Alphacoronavirus 1, and Chipapillomavirus 1. This is the first description of the intestinal virome of dogs in Northern Brazil and the first detection of canine norovirus GVII in the country. These results are important for helping to understand the viral groups that circulate in the canine population.

Exploring viral diversity and metagenomics in livestock: insights into disease emergence and spillover risks in cattle

Cattle have a significant impact on human societies in terms of both economics and health. Viral infections pose a relevant problem as they directly or indirectly disrupt the balance within cattle populations. This has negative consequences at the economic level for producers and territories, and also jeopardizes human health through the transmission of zoonotic diseases that can escalate into outbreaks or pandemics. To establish prevention strategies and control measures at various levels (animal, farm, region, or global), it is crucial to identify the viral agents present in animals. Various techniques, including virus isolation, serological tests, and molecular techniques like PCR, are typically employed for this purpose. However, these techniques have two major drawbacks: they are ineffective for non-culturable viruses, and they only detect a small fraction of the viruses present. In contrast, metagenomics offers a promising approach by providing a comprehensive and unbiased analysis for detecting all viruses in a given sample. It has the potential to identify rare or novel infectious agents promptly and establish a baseline of healthy animals. Nevertheless, the routine application of viral metagenomics for epidemiological surveillance and diagnostics faces challenges related to socioeconomic variables, such as resource availability and space dedicated to metagenomics, as well as the lack of standardized protocols and resulting heterogeneity in presenting results. This review aims to provide an overview of the current knowledge and prospects for using viral metagenomics to detect and identify viruses in cattle raised for livestock, while discussing the epidemiological and clinical implications.

Comparison of Three Viral Nucleic Acid Preamplification Pipelines for Sewage Viral Metagenomics

Viral metagenomics is a useful tool for detecting multiple human viruses in urban sewage. However, more refined protocols are required for its effective use in disease surveillance. In this study, we investigated the performance of three different preamplification pipelines (specific to RNA viruses, DNA viruses or both) for viral genome sequencing using spiked-in Phosphate Buffered Saline and sewage samples containing known concentrations of viruses. We found that compared to the pipeline targeting all genome types, the RNA pipeline performed better in detecting RNA viruses in both spiked and unspiked sewage samples, allowing the detection of various mammalian viruses including members from the Reoviridae, Picornaviridae, Astroviridae and Caliciviridae. However, the DNA-specific pipeline did not improve the detection of mammalian DNA viruses. We also measured viral recovery by quantitative reverse transcription polymerase chain reaction and assessed the impact of genetic background (non-viral genetic material) on viral coverage. Our results indicate that viral recoveries were generally lower in sewage (average of 11.0%) and higher in Phosphate Buffered Saline (average of 23.4%) for most viruses. Additionally, spiked-in viruses showed lower genome coverage in sewage, demonstrating the negative effect of genetic background on sequencing. Finally, correlation analysis revealed a relationship between virus concentration and genome normalized reads per million, indicating that viral metagenomic sequencing can be semiquantitative.

Efficient extraction of adventitious virus nucleic acid using commercially available methods

An essential step in pharmaceutical product development is screening for contamination with adventitious agents, and there is desire to develop highly sensitive assays to detect adventitious viral nucleic acid. This study sought to examine the nucleic acid extraction efficiency of three viral candidates in relevant background matrices using four different extraction methods. Three model adventitious viruses, Minute virus of Mice, Porcine Circovirus, and Feline Leukemia Virus, were diluted within a variety of background matrices relevant to pharmaceutical production methods. Upon extraction, the nucleic acid was quantified using droplet digital PCR methods. Four nucleic acid extraction methods were assessed, including commercially available kits and manual extraction methods. Each method recovered nucleic acid post-extraction for each of the model viruses within the tested background matrices. The silica-column based method recovered a greater amount of viral nucleic acid, compared to the other methods tested. Similar trends were observed when model virus was diluted in bioreactor supernatant, which replicates industry testing conditions and provides details on which extraction methods might be used in Next Generation Sequencing and PCR methods for detecting contamination within pharmaceutical products.

Sample and library preparation approaches for the analysis of the virome of irrigation water

BACKGROUND

The virome (i.e. community of mainly RNA and DNA eukaryotic viruses and bacteriophages) of waters is yet to be extensively explored. In particular, the virome of waters used for irrigation could therefore potentially carry viral pathogens that can contaminate fresh produce. One problem in obtaining viral sequences from irrigation waters is the relatively low amount of virus particles, as well as the presence of human, bacterial and protozoan cells. The present aimed study was to compare different processing, amplification, and sequencing approaches for virome characterization in irrigation waters.

RESULTS

Our analyses considered percentages of viral reads, values for diversity indices and number of families found in sequencing results. The results obtained suggest that enrichment protocols using two (bezonase and microccocal nuclease) or four enzymes at once (bezonase, microccocal nuclease, DNAse and RNase), regardless of an Amicon filtration step, are more appropriate than separated enzymatic treatments for virome characterization in irrigation water. The NetoVIR protocol combined with the ScriptSeq v2 RNA-Seq Library (P0-L20 protocol) showed the highest percentages of RNA viruses and identified the higher number of families.

CONCLUSION

Although virome characterization applied in irrigation waters is an important tool for protecting public health by informing on circulating human and zoonotic infections, optimized and standardized procedures should be followed to reduce the variability of results related to either the sample itself and the downstream bioinformatics analyses. Our results show that virome characterization can be an important tool in the discovery of pathogenic viruses in the environment and can be used to inform and optimize reference-based detection methods provided that appropriate and rigorous controls are included. © 2023 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Expanding known viral diversity in plants: virome of 161 species alongside an ancient canal

BACKGROUND

Since viral metagenomic approach was applied to discover plant viruses for the first time in 2006, many plant viruses had been identified from cultivated and non-cultivated plants. These previous researches exposed that the viral communities (virome) of plants have still largely uncharacterized. Here, we investigated the virome in 161 species belonging to 38 plant orders found in a riverside ecosystem.

RESULTS

We identified 245 distinct plant-associated virus genomes (88 DNA and 157 RNA viruses) belonging to 27 known viral families, orders, or unclassified virus groups. Some viral genomes were sufficiently divergent to comprise new species, genera, families, or even orders. Some groups of viruses were detected that currently are only known to infect organisms other than plants. It indicates a wider host range for members of these clades than previously recognized theoretically. We cannot rule out that some viruses could be from plant contaminating organisms, although some methods were taken to get rid of them as much as possible. The same viral species could be found in different plants and co-infections were common.

CONCLUSIONS

Our data describe a complex viral community within a single plant ecosystem and expand our understanding of plant-associated viral diversity and their possible host ranges.

Highly diverse ribonucleic acid viruses in the viromes of eukaryotic host species in Yunnan province, China

Background

The diversity in currently documented viruses and their morphological characteristics indicates the need for understanding the evolutionary characteristics of viruses. Notably, further studies are needed to obtain a comprehensive landscape of virome, the virome of host species in Yunnan province, China.

Materials and methods

We implemented the metagenomic next-generation sequencing strategy to investigate the viral diversity, which involved in 465 specimens collected from bats, pangolins, monkeys, and other species. The diverse RNA viruses were analyzed, especially focusing on the genome organization, genetic divergence and phylogenetic relationships.

Results

In this study, we investigated the viral composition of eight libraries from bats, pangolins, monkeys, and other species, and found several diverse RNA viruses, including the Alphacoronavirus from bat specimens. By characterizing the genome organization, genetic divergence, and phylogenetic relationships, we identified five Alphacoronavirus strains, which shared phylogenetic association with Bat-CoV-HKU8-related strains. The pestivirus-like virus related to recently identified Dongyang pangolin virus (DYPV) strains from dead pangolin specimens, suggesting that these viruses are evolving. Some genomes showed higher divergence from known species (e.g., calicivirus CS9-Cali-YN-CHN-2020), and many showed evidence of recombination events with unknown or known strains (e.g., mamastroviruses BF2-astro-YN-CHN-2020 and EV-A122 AKM5-YN-CHN-2020). The newly identified viruses showed extensive changes and could be assigned as new species, or even genus (e.g., calicivirus CS9-Cali-YN-CHN-2020 and iflavirus Ifla-YN-CHN-2020). Moreover, we identified several highly divergent RNA viruses and estimated their evolutionary characteristics among different hosts, providing data for further examination of their evolutionary dynamics.

Conclusion

Overall, our study emphasizes the close association between emerging viruses and infectious diseases, and the need for more comprehensive surveys.

Comprehensive Evaluation of RNA and DNA Viromic Methods Based on Species Richness and Abundance Analyses Using Marmot Rectal Samples

Viral metagenomics is the most powerful tool to profile viromic composition for a given sample. Different viromic methods, including amplification-free ones, have been developed, but choosing them for different purposes requires comprehensive benchmarks. Here, we assessed the performance of four routinely used methods, i.e., multiple displacement amplification (MDA), direct metagenomic sequencing (MTG), sequence-independent single-primer amplification (SIA), and metatranscriptomic sequencing (MTT), using marmot rectal samples as the templates spiked with five known viruses of different genome types. The obtained clean data were differently contaminated by host and bacterial genomes, resulting in MDA having the most, with ~72.1%, but MTT had only ~7.5% data, useful for follow-up viromic analysis. MDA showed a broader spectrum with higher efficiency to profile the DNA virome, and MTT captured almost all RNA viruses with extraordinary sensitivity; hence, they are advisable in richness-based viromic studies. MTG was weak in capturing single-stranded DNA viruses, and SIA could detect both RNA and DNA viruses but with high randomness. Due to biases to certain types of viruses, the four methods caused different alterations to species abundance compared to the initial virus composition. SIA and MDA introduced greater stochastic errors to relative abundances of species, genus, and family taxa, whereas the two amplification-free methods were more tolerant toward such errors and thus are recommendable in abundance-based analyses. In addition, genus taxon is a compromising analytic level that ensures technically supported and biologically and/or ecologically meaningful viromic conclusions. IMPORTANCE Viral metagenomics can be roughly divided into species richness-based studies and species abundance-based analyses. Viromic methods with different principles have been developed, but rational selection of these techniques according to different purposes requires comprehensive understanding of their properties. By assessing the four most widely used methods using template samples, we found that multiple displacement amplification (MDA) and metatranscriptomic sequencing (MTT) are advisable for species richness-based viromic studies, as they show excellent efficiency to detect DNA and RNA viruses. Meanwhile, metagenomic sequencing (MTG) and MTT are more compatible with stochastic errors of methods introduced into relative abundance of viromic taxa and hence are rational choices in species abundance-based analyses. This study also highlights that MTG needs to tackle host genome contamination and ameliorate the capacity to detect single-stranded DNA viruses in the future, and the MTT method requires an improvement in bacterial rRNA depletion prior to library preparation.

Contribution of single-cell omics to microbial ecology

Micro-organisms play key roles in various ecosystems, but many of their functions and interactions remain undefined. To investigate the ecological relevance of microbial communities, new molecular tools are being developed. Among them, single-cell omics assessing genetic diversity at the population and community levels and linking each individual cell to its functions is gaining interest in microbial ecology. By giving access to a wider range of ecological scales (from individual to community) than culture-based approaches and meta-omics, single-cell omics can contribute not only to micro-organisms' genomic and functional identification but also to the testing of concepts in ecology. Here, we discuss the contribution of single-cell omics to possible breakthroughs in concepts and knowledge on microbial ecosystems and ecoevolutionary processes.

The Diagnostic Value of Metagenomic Next-Generation Sequencing in Lower Respiratory Tract Infection

Lower respiratory tract infections are associated with high morbidity and mortality and significant clinical harm. Due to the limited ability of traditional pathogen detection methods, anti-infective therapy is mostly empirical. Therefore, it is difficult to adopt targeted drug therapy. In recent years, metagenomic next-generation sequencing (mNGS) technology has provided a promising means for pathogen-specific diagnosis and updated the diagnostic strategy for lower respiratory tract infections. This article reviews the diagnostic value of mNGS for lower respiratory tract infections, the impact of different sampling methods on the detection efficiency of mNGS, and current technical difficulties in the clinical application of mNGS.

New Parvoviruses and Picornavirus in Tissues and Feces of Foals with Interstitial Pneumonia

Six foals with interstitial pneumonia of undetermined etiology from Southern California were analyzed by viral metagenomics. Spleen, lung, and colon content samples obtained during necropsy from each animal were pooled, and nucleic acids from virus-like particles enriched for deep sequencing. The recently described equine copiparvovirus named eqcopivirus, as well as three previously uncharacterized viruses, were identified. The complete ORFs genomes of two closely related protoparvoviruses, and of a bocaparvovirus, plus the partial genome of a picornavirus were assembled. The parvoviruses were classified as members of new ungulate protoparvovirus and bocaparvovirus species in the Parvoviridae family. The picornavirus was classified as a new species in the Salivirus genus of the Picornaviridae family. Spleen, lung, and colon content samples from each foal were then tested for these viral genomes by nested PCR and RT-PCR. When present, parvoviruses were detected in both feces and spleen. The picornavirus, protoparvovirus, and eqcopivirus genomes were detected in the lungs of one animal each. Three foals were co-infected with the picornavirus and either a protoparvovirus, bocaparvovirus, or eqcopivirus. Two other foals were infected with a protoparvovirus only. No viral infection was detected in one animal. The complete ORFs of the first equine protoparvoviruses and bocaparvovirus, the partial ORF of the third equine picornavirus, and their detection in tissues of foals with interstitial pneumonia are described here. Testing the involvement of these viruses in fatal interstitial pneumonia or other equine diseases will require larger epidemiological and/or inoculation studies.

Ultrasensitive version of nucleic acid sequence-based amplification (NASBA) utilizing a nicking and extension chain reaction system

Nucleic acid sequence-based amplification (NASBA) is a transcription-based isothermal amplification technique especially designed for the detection of RNA targets. The NASBA basically relies on the linear production of T7 RNA promoter-containing double-stranded DNA (T7DNA), and thus the final amplification efficiency is not sufficiently high enough to achieve ultrasensitive detection. We herein ingeniously integrate a nicking and extension chain reaction system into the NASBA to establish an ultrasensitive version of NASBA, termed Nicking and Extension chain reaction System-Based Amplification (NESBA). By employing a NESBA primer set designed to contain an additional nicking site at the 5' end of a NASBA primer set, the T7DNA is exponentially amplified through continuously repeated nicking and extension chain reaction by the combined activities of nicking endonuclease (NE) and reverse transcriptase (RT). As a consequence, a much larger number of RNA amplicons would be produced through the transcription of the amplified T7DNA, greatly enhancing the final fluorescence signal from the molecular beacon (MB) probe binding to the RNA amplicon. Based on this unique design principle, we successfully identified the target respiratory syncytial virus A (RSV A) genomic RNA (gRNA) down to 1 aM under isothermal conditions, which is 100-fold more sensitive than regular NASBA.

Circovirus in Blood of a Febrile Horse with Hepatitis

Circoviruses infect vertebrates where they can result in a wide range of disease signs or in asymptomatic infections. Using viral metagenomics we analyzed a pool of five sera from four healthy and one sick horse. Sequences from parvovirus-H, equus anellovirus, and distantly related to mammalian circoviruses were recognized. PCR identified the circovirus reads as originating from a pregnant mare with fever and hepatitis. That horse's serum was also positive by real time PCR for equine parvovirus H and negative for the flavivirus equine hepacivirus. The complete circular genome of equine circovirus 1 strain Charaf (EqCV1-Charaf) was completed using PCR and Sanger sequencing. EqCV1 replicase showed 73-74% identity to those of their closest relatives, pig circoviruses 1/2, and elk circovirus. The closest capsid proteins were from the same ungulate circoviruses with 62-63% identity. The overall nucleotide identity of 72% to its closest relative indicates that EqCV1 is a new species in the Circovirus genus, the first reported in genus Equus. Whether EqCV1 alone or in co-infections can result in disease and its prevalence in different equine populations will require further studies now facilitated using EqCV1's genome sequence.

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.

Norovirus strains in patients with acute gastroenteritis in rural and low-income urban areas in northern Brazil

From 2010-2016, a total of 251 stool samples were screened for norovirus using next-generation sequencing (NGS) followed by phylogenetic analysis to investigate the genotypic diversity of noroviruses in rural and low-income urban areas in northern Brazil. Norovirus infection was detected in 19.9% (50/251) of the samples. Eight different genotypes were identified: GII.4_Sydney[P31] (64%, 32/50), GII.6[P7] (14%, 7/50), GII.17[P17] (6%, 3/50), GII.1[P33] (6%, 3/50), GII.3[P16] (4%, 2/50), GII.2[P16] (2%, 1/50), GII.2[P2] (2%, 1/50), and GII.4_New Orleans[P4] (2%, 1/50). Distinct GII.6[P7] variants were recognized, indicating the presence of different co-circulating strains. Elucidating norovirus genetic diversity will improve our understanding of their potential health burden, in particular for the GII.4_Sydney[P31] variant.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Astrovirus Outbreak in an Animal Shelter Associated With Feline Vomiting

An outbreak of cat vomiting was observed in an animal shelter. Testing for known enteric feline pathogens did not identify a causative agent. Viral metagenomics on four mini pools of feces from cases and controls housed in the same area revealed the presence of feline astrovirus in all pools. Also found with fewer reads in one pool each were rotavirus I, carnivore bocaparvovirus 3, norovirus (NoV) GVI, and a novel dependovirus. The genome of the highly prevalent astrovirus was sequenced and classified into mamastrovirus species two, also known as feline astrovirus. Real-time RT-PCR on longitudinally acquired fecal samples from 11 sick cases showed 10 (91%) to be shedding astrovirus for as long as 19 days. Affected cats were sick for an average of 9.8 days, with a median of 2.5 days (range = 1–31 days). Unaffected control cats housed in the same areas during the outbreak showed five out of nine (56%) to also be shedding astrovirus. Feline fecal samples collected from the same animal shelter ~1 year before (n = 8) and after (n = 10) showed none to be shedding astrovirus, indicating that this virus was temporarily associated with the vomiting outbreak and is not part of the commensal virome for cats in this shelter. Together with the absence of highly prevalent known pathogens, our results support a role for feline astrovirus infection, as well as significant asymptomatic shedding, in an outbreak of contagious feline vomiting.

Deep Impact of Random Amplification and Library Construction Methods on Viral Metagenomics Results

Clinical metagenomics is a broad-range agnostic detection method of pathogens, including novel microorganisms. A major limit is the low pathogen load compared to the high background of host nucleic acids. To overcome this issue, several solutions exist, such as applying a very high depth of sequencing, or performing a relative enrichment of viral genomes associated with capsids. At the end, the quantity of total nucleic acids is often below the concentrations recommended by the manufacturers of library kits, which necessitates to random amplify nucleic acids. Using a pool of 26 viruses representative of viral diversity, we observed a deep impact of the nature of sample (total nucleic acids versus RNA only), the reverse transcription, the random amplification and library construction method on virus recovery. We further optimized the two most promising methods and assessed their performance with fully characterized reference virus stocks. Good genome coverage and limit of detection lower than 100 or 1000 genome copies per mL of plasma, depending on the genome viral type, were obtained from a three million reads dataset. Our study reveals that optimized random amplification is a technique of choice when insufficient amounts of nucleic acid are available for direct libraries constructions.

Virome of Bat Guano from Nine Northern California Roosts

Bats are hosts to a large variety of viruses, including many capable of cross-species transmissions to other mammals, including humans. We characterized the virome in guano from five common bat species in 9 Northern California roosts and from a pool of 5 individual bats. Genomes belonging to 14 viral families known to infect mammals and 17 viral families infecting insects or of unknown tropism were detected. Nearly complete or complete genomes of a novel parvovirus, astrovirus, nodavirus, circular Rep-encoding single-stranded DNA (CRESS-DNA) viruses, and densoviruses, and more partial genomes of a novel alphacoronavirus and a bunyavirus were characterized. Lower numbers of reads with >90% amino acid identity to previously described calicivirus, circovirus, adenoviruses, hepatovirus, bocaparvoviruses, and polyomavirus in other bat species were also found, likely reflecting their wide distribution among different bats. Unexpectedly, a few sequence reads of canine parvovirus 2 and the recently described mouse kidney parvovirus were also detected and their presence confirmed by PCR; these possibly originated from guano contamination by carnivores and rodents. The majority of eukaryotic viral reads were highly divergent, indicating that numerous viruses still remain to be characterized, even from such a heavily investigated order as Chiroptera.IMPORTANCE Characterizing the bat virome is important for understanding viral diversity and detecting viral spillover between animal species. Using an unbiased metagenomics method, we characterize the virome in guano collected from multiple roosts of common Northern California bat species. We describe several novel viral genomes and report the detection of viruses with close relatives reported in other bat species, likely reflecting cross-species transmissions. Viral sequences from well-known carnivore and rodent parvoviruses were also detected, whose presence are likely the result of contamination from defecation and urination atop guano and which reflect the close interaction of these mammals in the wild.

Implementation of next-generation sequencing for virus identification in veterinary diagnostic laboratories

The value of next-generation sequencing (NGS)-based applications for testing purposes in human medicine is widely recognized. Although NGS-based metagenomic screening may be of interest in veterinary medicine, in particular for intensively farmed livestock species such as pigs, there is a lack of protocols tailored to veterinary requirements, likely because of the high diversity of species and samples. Therefore, we developed an NGS-based protocol for use in veterinary virology and present here different applications in porcine medicine. To develop the protocol, each step of sample preparation was optimized using porcine samples spiked with various RNA and DNA viruses. The resulting protocol was tested with clinical samples previously confirmed to be positive for specific viruses by a diagnostic laboratory. Additionally, we validated the protocol in an NGS viral metagenomics ring trial and tested the protocol on viral multiplex reference material (NIBSC, U.K.). We applied our ViroScreen protocol successfully for 1) virus identification, 2) virus characterization, and 3) herd screening. We identified torque teno sus virus and atypical porcine pestivirus in a neurologic case, determined the full-length genome sequence of swine influenza A virus in field samples, and screened pigs using pen floor fecal samples and chewing rope liquid.

Comparison of Nucleic Acid Extraction Methods for a Viral Metagenomics Analysis of Respiratory Viruses

Viral metagenomics next-generation sequencing (mNGS) is increasingly being used to characterize the human virome. The impact of viral nucleic extraction on virome profiling has been poorly studied. Here, we aimed to compare the sensitivity and sample and reagent contamination of three extraction methods used for viral mNGS: two automated platforms (eMAG; MagNA Pure 24, MP24) and the manual QIAamp Viral RNA Mini Kit (QIAamp). Clinical respiratory samples (positive for Respiratory Syncytial Virus or Herpes Simplex Virus), one mock sample (including five viruses isolated from respiratory samples), and a no-template control (NTC) were extracted and processed through an mNGS workflow. QIAamp yielded a lower proportion of viral reads for both clinical and mock samples. The sample cross-contamination was higher when using MP24, with up to 36.09% of the viral reads mapping to mock viruses in the NTC (vs. 1.53% and 1.45% for eMAG and QIAamp, respectively). The highest number of viral reads mapping to bacteriophages in the NTC was found with QIAamp, suggesting reagent contamination. Our results highlight the importance of the extraction method choice for accurate virome characterization.

Serendipitous Discovery of a Novel Murine Astrovirus Contaminating a Murine Helper T-cell Line and Incapable of Infecting Highly Immunodeficient Mice

The unexpected seroconversion of sentinel mice in our facility to murine T lymphotrophic virus (MTLV) positivity led to our identification of a novel murine astrovirus that we designated murine astrovirus 2 (MuAstV-2). During our investigation, MuAstV-2 was found to be a contaminant of the T helper cell line (D10. G4.1) that was used to generate the MTLV antigen that we included in the multiplex fluorometric immunoassay (MFIA) that we used for sentinel screening. We eventually determined that cross-reactivity with the astrovirus generated a positive result in the MTLV assay. A confirmatory immunofluorometric assay (IFA) using the same MTLV-infected cell line yielded a similar result. However, the use of antigen prepared from MTLV-infected neonatal mouse thymus did not reproduce a positive result, leading us to suspect that the seroreactivity we had observed was not due to infection with MTLV. A mouse antibody production test showed that mice inoculated with naïve D10. G4.1 cells and their contact sentinels tested positive for MTLV using cell-line generated antigen, but tested negative in assays using MTLV antigen produced in mice. Metagenomic analysis was subsequently used to identify MuAstV-2 in feces from 2 sentinel mice that had recently seroconverted to MTLV. Two closely related astrovirus sequences (99.6% capsid identity) were obtained and shared 95% capsid amino acid identity with the MuAstV-2 virus sequenced from the D10. G4.1 cell line. These viruses are highly divergent from previously identified murine astroviruses, displaying <30% capsid identity, yet were closely related to murine astrovirus 2 (85% capsid identity), which had recently been isolated from feral mice in New York City. A MuAstV-2 specific PCR assay was developed and used to eradicate MuAstV-2 from the infected colony using a test and cull strategy. The newly identified MuAstV2 readily transmits to immunocompetent mouse strains by fecal-oral exposure, but fails to infect NOD-Prkdcem26Cd52Il2rgem26Cd22/NjuCrl (NCG) mice, which have significantly impaired adaptive and innate immune systems. Neither immunocompetent nor immunodeficient mice showed any astrovirus-associated pathology. MuAstV-2 may provide a valuable model for the study of specific aspects of astrovirus pathogenesis and virus-host interactions.

Complex Virome in a Mesenteric Lymph Node from a Californian Sea Lion (Zalophus Californianus) with Polyserositis and Steatitis

An emaciated subadult free-ranging California sea lion (Csl or Zalophus californianus) died following stranding with lesions similar to 11 other stranded animals characterized by chronic disseminated granulomatous inflammation with necrotizing steatitis and vasculitis, involving visceral adipose tissues in the thoracic and peritoneal cavities. Histologically, affected tissues had extensive accumulations of macrophages with perivascular lymphocytes, plasma cells, and fewer neutrophils. Using viral metagenomics on a mesenteric lymph node six mammalian viruses were identified consisting of novel parvovirus, polyomavirus, rotavirus, anellovirus, and previously described Csl adenovirus 1 and Csl bocavirus 4. The causal or contributory role of these viruses to the gross and histologic lesions of this sea lion remains to be determined.

The virome in hematology-Stem cell transplantation and beyond

The virome has been recently studied in hematology and mostly in the setting of allogeneic hematopoietic stem cell transplantation. However, in hematology (as in the setting of nonhematological disorders) the study of the microbiome (that indeed includes the virome) is a growing field. The overall field is moving beyond species catalogue to the understanding of the complex ecological relationship that microbes have with each other and with their host. Here we review the existing literature on the virome in transplant recipients and in other settings, and discuss potential applications of the virome study in hematology.

Virome of a Feline Outbreak of Diarrhea and Vomiting Includes Bocaviruses and a Novel Chapparvovirus

An unexplained outbreak of feline diarrhea and vomiting, negative for common enteric viral and bacterial pathogens, was subjected to viral metagenomics and PCR. We characterized from fecal samples the genome of a novel chapparvovirus we named fechavirus that was shed by 8/17 affected cats and identified three different feline bocaviruses shed by 9/17 cats. Also detected were nucleic acids from attenuated vaccine viruses, members of the normal feline virome, viruses found in only one or two cases, and viruses likely derived from ingested food products. Epidemiological investigation of disease signs, time of onset, and transfers of affected cats between three facilities support a possible role for this new chapparvovirus in a highly contagious feline diarrhea and vomiting disease.

Semen virome of men with HIV on or off antiretroviral treatment

OBJECTIVES

Improving immune status of people living with HIV through antiretroviral therapy (ART) may also reduce shedding of other viruses in semen. We characterized the seminal fluid virome of men with HIV and tested potential associations between viruses present and CD4 T-cell count, HIV viremia, and antiretroviral therapy (ART) status.

DESIGN AND METHODS

Metagenomics was used to enrich and sequence viral nucleic acids from the seminal fluid of 55 semen samples from 42 men living with HIV from San Francisco with a median age of 33 (IQR, 28.7-45) and median CD4 T-cell counts of 837 cells/μl (IQR, 258-1571 cells/μl). All samples were collected between 2005 and 2015, and ART status was ascertained from medical records.

RESULTS

Anelloviruses, cytomegalovirus (CMV), and multiple genotypes of human papillomaviruses were detected. Participants shed from 0 to 4 distinct human viruses. Longitudinally collected seminal fluid samples showed changes in the viruses shed. Viruses were more frequently shed by individuals with detectable HIV viremia (43.7 vs. 15.4%, P = 0.042). A trend was seen for increased shedding by individuals who were not on ART (42.8 vs. 17.8%, P = 0.082) or with CD4 T-cell count less than 350 cells/μl (35.3 vs. 20%, P = 0.27).

CONCLUSION

Seminal fluid from men with HIV from San Francisco contains nucleic acids from three different DNA viral families. A greater number of viruses, particularly CMV, were shed by participants with detectable HIV viremia (18.9 vs. 0%, P = 0.022). Control of viremia through ART may lower shedding of other viruses in semen in addition to HIV.

Detection and characterization of Ilheus and Iguape virus genomes in historical mosquito samples from Southern Brazil

In Brazil, flaviviruses have caused massive outbreaks. Surveillance programs designed to monitor virus activity in vectors provides a system for mapping disease distribution and for identifying specific vector species for targeted control. The present study aimed to describe the detection, whole genome characterization and phylogenetic analysis of Ilheus virus (ILHV) and Iguape virus (IGUV) strains obtained from historical mosquito's samples. Twelve isolates of pooled mosquito specimens (inoculated in neonate mouse brain) collected in the state of São Paulo, Brazil, in 1993, 1994 and 1997 were investigated. Viral RNA was extracted and analyzed by qRT-PCR using Flavivirus genus-specific primers. Positive samples were sequenced and underwent phylogenetic analyses. Flavivirus was detected in 50% of the specimens. Positive samples were successfully Sanger sequenced. Three Anopholes cruzii pools collected in 1994 were positive for IGUV. One Culex sp. pool, one Anopheles triannulatus pool, and one Coquillettidia juxtamansonia pool, collected in 1994, were positive for ILHV. Metagenomic sequencing successfully characterize one ILHV and four IGUV full genomes, and revealed a high degree of homology between the Brazilian ILHV and IGUV strains and isolates available in GenBank. Phylogenetic analysis of partial ILHV NS5 gene revealed three distinct lineages (clades), an indication of genetic heterogeneity in strains circulating in Brazil. Nucleotide insertions and a high-level of nucleotide diversity were observed in the NS1 protein and capsid region of IGUV strains, respectively. Detection of ILHV and IGUV in mosquitoes from Southeastern Brazil confirms the historical circulation of these viruses in this area. Furthermore, this first evidence of ILHV in Anopheles triannulatus suggests the potential importance of Anopheles mosquitoes in the IGUV transmission cycle. Genomic and phylogenetic analysis of these viruses provided insights into their diversity and evolution, which are important for the emergence patterns of flaviviruses and their evolutionary trends in Brazil, an endemic country for several arbovirus. in In-depth studies of ILHV and IGUV including vector competence and molecular studies are needed to shed light on their epidemiology and potential risk of future emergence.

Next-generation sequencing of dsRNA is greatly improved by treatment with the inexpensive denaturing reagent DMSO

dsRNA is the genetic material of important viruses and a key component of RNA interference-based immunity in eukaryotes. Previous studies have noted difficulties in determining the sequence of dsRNA molecules that have affected studies of immune function and estimates of viral diversity in nature. DMSO has been used to denature dsRNA prior to the reverse-transcription stage to improve reverse transcriptase PCR and Sanger sequencing. We systematically tested the utility of DMSO to improve the sequencing yield of a dsRNA virus (Φ6) in a short-read next-generation sequencing platform. DMSO treatment improved sequencing read recovery by over two orders of magnitude, even when RNA and cDNA concentrations were below the limit of detection. We also tested the effects of DMSO on a mock eukaryotic viral community and found that dsRNA virus reads increased with DMSO treatment. Furthermore, we provide evidence that DMSO treatment does not adversely affect recovery of reads from a ssRNA viral genome (influenza A/California/07/2009). We suggest that up to 50 % DMSO treatment be used prior to cDNA synthesis when samples of interest are composed of or may contain dsRNA.

Metatranscriptomic virome assessment of Rhipicephalus microplus from Colombia

Ticks (Ixodida) are hematophagous ectoparasites that harbor and transmit diverse species of viruses, some of which cause serious diseases with worldwide veterinary and human health consequences. Rhipicephalus microplus is an important cattle tick in Colombia, where it causes significant economic losses. Despite the importance of this tick, its viral profile is unknown. RNA sequencing was used in this study as a surveillance method for virus detection in R. microplus. Most of the viral origin contigs were assigned to two putative viruses: one chuvirus (Wuhan tick virus 2) and one phlebovirus-like (Lihan tick virus). In addition, viral contigs corresponding to two jingmenviruses previously reported in R. microplus from China and Brazil were detected, as well as a novel putative tymovirus, named here as Antioquia tymovirus-like 1 (ATV-like 1). The presence of some of these viruses across numerous regions in the world could have several explanations, including i) a long-term association between those viruses and R. microplus and ii) a consequence of livestock historical trade. Our results shed new light on the virus diversity of this tick species and provide a basis for further studies on the evolutionary history and pathogenic potential of these interesting viruses.

Plasma virome of 781 Brazilians with unexplained symptoms of arbovirus infection include a novel parvovirus and densovirus

Plasma from patients with dengue-like symptoms was collected in 2013 to 2016 from the Brazilian states of Tocantins and Amapa. 781 samples testing negative for IgM against Dengue, Zika, and Chikungunya viruses and for flaviviruses, alphaviruses and enteroviruses RNA using RT-PCRs were analyzed using viral metagenomics. Viral particles-associated nucleic acids were enriched, randomly amplified, and deep sequenced in 102 mini-pools generating over 2 billion reads. Sequence data was analyzed for the presence of known and novel eukaryotic viral reads. Anelloviruses were detected in 80%, human pegivirus 1 in 19%, and parvovirus B19 in 17% of plasma pools. HIV and enteroviruses were detected in two pools each. Previously uncharacterized viral genomes were also identified, and their presence in single plasma samples confirmed by PCR. Chapparvovirus and ambidensovirus genomes, both in the Parvoviridae family, were partially characterized showing 33% and 34% identity in their NS1 sequences to their closest relative. Molecular surveillance using pre-existing plasma from febrile patients provides a readily scalable approach for the detection of novel, potentially emerging, viruses.

Nicotinamide activates latent HIV-1 ex vivo in ART suppressed individuals, revealing higher potency than the association of two methyltransferase inhibitors, chaetocin and BIX01294

Background

Latent HIV-1 is a major hurdle in obtaining HIV-1 sustained virological remission (SVR). Here we explored histone deacetylation inhibition property of nicotinamide (NAM; n = 17) for the first time in comparison to a combination of methyltransferase inhibitors (MTIs; Chaetocin and BIX01294; n = 25) to reactivate latent HIV ex vivo in CD8-depleted PBMCs from antiretroviral treated aviremic individuals.

Results

NAM reactivated HIV-1 from 13/17 (76.4%) samples compared to 20/25 (80.0%) using MTIs with mean viral load (VLs) of 4.32 and 3.22 log₁₀ RNA copies/mL, respectively (p = 0.004). Mean purging time after NAM and MTIs stimulation was 5.1 and 6.75 days, respectively (p = 0.73). Viral purging in autologous cultures exhibited blunted HIV recovery with fluctuating VLs followed by a complete viral extinction when expanded in allogenic system. Electron microscopy from five supernatants revealed anomalous viral particles, with lack of complete viral genomes when characterized by ultradeep sequencing through metagenomics approach (n = 4).

Conclusion

NAM alone was more potent HIV-1 activator than combination of MTIs, with potential of clinical use.

Performance of Metagenomic Next-Generation Sequencing for the Diagnosis of Viral Meningoencephalitis in a Resource-Limited Setting

Background

Meningoencephalitis is a devastating disease worldwide. Current diagnosis fails to establish the cause in ≥50% of patients. Metagenomic next-generation sequencing (mNGS) has emerged as pan-pathogen assays for infectious diseases diagnosis, but few studies have been conducted in resource-limited settings.

Methods

We assessed the performance of mNGS in the cerebrospinal fluid (CSF) of 66 consecutively treated adults with meningoencephalitis in a tertiary referral hospital for infectious diseases in Vietnam, a resource-limited setting. All mNGS results were confirmed by viral-specific polymerase chain reaction (PCR). As a complementary analysis, 6 viral PCR-positive samples were analyzed using MinION-based metagenomics.

Results

Routine diagnosis could identify a virus in 15 (22.7%) patients, including herpes simplex virus (HSV; n = 7) and varicella zoster virus (VZV; n = 1) by PCR, and mumps virus (n = 4), dengue virus (DENV; n = 2), and Japanese encephalitis virus (JEV; n = 1) by serological diagnosis. mNGS detected HSV, VZV, and mumps virus in 5/7, 1/1, and 1/4 of the CSF positive by routine assays, respectively, but it detected DENV and JEV in none of the positive CSF. Additionally, mNGS detected enteroviruses in 7 patients of unknown cause. Metagenomic MinION-Nanopore sequencing could detect a virus in 5/6 PCR-positive CSF samples, including HSV in 1 CSF sample that was negative by mNGS, suggesting that the sensitivity of MinION is comparable with that of mNGS/PCR.

Conclusions

In a single assay, metagenomics could accurately detect a wide spectrum of neurotropic viruses in the CSF of meningoencephalitis patients. Further studies are needed to determine the value that real-time sequencing may contribute to the diagnosis and management of meningoencephalitis patients, especially in resource-limited settings where pathogen-specific assays are limited in number.

Detection and complete genome characterization of a novel RNA virus related to members of the Hepe-Virga clade in bird species, hoopoe (Upupa epops)

Using viral metagenomics, next-generation sequencing and RT-PCR techniques a genetically divergent hepevirus-like RNA virus was identified and characterized from a faecal sample of wild bird species, hoopoe (Upupa epops) in Hungary. The complete viral genome sequence of hoopoe/BBanka01/2015/HUN (GenBank accession number MN852439) is 7052 nt long including a 54-nt 5' and an 18-nt 3' non-coding region without poly(A)-tail. Sequence analysis indicated that the hoopoe/BBanka01/2015/HUN genome has potentially three overlapping open reading frames (ORFs). The ORF1 (6558 nt/2185aa) encodes a long, non-structural polyprotein (replicase) including putative functional domains and conserved aa motifs of methyltransferase with domain Y, RNA helicase and RdRp and has <33% aa identity to the known hepe- and hepe-like viruses. The ORF2 (1446 nt/481aa) encodes a putative structural (capsid) protein overlapping with ORF1 but translated in different coding frame. The functions of the short ORF3 (426 nt/141aa) were not predictable. Similar virus sequences were not detected from samples from 21 further bird species. The taxonomic position of this novel virus is presently unknown.

Metagenomics revealing new virus species in farm and pet animals and aquaculture

Viral metagenomics is slowly taking over the traditional and widely used molecular techniques for the investigation of pathogenic viruses responsible for illness and inflicting great economic burden on the farm animal industry. Owing to the continued improvements in sequencing technologies and the dramatic reduction of per base costs of sequencing the use of next generation sequencing have been key factors in this progress. Discoveries linked to viral metagenomics are expected to be beneficial to the field of veterinary medicine starting from the development of better diagnostic assays to the design of new subunit vaccines with minimal investments. With these achievements the research has taken a giant leap even toward the better healthcare of animals and, as a result, the animal sector could be growing at an unprecedented pace.

Sensitivity and specificity of metatranscriptomics as an arbovirus surveillance tool

The ability to identify all the viruses within a sample makes metatranscriptomic sequencing an attractive tool to screen mosquitoes for arboviruses. Practical application of this technique, however, requires a clear understanding of its analytical sensitivity and specificity. To assess this, five dilutions (1:1, 1:20, 1:400, 1:8,000 and 1:160,000) of Ross River virus (RRV) and Umatilla virus (UMAV) isolates were spiked into subsamples of a pool of 100 Culex australicus mosquitoes. The 1:1 dilution represented the viral load of one RRV-infected mosquito in a pool of 100 mosquitoes. The subsamples underwent nucleic acid extraction, mosquito-specific ribosomal RNA depletion, and Illumina HiSeq sequencing. The viral load of the subsamples was also measured using reverse transcription droplet digital PCR (RT-ddPCR) and quantitative PCR (RT-qPCR). Metatranscriptomic sequencing detected both RRV and UMAV in the 1:1, 1:20 and 1:400 subsamples. A high specificity was achieved, with 100% of RRV and 99.6% of UMAV assembled contigs correctly identified. Metatranscriptomic sequencing was not as sensitive as RT-qPCR or RT-ddPCR; however, it recovered whole genome information and detected 19 other viruses, including four first detections for Australia. These findings will assist arbovirus surveillance programs in utilising metatranscriptomics in routine surveillance activities to enhance arbovirus detection.

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.

Viruses in Horses with Neurologic and Respiratory Diseases

Metagenomics was used to identify viral sequences in the plasma and CSF (cerobrospinal fluid) of 13 horses with unexplained neurological signs and in the plasma and respiratory swabs of 14 horses with unexplained respiratory signs. Equine hepacivirus and two copiparvoviruses (horse parvovirus-CSF and a novel parvovirus) were detected in plasma from neurological cases. Plasma from horses with respiratory signs contained the same two copiparvoviruses plus equine pegivirus D and respiratory swabs contained equine herpes virus 2 and 5. Based on genetic distances the novel copiparvovirus qualified as a member of a new parvovirus species we named Eqcopivirus. These samples plus another 41 plasma samples from healthy horses were tested by real-time PCRs for multiple equine parvoviruses and hepacivirus. Over half the samples tested were positive for one to three viruses with eqcopivirus DNA detected in 20.5%, equine hepacivirus RNA and equine parvovirus-H DNA in 16% each, and horse parvovirus-CSF DNA in 12% of horses. Comparing viral prevalence in plasma none of the now three genetically characterized equine parvoviruses (all in the copiparvovirus genus) was significantly associated with neurological and respiratory signs in this limited sampling.

Effect of Geographic Isolation on the Nasal Virome of Indigenous Children

The influence of living in small remote villages on the diversity of viruses in the nasal mucosa was investigated in three Colombian villages with very different levels of geographic isolation. Metagenomic analysis was used to characterize viral nucleic acids in nasal swabs from 63 apparently healthy young children. Sequences from human virus members of the families Anelloviridae, Papillomaviridae, Picornaviridae, Herpesviridae, Polyomaviridae, Adenoviridae, and Paramyxoviridae were detected in decreasing proportions of children. The number of papillomavirus infections detected was greater among Hispanic children most exposed to outside contacts, while anellovirus infections were more common in the isolated indigenous villages. The diversity of the other human viruses detected did not differ among the villages. Closely related variants of rhinovirus A or B were identified in 2 to 4 children from each village, reflecting ongoing transmission clusters. Genomes of viruses not currently known to infect humans, including members of the families Parvoviridae, Partitiviridae, Dicistroviridae, and Iflaviridae and circular Rep-encoding single-stranded DNA (CRESS-DNA) virus, were also detected in nasal swabs, possibly reflecting environmental contamination from insect, fungal, or unknown sources. Despite the high levels of geographic and cultural isolation, the overall diversity of human viruses in the nasal passages of children was not reduced in highly isolated indigenous villages, indicating ongoing exposure to globally circulating viruses.IMPORTANCE Extreme geographic and cultural isolation can still be found in some indigenous South American villages. Such isolation may be expected to limit the introduction of otherwise common and widely distributed viruses. Very small population sizes may also result in rapid local viral extinction due to a lack of seronegative subjects to maintain transmission chains for rapidly cleared viruses. We compared the viruses in the nasal passages of young children in three villages with increasing levels of geographic isolation. We found that isolation did not reduce the overall diversity of viral infections. Multiple infections with nearly identical rhinoviruses could be detected within each village, likely reflecting recent viral introductions and transmission clusters among epidemiologically linked members of these very small communities. We conclude that, despite their geographic isolation, remote indigenous villages show evidence of ongoing exposure to globally circulating viruses.

The rare enterovirus c99 and echovirus 29 strains in Brazil: potential risks associated to silent circulation

Human enteroviruses (EVs) are associated with a wide spectrum of human diseases. Here we report the complete genome sequences of one EV-C99 strain and one E29 strain obtained from children suffering from acute gastroenteritis, without symptoms of enteroviral syndromes. This is the first report of EV-C99 in South America, and the second E29 genome described worldwide. Continuous surveillance on EVs is vital to provide further understanding of the circulation of new or rare EV serotypes in the country. The present study also highlights the capacity of EVs to remain in silent circulation in populations.

Lyon-IARC Polyomavirus DNA in Feces of Diarrheic Cats

A viral metagenomic analysis of feces from an unexplained outbreak of feline diarrhea revealed the presence of Lyon-IARC polyomavirus (LIPyV) DNA. LIPyV, whose genome was originally sequenced from swabs of human skin, was fecally shed by three out of five diarrheic cats.

A viral metagenomic analysis of feces from an unexplained outbreak of feline diarrhea revealed the presence of Lyon-IARC polyomavirus (LIPyV) DNA. LIPyV, whose genome was originally sequenced from swabs of human skin, was fecally shed by three out of five diarrheic cats.

Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics

The Tasmanian devil is an endangered carnivorous marsupial threatened by devil facial tumor disease (DFTD). While research on DFTD has been extensive, little is known about viruses in devils and whether any are of potential conservation relevance for this endangered species. Using both metagenomics based on virion enrichment and sequence-independent amplification (virion-enriched metagenomics) and metatranscriptomics based on bulk RNA sequencing, we characterized and compared the fecal viromes of captive and wild devils. A total of 54 fecal samples collected from two captive and four wild populations were processed for virome characterization using both approaches. In total, 24 novel marsupial-related viruses, comprising a sapelovirus, astroviruses, rotaviruses, picobirnaviruses, parvoviruses, papillomaviruses, polyomaviruses, and a gammaherpesvirus, were identified, as well as known mammalian pathogens such as rabbit hemorrhagic disease virus 2. Captive devils showed significantly lower viral diversity than wild devils. Comparison of the two virus discovery approaches revealed substantial differences in the number and types of viruses detected, with metatranscriptomics better suited for RNA viruses and virion-enriched metagenomics largely identifying more DNA viruses. Thus, the viral communities revealed by virion-enriched metagenomics and metatranscriptomics were not interchangeable and neither approach was able to detect all viruses present. An integrated approach using both virion-enriched metagenomics and metatranscriptomics constitutes a powerful tool for obtaining a complete overview of both the taxonomic and functional profiles of viral communities within a sample.IMPORTANCE The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

Recombinant Strains of Human Parechovirus in Rural Areas in the North of Brazil

We characterized the 24 nearly full-length genomes of human parechoviruses (PeV) from children in the north of Brazil. The initial phylogenetic analysis indicated that 17 strains belonged to genotype 1, 5 to genotype 4, and 1 to genotype 17. A more detailed analysis revealed a high frequency of recombinant strains (58%): A total of 14 of our PeV-As were chimeric, with four distinct recombination patterns identified. Five strains were composed of genotypes 1 and 5 (Rec1/5); five strains shared a complex mosaic pattern formed by genotypes 4, 5, and 17 (Rec4/17/5); two strains were composed of genotypes 1 and 17 (Rec1/17); and two strains were composed of genotype 1 and an undetermined strain (Rec1/und). Coalescent analysis based on the Vp1 gene, which is free of recombination, indicated that the recombinant strains most likely arose in this region approximately 30 years ago. They are present in high frequencies and are circulating in different small and isolated cities in the state of Tocantins. Further studies will be needed to establish whether the detected recombinant strains have been replacing parental strains or if they are co-circulating in distinct frequencies in Tocantins.

The first reptilian circovirus identified infects gut and liver tissues of black-headed pythons

Viral metagenomic analysis of the liver of a black headed python (Aspidites melanocephalus) euthanized for a proliferative spinal lesion of unknown etiology yielded the first characterized genome of a reptile-infecting circovirus (black-headed python circovirus or BhPyCV). BhPyCV-specific in situ hybridization (ISH) showed that viral nucleic acids were strongly expressed in the intestinal lining and mucosa and multifocally in the liver. To investigate the presence of this virus in other snakes and its possible pathogenicity, 17 snakes in the python family with spinal disease were screened with ISH yielding a second BhP positive in intestinal tissue, and a Boelen’s python (Morelia boeleni) positive in the liver. BhPyCV specific PCR was used to screen available frozen tissues from 13 of these pythons, four additional deceased pythons with and without spinal disease, and fecal samples from 37 live snakes of multiple species with unknown disease status. PCR detected multiple positive tissues in both of the ISH positive BhP and in the feces of another two live BhP and two live annulated tree boas (Corallus annulatus). Preliminary analysis indicates this circovirus can infect BhPs where it was found in 4/5 BhPs tested (2/2 with spinal disease, 2/3 live with unknown status), Boelen’s python (1/2 with spinal disease), and annulated tree boa (2/6 live with unknown status) but was not detected in other python species with the same spinal lesions. This circovirus’ causal or contributory role in spinal disease remains speculative and not well supported by these initial data.

Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics

The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

The Tasmanian devil is an endangered carnivorous marsupial threatened by devil facial tumor disease (DFTD). While research on DFTD has been extensive, little is known about viruses in devils and whether any are of potential conservation relevance for this endangered species. Using both metagenomics based on virion enrichment and sequence-independent amplification (virion-enriched metagenomics) and metatranscriptomics based on bulk RNA sequencing, we characterized and compared the fecal viromes of captive and wild devils. A total of 54 fecal samples collected from two captive and four wild populations were processed for virome characterization using both approaches. In total, 24 novel marsupial-related viruses, comprising a sapelovirus, astroviruses, rotaviruses, picobirnaviruses, parvoviruses, papillomaviruses, polyomaviruses, and a gammaherpesvirus, were identified, as well as known mammalian pathogens such as rabbit hemorrhagic disease virus 2. Captive devils showed significantly lower viral diversity than wild devils. Comparison of the two virus discovery approaches revealed substantial differences in the number and types of viruses detected, with metatranscriptomics better suited for RNA viruses and virion-enriched metagenomics largely identifying more DNA viruses. Thus, the viral communities revealed by virion-enriched metagenomics and metatranscriptomics were not interchangeable and neither approach was able to detect all viruses present. An integrated approach using both virion-enriched metagenomics and metatranscriptomics constitutes a powerful tool for obtaining a complete overview of both the taxonomic and functional profiles of viral communities within a sample.

IMPORTANCE The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

Human sapovirus GI.2 and GI.3 from children with acute gastroenteritis in northern Brazil

Human sapoviruses (HSaV) are considered important causative agents of acute gastroenteritis in humans worldwide. However, knowledge of the genetic characteristics of the whole genome of HSaV in Brazil is limited. Here we report the complete genome sequences of six HSaVs GI.2 and two GI.3 strains obtained from children with acute gastroenteritis in the Northern region of Brazil. Next generation sequencing was used to obtain the full genome and molecular characterization of the genome was performed. Phylogenetic analysis of the genome was also performed. Only one complete HSaV GI.2 genome characterization in the country precedes that of the present study. This is the first complete genome sequence of genotype GI.3 in Brazil. The data obtained in this investigation can contribute to the augmentation of the database on the molecular diversity of HSaVs strains circulating in Brazil, and to the improvement of current typing protocols.

A Highly Divergent Picornavirus Infecting the Gut Epithelia of Zebrafish (Danio rerio) in Research Institutions Worldwide

Zebrafish have been extensively used as a model system for research in vertebrate development and pathogen-host interactions. We describe the complete genome of a novel picornavirus identified during a viral metagenomics analysis of zebrafish gut tissue. The closest relatives of this virus showed identity of <20% in their P1 capsids and <36% in their RdRp qualifying zebrafish picornavirus-1 (ZfPV-1) as member of a novel genus with a proposed name of Cyprivirus. Reverse transcription (RT)-PCR testing of zebrafish from North America, Europe, and Asia showed ZfPV-1 to be globally distributed, being detected in 23 of 41 (56%) institutions tested. In situ hybridization of whole zebrafish showed viral RNA was restricted to a subset of enterocytes and cells in the subjacent lamina propria of the intestine and the intestinal mucosa. This naturally occurring and apparently asymptomatic infection (in wild-type zebrafish lineage AB) provides a natural infection system to study picornavirus-host interactions in an advanced vertebrate model organism. Whether ZfPV-1 infection affects any immunological, developmental, or other biological processes in wild-type or mutant zebrafish lineages remains to be determined.

Illumina sequencing of clinical samples for virus detection in a public health laboratory

High-throughput sequencing (HTS) provides the opportunity, once a diagnostic result is obtained, to extract additional information from a virus-containing sample. Hence, it offers advantages over established quantitative amplification technology, such as quantitative PCR, particularly in a public health environment. At this early stage of its clinical application, there have been limited studies comparing HTS performance to that of the more established quantitative PCR technology for direct detection of viruses. In this pilot-scale study, we tested HTS with a range of viruses and sample types routinely encountered in a public health virology laboratory. In comparison with quantitative PCR, our HTS method was able to sensitively (92%) detect all viruses in any sample type with the exception of certain tissues. Moreover, sufficient nucleotide sequence information was obtained to enable genotyping of strains detected, thus providing additional useful epidemiological information. While HTS sensitivity may not yet match that of PCR, the added value through enhanced epidemiological data has considerable potential to enable real-time surveillance of circulating strains so as to facilitate rapid and appropriate response to outbreaks and virus zoonotic spillover events.

Discovery of Cucumis melo endornavirus by deep sequencing of human stool samples in Brazil

The nearly complete genome sequences of two Cucumis melo endornavirus (CmEV) strains were obtained using deep sequencing while investigating fecal samples for the presence of gastroenteritis viruses. The Brazilian CmEV BRA/TO-23 (aa positions 116-5027) and BRA/TO-74 (aa positions 26-5057) strains were nearly identical to the reference CmEV CL-01 (USA) and SJ1 (South Korea) strains, showing 97% and 98% of nucleotide and amino acid identity, respectively. Endornaviruses are not known to be associated with human disease and their presence may simply reflect recent dietary consumption. Metagenomic analyses offered an opportunity to identify for the first time in Brazil a newly described endornavirus species.

Viral metagenomics reveals significant viruses in the genital tract of apparently healthy dairy cows

The virome in genital tract secretion samples collected from 80 dairy cattle in Shanghai, China, was characterized. Viruses detected included members of the families Papillomaviridae, Polyomaviridae, Hepeviridae, Parvoviridae, Astroviridae, Picornaviridae, and Picobirnaviridae. A member of a new species within the genus Dyoxipapillomavirus and six circular Rep-encoding single-stranded DNA (ssDNA) (CRESS-DNA) viral genomes were fully sequenced and phylogenetically analyzed. The prevalence of bovine polyomaviruses 1 and 2 was measured by PCR to be 10% (8/80) and 6.25% (5/80), respectively. PCR screening also indicated that the novel papillomavirus ujs-21015 and bovine herpesvirus 6 were present in three and two out of the 80 samples, respectively.