Inter- and Intra-Host Nucleotide Variations in Hepatitis A Virus in Culture and Clinical Samples Detected by Next-Generation Sequencing

Whole-Genome Sequencing-Based Confirmatory Methods on RT-qPCR Results for the Detection of Foodborne Viruses in Frozen Berries

Accurate detection, identification, and subsequent confirmation of pathogens causing foodborne illness are essential for the prevention and investigation of foodborne outbreaks. This is particularly true when the causative agent is an enteric virus that has a very low infectious dose and is likely to be present at or near the limit of detection. In this study, whole-genome sequencing (WGS) was combined with either of two non-targeted pre-amplification methods (SPIA and SISPA) to investigate their utility as a confirmatory method for RT-qPCR positive results of foods contaminated with enteric viruses. Frozen berries (raspberries, strawberries, and blackberries) were chosen as the food matrix of interest due to their association with numerous outbreaks of foodborne illness. The hepatitis A virus (HAV) and human norovirus (HuNoV) were used as the contaminating agents. The non-targeted WGS strategy employed in this study could detect and confirm HuNoV and HAV at genomic copy numbers in the single digit range, and in a few cases, identified viruses present in samples that had been found negative by RT-qPCR analyses. However, some RT-qPCR-positive samples could not be confirmed using the WGS method, and in cases with very high Ct values, only a few viral reads and short sequences were recovered from the samples. WGS techniques show great potential for confirmation and identification of virally contaminated food items. The approaches described here should be further optimized for routine application to confirm the viral contamination in berries.

A combined approach of rolling-circle amplification-single site restriction endonuclease digestion followed by next generation sequencing to characterize the whole genome and intra-host variants of human Torque teno virus

Torque Teno Virus (TTV) was initially associated with post-transfusion hepatitis, but growing evidence of its ubiquity in humans is compatible to no apparent clinical significance. TTV is a small non-enveloped virus with a circular single-negative-stranded DNA genome, belonging to the Anelloviridae family. Currently, TTVs are divided in seven phylogenetic groups and are further classified into 21 species. Studies about diversity of TTV in different conditions are receiving increasing interest and in this sense, sequencing of whole genomes for better genetic characterization becomes even more important. Since its discovery in 1997, few TTV complete genomes have been reported worldwide. This is probably due, among other reasons, to the great genetic heterogeneity among TTV strains that prevents its amplification and sequencing by conventional PCR and cloning methods. In addition, although metagenomics approach is useful in these cases, it remains a challenging tool for viromic analysis. With the aim of contributing to the expansion of the TTV whole genomes dataset and to study intra-host variants, we employed a methodology that combined a rolling-circle amplification approach followed by EcoRI digestion, generating a DNA fragment of ∼4Kb consistent with TTV genome length which was sequenced by Illumina next generation sequencing. A genogroup 3 full-length consensus TTV genome was obtained and co-infection with other species (at least those with a single EcoRI cleavage site) was not identified. Additionally, bioinformatics analysis allowed to identify the spectrum of TTV intra-host variants which provides evidence of a complex evolution dynamics of these DNA circular viruses, similarly to what occurs with RNA viruses.

Recovery of Recombinant Canine Distemper Virus That Expresses CPV-2a VP2: Uncovering the Mutation Profile of Recombinant Undergoing 50 Serial Passages In Vitro

Canine distemper and canine parvoviral enteritis are infections caused by the canine distemper virus (CDV) and canine parvovirus type 2 (CPV-2), respectively. They are two common infectious diseases that cause high morbidity and mortality in affected dogs. Combination vaccines have been broadly used to protect dogs from infections of CDV, CPV-2, and other viruses. VP2 is the most abundant protein of the CPV-2 capsid. It elicits potent immunity in animals and, therefore, is widely used for designing subunit antigen-based vaccines. In this study, we rescued a recombinant CDV (QN vaccine strain) using reverse genetics. The recombinant CDV (rCDV-VP2) was demonstrated to express stably the VP2 in cells for at least 33 serial passages in vitro. Unfortunately, a nonsense mutation was initially identified in the VP2 open reading frame (ORF) at passage-34 (P34) and gradually became predominant in rCDV-VP2 quasispecies with passaging. Neither test strip detection nor indirect immunofluorescence assay demonstrated the expression of the VP2 at P50. The P50 rCDV-VP2 was subjected to next-generation sequencing, which totally identified 17 single-nucleotide variations (SNVs), consisting of 11 transitions and 6 transversions. Out of the 17 SNVs, 1 and 9 were identified as nonsense and missense mutations, respectively. Since the nonsense mutation arose in the VP2 ORF as early as P34, an earlier rCDV-VP2 progeny should be selected for the vaccination of animals in future experiments.

Rescuing eGFP-Tagged Canine Distemper Virus for 40 Serial Passages Separately in Ribavirin- and Non-Treated Cells: Comparative Analysis of Viral Mutation Profiles

Due to lacking a proofreading mechanism in their RNA-dependent RNA polymerases (RdRp), RNA viruses generally possess high mutation frequencies, making them evolve rapidly to form viral quasispecies during serial passages in cells, especially treated with mutagens, like ribavirin. Canine distemper virus (CDV) belongs to the genus Morbillivirus. Its L protein functions as an RdRp during viral replication. In this study, a recombinant enhanced green fluorescence protein-tagged CDV (rCDV-eGFP) was rescued from its cDNA clone, followed by viral identification and characterization at passage-7 (P7). This recombinant was independently subjected to extra 40 serial passages (P8 to 47) in ribavirin- and non-treated cells. Two viral progenies, undergoing passages in ribavirin- and non-treated VDS cells, were named rCDV-eGFP-R and -N, respectively. Both progenies were simultaneously subjected to next-generation sequencing (NGS) at P47 for comparing their quasispecies diversities with each other. The rCDV-eGFP-R and -N showed 62 and 23 single-nucleotide mutations (SNMs) in individual antigenomes, respectively, suggesting that the ribavirin conferred a mutagenic effect on the rCDV-eGFP-R. The spectrum of 62 SNMs contained 26 missense and 36 silent mutations, and that of 23 SNMs was composed of 17 missense and 6 silent mutations. Neither the rCDV-eGFP-R nor -N exhibited nonsense mutation in individual antigenomes. We speculate that the rCDV-eGFP-R may contain at least one P47 sub-progeny characterized by high-fidelity replication in cells. If such a sub-progeny can be purified from the mutant swarm, its L protein would elucidate a molecular mechanism of CDV high-fidelity replication.

Mutation Profiles of eGFP-Tagged Small Ruminant Morbillivirus During 45 Serial Passages in Ribavirin-Treated Cells

Small ruminant morbillivirus (SRMV), formerly known as peste-des-petits-ruminants virus, classified into the genus Morbillivirus in the family Paramyxoviridae. Its L protein functions as the RNA-dependent RNA polymerases (RdRp) during viral replication. Due to the absence of efficient proofreading activity in their RdRps, various RNA viruses reveal high mutation frequencies, making them evolve rapidly during serial passages in cells, especially treated with a certain mutagen, like ribavirin. We have previously rescued a recombinant enhanced green fluorescence protein-tagged SRMV (rSRMV-eGFP) using reverse genetics. In this study, the rSRMV-eGFP was subjected to serial passages in ribavirin-treated cells. Due to the ribavirin-exerted selective pressure, it was speculated that viral progenies would form quasispecies after dozens of passages. Viral progenies at passage-10, -20, -30, -40, and -50 were separately subjected to next-generation sequencing (NGS), consequently revealing a total of 34 single-nucleotide variations, including five synonymous, 21 missense, and one non-sense mutations. The L sequence was found to harbor eight missense mutations during serial passaging. It was speculated that at least one high-fidelity variant was present in viral quasispecies at passage-50. If purified from the population of viral progenies, this putative variant would contribute to clarifying a molecular mechanism in viral high-fidelity replication in vitro.

Near-Complete Genome Sequence of a Human Norovirus GII.1[Pg] Strain Associated with Acute Gastroenteritis, Determined Using Long-Read Sequencing

High-throughput sequencing is one of the approaches used for the detection of foodborne pathogens such as noroviruses. Long-read sequencing has advantages over short-read sequencing in speed, read length, and lower fragmentation bias, which makes it a potential powerful tool for the fast detection and identification of viruses. Using Nanopore sequencing technology, we were able to successfully recover a nearly complete genome sequence of a human norovirus GII.1[Pg] strain in a single long read from a sample from a patient with norovirus gastroenteritis.

Application of whole-genome sequencing for norovirus outbreak tracking and surveillance efforts in Orange County, CA

Noroviruses are the leading cause of acute gastroenteritis and foodborne illness in the United States. Traditional Sanger sequencing of short genomic regions (~300-600 bp) is the primary method for differentiation of this pathogen; however, whole-genome sequencing (WGS) offers a valuable approach to further characterize strains of this virus. The objective of this study was to investigate the ability of WGS compared to Sanger sequencing to differentiate norovirus strains and enhance outbreak investigation and surveillance efforts. WGS results for 41 norovirus-positive stool samples from 15 different outbreaks occurring from 2012 to 2019 in Orange County, CA, were analyzed for this study. All samples were genotyped with both WGS and Sanger sequencing based on the B-C region. WGS generated nearly full-length viral genome sequences (7029-7768 bp) with 4x to 35,378x coverage. Phylogenetic analysis of WGS data enabled differentiation of genotypically similar strains from separate outbreaks. Single nucleotide variation (SNV) analysis on a subset of strains revealed nucleotide variations (15-79 nt) among isolates from multiple outbreaks of GII.4 Sydney_2015[P31] and GII.17[P17]. Overall, the results demonstrated that coupling norovirus genotype identification with WGS enables enhanced genetic differentiation of strains and provides valuable information for outbreak investigation and surveillance efforts.

Rescue of Senecavirus A to uncover mutation profiles of its progenies during 80 serial passages in vitro

Senecavirus A (SVA), also known as Seneca Valley virus, belongs to the genus Senecavirus in the family Picornaviridae. In this study, a China SVA isolate (CH-LX-01-2016) was rescued from its cDNA clone, and then identified by RT-PCR, indirect immunofluorescence assay and mass spectrometry. The rescued SVA could separately induce typical plaque formations and cytopathic effects in cell monolayers. In order to uncover its evolutionary dynamics, the SVA was subjected to eighty serial passages in vitro. Its progenies per ten passages were analyzed by next-generation sequencing (NGS). The NGS analyses showed that neither sequence-deleting nor -inserting phenotype was detectable in eight progenies, within which a total of forty-one intra-host single-nucleotide variations (SNVs) arose with passaging. Almost all SNVs were identified as the single-nucleotide polymorphism with mixture of two nucleotides. SNVs led to eighteen nonsynonymous mutations, out of which sixteen could directly reflect their own frequencies of amino acid mutation, due to only one SNV occurring in their individual codons. Compared with its parental virus without passaging, the passage-80 SVA progeny had formed a viral quasispecies, as evidenced by a total of twenty-eight SNVs identified in it.