Looking into sewage: how far can metagenomics help to detect human enteric viruses?

The impact of human sewage on environmental and food contamination constitutes an important safety issue. Indeed, human sewage reflects the microbiome of the local population, and a variety of human viruses can be detected in wastewater samples. Being able to describe the diversity of viruses present in sewage will provide information on the health of the surrounding population health and will help to prevent further transmission. Metagenomic developments, allowing the description of all the different genomes present in a sample, are very promising tools for virome analysis. However, looking for human enteric viruses with short RNA genomes which are present at low concentrations is challenging. In this study we demonstrate the benefits of performing technical replicates to improve viral identification by increasing contig length, and the set-up of quality criteria to increase confidence in results. Our approach was able to effectively identify some virus sequences and successfully describe the viral diversity. The method yielded full genomes either for norovirus, enterovirus and rotavirus, even if, for these segmented genomes, combining genes remain a difficult issue. Developing reliable viromic methods is important as wastewater sample analysis provides an important tool to prevent further virus transmission by raising alerts in case of viral outbreaks or emergence.

Human Sapovirus Replication in Human Intestinal Enteroids.

Human sapoviruses (HuSaVs), like human noroviruses (HuNoV), belong to the Caliciviridae family and cause acute gastroenteritis in humans. Since their discovery in 1976, numerous attempts to grow HuSaVs in vitro were unsuccessful until 2020, when these viruses were reported to replicate in a duodenal cancer cell-derived line. Physiological cellular models allowing viral replication are essential to investigate HuSaV biology and replication mechanisms such as genetic susceptibility, restriction factors, and immune responses to infection. In this study, we demonstrate replication of two HuSaV strains in human intestinal enteroids (HIEs) known to support the replication of HuNoV and other human enteric viruses. HuSaVs replicated in differentiated HIEs originating from jejunum, duodenum and ileum, but not from the colon, and bile acids were required. Between 2h and 3 to 6 days postinfection, viral RNA levels increased up from 0.5 to 1.8 log10-fold. Importantly, HuSaVs were able to replicate in HIEs independent of their secretor status and histo-blood group antigen expression. The HIE model supports HuSaV replication and allows a better understanding of host-pathogen mechanisms such as cellular tropism and mechanisms of viral replication. IMPORTANCE Human sapoviruses (HuSaVs) are a frequent but overlooked cause of acute gastroenteritis, especially in children. Little is known about this pathogen, whose successful in vitro cultivation was reported only recently, in a cancer cell-derived line. Here, we assessed the replication of HuSaV in human intestinal enteroids (HIEs), which are nontransformed cultures originally derived from human intestinal stem cells that can be grown in vitro and are known to allow the replication of other enteric viruses. Successful infection of HIEs with two strains belonging to different genotypes of the virus allowed discovery that the tropism of these HuSaVs is restricted to the small intestine, does not occur in the colon, and replication requires bile acid but is independent of the expression of histo-blood group antigens. Thus, HIEs represent a physiologically relevant model to further investigate HuSaV biology and a suitable platform for the future development of vaccines and antivirals.

SARS-CoV-2 Whole-Genome Sequencing Using Oxford Nanopore Technology for Variant Monitoring in Wastewaters.

Since the beginning of the Coronavirus Disease-19 (COVID-19) pandemic, multiple Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) mutations have been reported and led to the emergence of variants of concern (VOC) with increased transmissibility, virulence or immune escape. In parallel, the observation of viral fecal shedding led to the quantification of SARS-CoV-2 genomes in wastewater, providing information about the dynamics of SARS-CoV-2 infections within a population including symptomatic and asymptomatic individuals. Here, we aimed to adapt a sequencing technique initially designed for clinical samples to apply it to the challenging and mixed wastewater matrix, and hence identify the circulation of VOC at the community level. Composite raw sewage sampled over 24 h in two wastewater-treatment plants (WWTPs) from a city in western France were collected weekly and SARS-CoV-2 quantified by RT-PCR. Samples collected between October 2020 and May 2021 were submitted to whole-genome sequencing (WGS) using the primers and protocol published by the ARTIC Network and a MinION Mk1C sequencer (Oxford Nanopore Technologies, Oxford, United Kingdom). The protocol was adapted to allow near-full genome coverage from sewage samples, starting from ∼5% to reach ∼90% at depth 30. This enabled us to detect multiple single-nucleotide variant (SNV) and assess the circulation of the SARS-CoV-2 VOC Alpha, Beta, Gamma, and Delta. Retrospective analysis of sewage samples shed light on the emergence of the Alpha VOC with detection of first co-occurring signature mutations in mid-November 2020 to reach predominance of this variant in early February 2021. In parallel, a mutation-specific qRT-PCR assay confirmed the spread of the Alpha VOC but detected it later than WGS. Altogether, these data show that SARS-CoV-2 sequencing in sewage can be used for early detection of an emerging VOC in a population and confirm its ability to track shifts in variant predominance.

Development of passive samplers for the detection of SARS-CoV-2 in sewage and seawater: Application for the monitoring of sewage.

Recent studies have shown that passive sampling is a promising tool for SARS-CoV-2 detection for wastewater-based epidemiology (WBE) application. We have previously developed passive sampling of viruses using polymer membranes in seawater. Even though SARS-CoV-2 was not detected yet in seawater, passive sampling could be optimized for future application in coastal areas close to wastewater treatment plant (WWTP). The aim of this study was to optimize passive sampling of SARS-CoV-2 in sewage and seawater by selecting a suitable membrane, to determine whether the quantities of virus increase over time, and then to determine if passive sampling and traditional sampling are correlated when conducted in a wastewater treatment plant. Nylon and Zetapor allowed the detection of heat inactivated SARS-CoV-2 and of the Porcine Epidemic Diarrhea Virus (PEDV), a coronavirus surrogate, in wastewater and seawater spiked with these 2 viruses, showing an increase in detection between 4 h and 24 h of immersion and significantly higher recoveries of both viruses with nylon in seawater (15%) compared to wastewater (4%). On wastewater samples, both membranes detected the virus, the recovery rate was of about 3% for freshly collected samples, and no significant difference was found between SARS-CoV-2 genome concentration on Zetapor and that in water. In sewage spiked seawater, similar concentrations of genome were found on both membranes, with a mean recovery rate of 16% and 11% respectively for nylon and Zetapor. A 3-weeks monitoring with passive sampler allowed the detection of viruses in the influent of a WWTP with a frequency of 100% and 76% for SARS-CoV-2 and norovirus GII respectively. Passive and traditional sampling gave the same evolution of the SARS-CoV-2 concentration over time. All these results confirmed the interest of passive sampling for virus detection and its potential application for monitoring in the wastewater system for targeted public health actions.

Human and Animal RNA Virus Diversity Detected by Metagenomics in Cameroonian Clams.

Many recent pandemics have been recognized as zoonotic viral diseases. While their origins remain frequently unknown, environmental contamination may play an important role in emergence. Thus, being able to describe the viral diversity in environmental samples contributes to understand the key issues in zoonotic transmission. This work describes the use of a metagenomic approach to assess the diversity of eukaryotic RNA viruses in river clams and identify sequences from human or potentially zoonotic viruses. Clam samples collected over 2years were first screened for the presence of norovirus to verify human contamination. Selected samples were analyzed using metagenomics, including a capture of sequences from viral families infecting vertebrates (VirCapSeq-VERT) before Illumina NovaSeq sequencing. The bioinformatics analysis included pooling of data from triplicates, quality filtering, elimination of bacterial and host sequences, and a deduplication step before de novo assembly. After taxonomic assignment, the viral fraction represented 0.8-15% of reads with most sequences (68-87%) remaining un-assigned. Yet, several mammalian RNA viruses were identified. Contigs identified as belonging to the Astroviridae were the most abundant, with some nearly complete genomes of bastrovirus identified. Picobirnaviridae sequences were related to strains infecting bats, and few others to strains infecting humans or other hosts. Hepeviridae sequences were mostly related to strains detected in sponge samples but also strains from swine samples. For Caliciviridae and Picornaviridae, most of identified sequences were related to strains infecting bats, with few sequences close to human norovirus, picornavirus, and genogroup V hepatitis A virus. Despite a need to improve the sensitivity of our method, this study describes a large diversity of RNA virus sequences from clam samples. To describe all viral contaminants in this type of food, and being able to identify the host infected by viral sequences detected, may help to understand some zoonotic transmission events and alert health authorities of possible emergence.

Metagenomics-Based Proficiency Test of Smoked Salmon Spiked with a Mock Community.

An inter-laboratory proficiency test was organized to assess the ability of participants to perform shotgun metagenomic sequencing of cold smoked salmon, experimentally spiked with a mock community composed of six bacteria, one parasite, one yeast, one DNA, and two RNA viruses. Each participant applied its in-house wet-lab workflow(s) to obtain the metagenomic dataset(s), which were then collected and analyzed using MG-RAST. A total of 27 datasets were analyzed. Sample pre-processing, DNA extraction protocol, library preparation kit, and sequencing platform, influenced the abundance of specific microorganisms of the mock community. Our results highlight that despite differences in wet-lab protocols, the reads corresponding to the mock community members spiked in the cold smoked salmon, were both detected and quantified in terms of relative abundance, in the metagenomic datasets, proving the suitability of shotgun metagenomic sequencing as a genomic tool to detect microorganisms belonging to different domains in the same food matrix. The implementation of standardized wet-lab protocols would highly facilitate the comparability of shotgun metagenomic sequencing dataset across laboratories and sectors. Moreover, there is a need for clearly defining a sequencing reads threshold, to consider pathogens as detected or undetected in a food sample.

A Targeted Metagenomics Approach to Study the Diversity of Norovirus GII in Shellfish Implicated in Outbreaks.

Human noroviruses (NoV) cause epidemics of acute gastroenteritis (AGE) worldwide and can be transmitted through consumption of contaminated foods. Fresh products such as shellfish can be contaminated by human sewage during production, which results in the presence of multiple virus strains, at very low concentrations. Here, we tested a targeted metagenomics approach by deep-sequencing PCR amplicons of the capsid (VP1) and polymerase (RdRp) viral genes, on a set of artificial samples and on shellfish samples associated to AGE outbreaks, to evaluate its advantages and limitations in the identification of strains from the NoV genogroup (G) II. Using artificial samples, the method allowed the sequencing of most strains, but not all, and displayed variability between replicates especially with lower viral concentrations. Using shellfish samples, targeted metagenomics was compared to Sanger-sequencing of cloned amplicons and was able to identify a higher diversity of NoV GII and GIV strains. It allowed phylogenetic analyses of VP1 sequences and the identification, in most samples, of GII.17[P17] strains, also identified in related clinical samples. Despite several limitations, combining RdRp- and VP1-targeted metagenomics is a sensitive approach allowing the study NoV diversity in low-contaminated foods and the identification of NoV strains implicated in outbreaks.

Contamination of Clams with Human Norovirus and a Novel Hepatitis A Virus in Cameroon.

Shellfish constitute an important protein source but may be contaminated by viruses from various origins. A study performed on clams collected in Cameroon showed a high prevalence of norovirus and hepatitis A virus. After sequencing, the hepatitis A virus showed similarities with the genotype V simian strains.

Metagenomic to evaluate norovirus genomic diversity in oysters: Impact on hexamer selection and targeted capture-based enrichment.

Human virus transmission through food consumption has been identified since many years and the international trade increases the risk of dissemination of viral pathogens. The development of metagenomic approach holds many promises for the surveillance of viruses in food and water. This work aimed to analyze norovirus diversity and to evaluate strain-dependent accumulation patterns in three oyster types by using a metagenomic approach. Different hexamer sets to prime cDNA were evaluated before capture-based approach to enhance virus reads recovery during deep sequencing. The study includes the use of technical replicates of artificially contaminated oysters and the analysis of multiple negatives controls. Results showed a clear impact of the hexamer set used for cDNA synthesis. A set of In-house designed (I-HD) hexamers, selected to lower mollusk amplification, gave promising results in terms of viral reads abundancy. However, the best correlation between CT values, thus concentrations, and number of reads was observed using random hexamers. Random hexamers also provided the highest numbers of reads and allowed the identification of sequence of different human enteric viruses. Regarding human norovirus, different genogroups and genotypes were identified among contigs longer than 500 bp. Two full genomes and six sequences longer than 3600 bases were obtained allowing a precise strain identification. The use of technical triplicates was found valuable to increase the chances to sequence viral strains present at low concentrations. Analyzing viral contamination in shellfish samples is quite challenging, however this work demonstrates that the recovery of full genome or long contigs, allowing clear identification of viral strains is possible.

Complete Genome Sequence of Campylobacter armoricus CA639, Which Carries Two Plasmids, Compiled Using Oxford Nanopore and Illumina Sequencing Technologies.

As determined by a hybrid approach combining Oxford Nanopore MinION and Illumina MiniSeq sequence data, Campylobacter armoricus strain CA639 harbored a circular chromosome of 1,688,169 bp with a G+C content of 28.47% and two plasmids named pCA639-1 and pCA639-2, with lengths of 51,123 and 28,139 bp, and G+C contents of 26.5% and 28.45%, respectively.

Metavirome Sequencing to Evaluate Norovirus Diversity in Sewage and Related Bioaccumulated Oysters.

Metagenomic sequencing is a promising method to determine the virus diversity in environmental samples such as sewage or shellfish. However, to identify the short RNA genomes of human enteric viruses among the large diversity of nucleic acids present in such complex matrices, method optimization is still needed. This work presents methodological developments focused on norovirus, a small ssRNA non-enveloped virus known as the major cause of human gastroenteritis worldwide and frequently present in human excreta and sewage. Different elution protocols were applied and Illumina MiSeq technology were used to study norovirus diversity. A double approach, agnostic deep sequencing and a capture-based approach (VirCapSeq-VERT) was used to identify norovirus in environmental samples. Family-specific viral contigs were classified and sorted by SLIM and final norovirus contigs were genotyped using the online Norovirus genotyping tool v2.0. From sewage samples, 14 norovirus genogroup I sequences were identified of which six were complete genomes. For norovirus genogroup II, nine sequences were identified and three of them comprised more than half of the genome. In oyster samples bioaccumulated with these sewage samples, only the use of an enrichment step during library preparation allowed successful identification of nine different sequences of norovirus genogroup I and four for genogroup II (>500 bp). This study demonstrates the importance of method development to increase virus recovery, and the interest of a capture-based approach to be able to identify viruses present at low concentrations.

Characterization of Norovirus and Other Human Enteric Viruses in Sewage and Stool Samples Through Next-Generation Sequencing.

This study aimed to optimize a method to identify human enteric viruses in sewage and stool samples using random primed next-generation sequencing. We tested three methods, two employed virus enrichment based on the binding properties of the viral capsid using pig-mucin capture or by selecting viral RNA prior to library preparation through a capture using the SureSelect target enrichment. The third method was based on a non-specific biophysical precipitation with polyethylene glycol. Full genomes of a number of common human enteric viruses including norovirus, rotavirus, husavirus, enterovirus and astrovirus were obtained. In stool samples full norovirus genome were detected as well as partial enterovirus genome. A variety of norovirus sequences was detected in sewage samples, with genogroup II being more prevalent. Interestingly, the pig-mucin capture enhanced not only the recovery of norovirus and rotavirus but also recovery of astrovirus, sapovirus and husavirus. Documenting sewage virome using these methods provides information for molecular epidemiology and may be useful in developing strategies to prevent further spread of viruses.

Improving the efficacy of sewage treatment decreases norovirus contamination in oysters.

As human population increases worldwide, water quality will become increasingly problematic, and food consumed raw may be of higher risk. This is already evident for oysters grown in coastal areas - despite regulations based on bacterial indicators, oysters are still implicated in food-borne outbreaks worldwide. The pathogens most frequently detected are human noroviruses, which are shed at high concentrations in human excreta and are very resistant to environmental conditions. Sewage treatment plants usually apply a variety of steps such as activated sludge treatment, chlorine or UV disinfection to eliminate contaminants, these processes have variable efficacy. This study demonstrates the impact of replacing an old lagoon-based sewage treatment plant with a new membrane bioreactor sewage treatment plant on human norovirus levels in treated sewage and oysters. While comparable norovirus concentrations were detected in the influent samples, a clear difference was observed in effluent quality, as norovirus was only detected in one sample after treatment in the new membrane bioreactor system, confirming the efficiency of such technology. As a direct impact, oysters located close to the membrane bioreactor sewage outfall were less frequently contaminated by norovirus, and showed lower concentrations compared to the first period of the study when they were exposed to sewage effluent from the lagoon outfall. Shellfish located upstream showed comparable contamination levels suggesting that there are also other sources of norovirus contamination in the estuary. Considering the health benefits of shellfish consumption, improving wastewater quality will make an important contribution to enhancing the safety of shellfish and international food security.

Infectivity and RNA Persistence of a Norovirus Surrogate, the Tulane Virus, in Oysters.

Oysters, being filter feeders, can accumulate some human pathogens such as norovirus, a highly infectious calicivirus, most common cause of acute gastroenteritis worldwide. Accumulated virus decays over a period of days to weeks, possibly rendering contaminated oysters safe again. Sensitive molecular methods have been set up for shellfish analysis but without answering the question of infectious virus detection. Using the Tulane virus (TV), a norovirus surrogate that recognizes the same ligand as human norovirus in oyster tissues, the genome and infectious virus decay rates were estimated using inverse linear regression in a Bayesian framework for genome copies. Infectivity decreased faster than genome copies but infectious viruses were detected for several days. Quantifying the decrease in viral infectivity and genome detection in oysters over such a long period may help local authorities to manage production areas implicated in shellfish-borne outbreaks, and thus protect consumers.

Virus Type-Specific Removal in a Full-Scale Membrane Bioreactor Treatment Process.

We investigated removal of noroviruses, sapoviruses, and rotaviruses in a full-scale membrane bioreactor (MBR) plant by monitoring virus concentrations in wastewater samples during two gastroenteritis seasons and evaluating the adsorption of viruses to mixed liquor suspended solids (MLSS). Sapoviruses and rotaviruses were detected in 25% of MBR effluent samples with log reduction values of 3- and 2-logs in geometric mean concentrations, respectively, while noroviruses were detected in only 6% of the samples. We found that norovirus and sapovirus concentrations in the solid phase of mixed liquor samples were significantly higher than in the liquid phase (P < 0.01, t test), while the concentration of rotaviruses was similar in both phases. The efficiency of adsorption of the rotavirus G1P[8] strain to MLSS was significantly less than norovirus GI.1 and GII.4 and sapovirus GI.2 strains (P < 0.01, t test). Differences in the adsorption of viruses to MLSS may cause virus type-specific removal during the MBR treatment process as shown by this study.

Digital PCR for Quantifying Norovirus in Oysters Implicated in Outbreaks, France.

Using samples from oysters clearly implicated in human disease, we quantified norovirus levels by using digital PCR. Concentrations varied from 43 to 1,170 RNA copies/oyster. The analysis of frozen samples from the production area showed the presence of norovirus 2 weeks before consumption.

Tulane Virus as a Potential Surrogate To Mimic Norovirus Behavior in Oysters.

Oyster contamination by noroviruses is an important health and economic problem. The present study aimed to compare the behaviors of Norwalk virus (the prototype genogroup I norovirus) and two culturable viruses: Tulane virus and mengovirus. After bioaccumulation, tissue distributions were quite similar for Norwalk virus and Tulane virus, with the majority of viral particles detected in digestive tissues, while mengovirus was detected in large amounts in the gills and mantle as well as in digestive tissues. The levels of persistence of all three viruses over 8 days were comparable, but clear differences were observed over longer periods, with Norwalk and Tulane viruses displaying rather similar half-lives, unlike mengovirus, which was cleared more rapidly. These results indicate that Tulane virus may be a good surrogate for studying norovirus behavior in oysters, and they confirm the prolonged persistence of Norwalk virus in oyster tissues.

Bioaccumulation efficiency, tissue distribution, and environmental occurrence of hepatitis E virus in bivalve shellfish from France.

Hepatitis E virus (HEV), an enteric pathogen of both humans and animals, is excreted by infected individuals and is therefore present in wastewaters and coastal waters. As bivalve molluscan shellfish are known to concentrate viral particles during the process of filter feeding, they may accumulate this virus. The bioaccumulation efficiencies of oysters (Crassostrea gigas), flat oysters (Ostrea edulis), mussels (Mytilus edulis), and clams (Ruditapes philippinarum) were compared at different time points during the year. Tissue distribution analysis showed that most of the viruses were concentrated in the digestive tissues of the four species. Mussels and clams were found to be more sensitive to sporadic contamination events, as demonstrated by rapid bioaccumulation in less than 1 h compared to species of oysters. For oysters, concentrations increased during the 24-h bioaccumulation period. Additionally, to evaluate environmental occurrence of HEV in shellfish, an environmental investigation was undertaken at sites potentially impacted by pigs, wild boars, and human waste. Of the 286 samples collected, none were contaminated with hepatitis E virus, despite evidence that this virus is circulating in some French areas. It is possible that the number of hepatitis E viral particles discharged into the environment is too low to detect or that the virus may have a very short period of persistence in pig manure and human waste.