Time-series sewage metagenomics distinguishes seasonal, human-derived and environmental microbial communities potentially allowing source-attributed surveillance.

Sewage metagenomics has risen to prominence in urban population surveillance of pathogens and antimicrobial resistance (AMR). Unknown species with similarity to known genomes cause database bias in reference-based metagenomics. To improve surveillance, we seek to recover sewage genomes and develop a quantification and correlation workflow for these genomes and AMR over time. We use longitudinal sewage sampling in seven treatment plants from five major European cities to explore the utility of catch-all sequencing of these population-level samples. Using metagenomic assembly methods, we recover 2332 metagenome-assembled genomes (MAGs) from prokaryotic species, 1334 of which were previously undescribed. These genomes account for ~69% of sequenced DNA and provide insight into sewage microbial dynamics. Rotterdam (Netherlands) and Copenhagen (Denmark) show strong seasonal microbial community shifts, while Bologna, Rome, (Italy) and Budapest (Hungary) have occasional blooms of Pseudomonas-dominated communities, accounting for up to ~95% of sample DNA. Seasonal shifts and blooms present challenges for effective sewage surveillance. We find that bacteria of known shared origin, like human gut microbiota, form communities, suggesting the potential for source-attributing novel species and their ARGs through network community analysis. This could significantly improve AMR tracking in urban environments.

Rise and fall of SARS-CoV-2 variants in Rotterdam: Comparison of wastewater and clinical surveillance.

Monitoring of SARS-CoV-2 in wastewater (WW) is a promising tool for epidemiological surveillance, correlating not only viral RNA levels with the infection dynamics within the population, but also to viral diversity. However, the complex mixture of viral lineages in WW samples makes tracking of specific variants or lineages circulating in the population a challenging task. We sequenced sewage samples of 9 WW-catchment areas within the city of Rotterdam, used specific signature mutations from individual SARS-CoV-2 lineages to estimate their relative abundances in WW and compared them against those observed in clinical genomic surveillance of infected individuals between September 2020 and December 2021. We showed that especially for dominant lineages, the median of the frequencies of signature mutations coincides with the occurrence of those lineages in Rotterdam's clinical genomic surveillance. This, along with digital droplet RT-PCR targeting signature mutations of specific variants of concern (VOCs), showed that several VOCs emerged, became dominant and were replaced by the next VOC in Rotterdam at different time points during the study. In addition, single nucleotide variant (SNV) analysis provided evidence that spatio-temporal clusters can also be discerned from WW samples. We were able to detect specific SNVs in sewage, including one resulting in the Q183H amino acid change in the Spike gene, that was not captured by clinical genomic surveillance. Our results highlight the potential use of WW samples for genomic surveillance, increasing the set of epidemiological tools to monitor SARS-CoV-2 diversity.

Profiling of humoral immune responses to norovirus in children across Europe.

Norovirus is a leading cause of epidemic acute gastroenteritis. More than 30 genotypes circulate in humans, some are common, and others are only sporadically detected. Here, we investigated whether serology can be used to determine which genotypes infect children. We established a multiplex protein microarray with structural and non-structural norovirus antigens that allowed simultaneous antibody testing against 30 human GI and GII genotypes. Antibody responses of sera obtained from 287 children aged < 1 month to 5.5 years were profiled. Most specific IgG and IgA responses were directed against the GII.2, GII.3, GII.4, and GII.6 capsid genotypes. While we detected antibody responses against rare genotypes, we found no evidence for wide circulation. We also detected genotype-specific antibodies against the non-structural proteins p48 and p22 in sera of older children. In this study, we show the age-dependent antibody responses to a broad range of norovirus capsid and polymerase genotypes, which will aid in the development of vaccines.

Clinical and In Vitro Evidence Favoring Immunoglobulin Treatment of a Chronic Norovirus Infection in a Patient With Common Variable Immunodeficiency.

BACKGROUND: Immunocompromised individuals can become chronically infected with norovirus, but effective antiviral therapies are not yet available. METHODS: Treatments with nitazoxanide, ribavirin, interferon alpha-2a, and nasoduodenally administered immunoglobulins were evaluated sequentially in an immunocompromised patient chronically infected with norovirus. In support, these components were also applied to measure norovirus inhibition in intestinal enteroid cultures in vitro. Viral RNA levels were determined in fecal and plasma samples during each treatment and viral genomes were sequenced. RESULTS: None of the antivirals resulted in a reduction of viral RNA levels in feces or plasma. However, during ribavirin treatment, there was an increased accumulation of virus genome mutations. In vitro, an effect of interferon alpha-2a on virus replication was observed and a genetically related strain was neutralized effectively in vitro using immunoglobulins and post-norovirus-infection antiserum. In agreement, after administration of immunoglobulins, the patient cleared the infection. CONCLUSIONS: Intestinal enteroid cultures provide a relevant system to evaluate antivirals and the neutralizing potential of immunoglobulins. We successfully treated a chronically infected patient with immunoglobulins, despite varying results reported by others. This case study provides in-depth, multifaceted exploration of norovirus treatment that can be used as a guidance for further research towards norovirus treatments.

Microcalorimetry: A Novel Application to Measure In Vitro Phage Susceptibility of Staphylococcus aureus in Human Serum.

Infections involving antibiotic resistant Staphylococcus aureus (S. aureus) represent a major challenge to successful treatment. Further, although bacteriophages (phages) could be an alternative to antibiotics, there exists a lack of correlation in phage susceptibility results between conventional in vitro and in vivo assays. This discrepancy may hinder the potential implementation of bacteriophage therapy. In this study, the susceptibility of twelve S. aureus strains to three commercial phage cocktails and two single phages was assessed. These S. aureus strains (including ten clinical isolates, five of which were methicillin-resistant) were compared using four assays: the spot test, efficiency of plating (EOP), the optical density assay (all in culture media) and microcalorimetry in human serum. In the spot test, EOP and optical density assay, all cocktails and single phages lysed both methicillin susceptible and methicillin resistant S. aureus strains. However, there was an absence of phage-mediated lysis in high concentrations of human serum as measured using microcalorimetry. As this microcalorimetry-based assay more closely resembles in vivo conditions, we propose that microcalorimetry could be included as a useful addition to conventional assays, thereby facilitating more accurate predictions of the in vivo susceptibility of S. aureus to phages during phage selection for therapeutic purposes.

Mutagenesis of the Varicella-Zoster Virus Genome Demonstrates That VLT and VLT-ORF63 Proteins Are Dispensable for Lytic Infection.

Primary varicella-zoster virus (VZV) infection leads to varicella and the establishment of lifelong latency in sensory ganglion neurons. Reactivation of latent VZV causes herpes zoster, which is frequently associated with chronic pain. Latent viral gene expression is restricted to the VZV latency-associated transcript (VLT) and VLT-ORF63 (VLT63) fusion transcripts. Since VLT and VLT63 encode proteins that are expressed during lytic infection, we investigated whether pVLT and pVLT-ORF63 are essential for VZV replication by performing VZV genome mutagenesis using CRISPR/Cas9 and BAC technologies. We first established that CRISPR/Cas9 can efficiently mutate VZV genomes in lytically VZV-infected cells through targeting non-essential genes ORF8 and ORF11 and subsequently show recovery of viable mutant viruses. By contrast, the VLT region was markedly resistant to CRISPR/Cas9 editing. Whereas most mutants expressed wild-type or N-terminally altered versions of pVLT and pVLT-ORF63, only a minority of the resulting mutant viruses lacked pVLT and pVLT-ORF63 coding potential. Growth curve analysis showed that pVLT/pVLT-ORF63 negative viruses were viable, but impaired in growth in epithelial cells. We confirmed this phenotype independently using BAC-derived pVLT/pVLT-ORF63 negative and repaired viruses. Collectively, these data demonstrate that pVLT and/or pVLT-ORF63 are dispensable for lytic VZV replication but promote efficient VZV infection in epithelial cells.

Monitoring SARS-CoV-2 Circulation and Diversity through Community Wastewater Sequencing, the Netherlands and Belgium.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly become a major global health problem, and public health surveillance is crucial to monitor and prevent virus spread. Wastewater-based epidemiology has been proposed as an addition to disease-based surveillance because virus is shed in the feces of ≈40% of infected persons. We used next-generation sequencing of sewage samples to evaluate the diversity of SARS-CoV-2 at the community level in the Netherlands and Belgium. Phylogenetic analysis revealed the presence of the most prevalent clades (19A, 20A, and 20B) and clustering of sewage samples with clinical samples from the same region. We distinguished multiple clades within a single sewage sample by using low-frequency variant analysis. In addition, several novel mutations in the SARS-CoV-2 genome were detected. Our results illustrate how wastewater can be used to investigate the diversity of SARS-CoV-2 viruses circulating in a community and identify new outbreaks.

Human Noroviruses Attach to Intestinal Tissues of a Broad Range of Animal Species.

Human noroviruses are the most common nonbacterial cause of gastroenteritis outbreaks, with new variants and genotypes frequently emerging. The origin of these new viruses is unknown; however, animals have been proposed as a potential source, as human noroviruses have been detected in animal species. Here, we investigated the potential of animals to serve as a reservoir of human noroviruses by testing norovirus attachment to formalin-fixed intestinal tissues of a range of potential reservoir animals. We set up a novel method to study norovirus binding using fluorescein isothiocyanate (FITC)-labeled virus-like particles (VLPs). In humans, noroviruses interact with histo-blood group antigens (HBGAs), carbohydrates that are expressed, among others, on the epithelial lining of the gastrointestinal tract. In animals, this interaction is not well understood. To test if virus binding depends on HBGAs, we characterized the HBGA phenotype in animal tissues by immunohistochemistry. With the exception of the black-headed gull and the straw-colored fruitbat, we observed the attachment of several human norovirus genotypes to the intestinal epithelium of all tested animal species. However, we did not find an association between the expression of a specific HBGA phenotype and virus-like particle (VLP) attachment. We show that selected human noroviruses can attach to small-intestinal tissues across species, supporting the hypothesis that human noroviruses can reside in an animal reservoir. However, whether this attachment can subsequently lead to infection needs to be further assessed.IMPORTANCE Noroviruses are a major cause of acute gastroenteritis in humans. New norovirus variants and recombinants (re)emerge regularly in the human population. From animal experiments and surveillance studies, it has become clear that at least seven animal models are susceptible to infection with human strains and that domesticated and wild animals shed human noroviruses in their feces. As virus attachment is an important first step for infection, we used a novel method utilizing FITC-labeled VLPs to test for norovirus attachment to intestinal tissues of potential animal hosts. We further characterized these tissues with regard to their HBGA expression, a well-studied norovirus susceptibility factor in humans. We found attachment of several human strains to a variety of animal species independent of their HBGA phenotype. This supports the hypothesis that human strains could reside in an animal reservoir.

Norovirus outbreak in a natural playground: A One Health approach.

Norovirus constitutes the most frequently identified infectious cause of disease outbreaks associated with untreated recreational water. When investigating outbreaks related to surface water, a One Health approach is insightful. Historically, there has been a focus on potential contamination of recreational water by bird droppings and a recent publication demonstrating human noroviruses in bird faeces suggested this should be investigated in future water-related norovirus outbreaks. Here, we describe a One Health approach investigating a norovirus outbreak in a natural playground. On social media, a large amount of waterfowl were reported to defecate near these playground premises leading to speculations about their potential involvement. Surface water, as well as human and bird faecal specimens, was tested for human noroviruses. Norovirus was found to be the most likely cause of the outbreak but there was no evidence for transmission via waterfowl. Cases had become known on social media prior to notification to the public health service underscoring the potential of online media as an early warning system. In view of known risk factors, advice was given for future outbreak investigations and natural playground design.

Comparison of sequencing methods and data processing pipelines for whole genome sequencing and minority single nucleotide variant (mSNV) analysis during an influenza A/H5N8 outbreak.

As high-throughput sequencing technologies are becoming more widely adopted for analysing pathogens in disease outbreaks there needs to be assurance that the different sequencing technologies and approaches to data analysis will yield reliable and comparable results. Conversely, understanding where agreement cannot be achieved provides insight into the limitations of these approaches and also allows efforts to be focused on areas of the process that need improvement. This manuscript describes the next-generation sequencing of three closely related viruses, each analysed using different sequencing strategies, sequencing instruments and data processing pipelines. In order to determine the comparability of consensus sequences and minority (sub-consensus) single nucleotide variant (mSNV) identification, the biological samples, the sequence data from 3 sequencing platforms and the *.bam quality-trimmed alignment files of raw data of 3 influenza A/H5N8 viruses were shared. This analysis demonstrated that variation in the final result could be attributed to all stages in the process, but the most critical were the well-known homopolymer errors introduced by 454 sequencing, and the alignment processes in the different data processing pipelines which affected the consistency of mSNV detection. However, homopolymer errors aside, there was generally a good agreement between consensus sequences that were obtained for all combinations of sequencing platforms and data processing pipelines. Nevertheless, minority variant analysis will need a different level of careful standardization and awareness about the possible limitations, as shown in this study.

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.

Complete Genome Sequences of Six Measles Virus Strains.

Genetic characterization of wild-type measles virus (MV) strains is a critical component of measles surveillance and molecular epidemiology. We have obtained complete genome sequences of six MV strains belonging to different genotypes, using random-primed next generation sequencing.

Whole-Genome Next-Generation Sequencing to Study Within-Host Evolution of Norovirus (NoV) Among Immunocompromised Patients With Chronic NoV Infection.

Background: The genus Norovirus comprises large genetic diversity, and new GII.4 variants emerge every 2-3 years. It is unknown in which host these new variants originate. Here we study whether prolonged shedders within the immunocompromised population could be a reservoir for newly emerging strains. Methods: Sixty-five fecal samples from 16 immunocompromised patients were retrospectively selected. Isolated viral RNA was enriched by hybridization with a custom norovirus whole-genome RNA bait set and deep sequenced on the Illumina MiSeq platform. Results: Patients shed virus for average 352 days (range, 76-716 days). Phylogenetic analysis showed distinct GII.4 variants in 3 of 13 patients (23%). The viral mutation rates were variable between patients but did not differ between various immune status groups. All within-host GII.4 viral populations showed amino acid changes at blocking epitopes over time, and the majority of VP1 amino acid mutations were located at the capsid surface. Conclusions: This study found viruses in immunocompromised hosts that are genetically distinct from viruses circulating in the general population, and these patients therefore may contain a reservoir for newly emerging strains. Future studies need to determine whether these new strains are of risk to other immunocompromised patients and the general population.

Detection of Circovirus in Foxes with Meningoencephalitis, United Kingdom, 2009-2013.

A fox circovirus was identified in serum samples from foxes with unexplained neurologic signs by using viral metagenomics. Fox circovirus nucleic acid was localized in histological lesions of the cerebrum by in situ hybridization. Viruses from the family Circoviridae may have neurologic tropism more commonly than previously anticipated.

Exploring the potential of next-generation sequencing in detection of respiratory viruses.

Efficient detection of human respiratory viral pathogens is crucial in the management of patients with acute respiratory tract infection. Sequence-independent amplification of nucleic acids combined with next-generation sequencing technology and bioinformatics analyses is a promising strategy for identifying pathogens in clinical and public health settings. It allows the characterization of hundreds of different known pathogens simultaneously and of novel pathogens that elude conventional testing. However, major hurdles for its routine use exist, including cost, turnaround time, and especially sensitivity of the assay, as the detection limit is dependent on viral load, host genetic material, and sequencing depth. To obtain insights into these aspects, we analyzed nasopharyngeal aspirates from a cohort of 81 Thai children with respiratory disease for the presence of respiratory viruses using a sequence-independent next-generation sequencing approach and routinely used diagnostic real-time reverse transcriptase PCR (real-time RT-PCR) assays. With respect to the detection of rhinovirus and human metapneumovirus, the next-generation sequencing approach was at least as sensitive as diagnostic real-time RT-PCR in this small cohort, whereas for bocavirus and enterovirus, next-generation sequencing was less sensitive than real-time RT-PCR. The advantage of the sequencing approach over real-time RT-PCR was the immediate availability of virus-typing information. Considering the development of platforms capable of generating more output data at declining costs, next-generation sequencing remains of interest for future virus diagnosis in clinical and public health settings and certainly as an additional tool when screening results from real-time RT-PCR are negative.

Metagenomic survey for viruses in Western Arctic caribou, Alaska, through iterative assembly of taxonomic units.

Pathogen surveillance in animals does not provide a sufficient level of vigilance because it is generally confined to surveillance of pathogens with known economic impact in domestic animals and practically nonexistent in wildlife species. As most (re-)emerging viral infections originate from animal sources, it is important to obtain insight into viral pathogens present in the wildlife reservoir from a public health perspective. When monitoring living, free-ranging wildlife for viruses, sample collection can be challenging and availability of nucleic acids isolated from samples is often limited. The development of viral metagenomics platforms allows a more comprehensive inventory of viruses present in wildlife. We report a metagenomic viral survey of the Western Arctic herd of barren ground caribou (Rangifer tarandus granti) in Alaska, USA. The presence of mammalian viruses in eye and nose swabs of 39 free-ranging caribou was investigated by random amplification combined with a metagenomic analysis approach that applied exhaustive iterative assembly of sequencing results to define taxonomic units of each metagenome. Through homology search methods we identified the presence of several mammalian viruses, including different papillomaviruses, a novel parvovirus, polyomavirus, and a virus that potentially represents a member of a novel genus in the family Coronaviridae.

Viral metagenomic analysis of feces of wild small carnivores.

BACKGROUND: Recent studies have clearly demonstrated the enormous virus diversity that exists among wild animals. This exemplifies the required expansion of our knowledge of the virus diversity present in wildlife, as well as the potential transmission of these viruses to domestic animals or humans. METHODS: In the present study we evaluated the viral diversity of fecal samples (n = 42) collected from 10 different species of wild small carnivores inhabiting the northern part of Spain using random PCR in combination with next-generation sequencing. Samples were collected from American mink (Neovison vison), European mink (Mustela lutreola), European polecat (Mustela putorius), European pine marten (Martes martes), stone marten (Martes foina), Eurasian otter (Lutra lutra) and Eurasian badger (Meles meles) of the family of Mustelidae; common genet (Genetta genetta) of the family of Viverridae; red fox (Vulpes vulpes) of the family of Canidae and European wild cat (Felis silvestris) of the family of Felidae. RESULTS: A number of sequences of possible novel viruses or virus variants were detected, including a theilovirus, phleboviruses, an amdovirus, a kobuvirus and picobirnaviruses. CONCLUSIONS: Using random PCR in combination with next generation sequencing, sequences of various novel viruses or virus variants were detected in fecal samples collected from Spanish carnivores. Detected novel viruses highlight the viral diversity that is present in fecal material of wild carnivores.

New viruses in idiopathic human diarrhea cases, the Netherlands.

Emerging viral infections can be identified by using a viral metagenomics approach for clinical human material. Diarrhea samples of patients with unexplained gastroenteritis from the Netherlands were analyzed by using viral metagenomics. Novel circular DNA viruses, bufaviruses, and genogroup III picobirnaviruses were identified. These data expand our knowledge of the human virome.

Novel B19-like parvovirus in the brain of a harbor seal.

Using random PCR in combination with next-generation sequencing, a novel parvovirus was detected in the brain of a young harbor seal (Phoca vitulina) with chronic non-suppurative meningo-encephalitis that was rehabilitated at the Seal Rehabilitation and Research Centre (SRRC) in the Netherlands. In addition, two novel viruses belonging to the family Anelloviridae were detected in the lungs of this animal. Phylogenetic analysis of the coding sequence of the novel parvovirus, tentatively called Seal parvovirus, indicated that this virus belonged to the genus Erythrovirus, to which human parvovirus B19 also belongs. Although no other seals with similar signs were rehabilitated in SRRC in recent years, a prevalence study of tissues of seals from the same area collected in the period 2008-2012 indicated that the Seal parvovirus has circulated in the harbor seal population at least since 2008. The presence of the Seal parvovirus in the brain was confirmed by real-time PCR and in vitro replication. Using in situ hybridization, we showed for the first time that a parvovirus of the genus Erythrovirus was present in the Virchow-Robin space and in cerebral parenchyma adjacent to the meninges. These findings showed that a parvovirus of the genus Erythrovirus can be involved in central nervous system infection and inflammation, as has also been suspected but not proven for human parvovirus B19 infection.

Novel cyclovirus in human cerebrospinal fluid, Malawi, 2010-2011.

To identify unknown human viruses, we analyzed serum and cerebrospinal fluid samples from patients with unexplained paraplegia from Malawi by using viral metagenomics. A novel cyclovirus species was identified and subsequently found in 15% and 10% of serum and cerebrospinal fluid samples, respectively. These data expand our knowledge of cyclovirus diversity and tropism.