Metagenomic Surveillance of Viral Gastroenteritis in a Public Health Setting.

Norovirus is the primary cause of viral gastroenteritis (GE). To investigate norovirus epidemiology, there is a need for whole-genome sequencing and reference sets consisting of complete genomes. To investigate the potential of shotgun metagenomic sequencing on the Illumina platform for whole-genome sequencing, 71 reverse transcriptase quantitative PCR (RT-qPCR) norovirus positive-feces (threshold cycle [CT], <30) samples from norovirus surveillance within The Netherlands were subjected to metagenomic sequencing. Data were analyzed through an in-house next-generation sequencing (NGS) analysis workflow. Additionally, we assessed the potential of metagenomic sequencing for the surveillance of off-target viruses that are of importance for public health, e.g., sapovirus, rotavirus A, enterovirus, parechovirus, aichivirus, adenovirus, and bocaparvovirus. A total of 60 complete and 10 partial norovirus genomes were generated, representing 7 genogroup I capsid genotypes and 12 genogroup II capsid genotypes. In addition to the norovirus genomes, the metagenomic approach yielded partial or complete genomes of other viruses for 39% of samples from children and 6.7% of samples from adults, including adenovirus 41 (N = 1); aichivirus 1 (N = 1); coxsackievirus A2 (N = 2), A4 (N = 2), A5 (N = 1), and A16 (N = 1); bocaparvovirus 1 (N = 1) and 3 (N = 1); human parechovirus 1 (N = 2) and 3 (N = 1); Rotavirus A (N = 1); and a sapovirus GI.7 (N = 1). The sapovirus GI.7 was initially not detected through RT-qPCR and warranted an update of the primer and probe set. Metagenomic sequencing on the Illumina platform robustly determines complete norovirus genomes and may be used to broaden gastroenteritis surveillance by capturing off-target enteric viruses. IMPORTANCE Viral gastroenteritis results in significant morbidity and mortality in vulnerable individuals and is primarily caused by norovirus. To investigate norovirus epidemiology, there is a need for whole-genome sequencing and reference sets consisting of full genomes. Using surveillance samples sent to the Dutch National Institute for Public Health and the Environment (RIVM), we compared metagenomics against conventional techniques, such as RT-qPCR and Sanger-sequencing, with norovirus as the target pathogen. We determined that metagenomics is a robust method to generate complete norovirus genomes, in parallel to many off-target pathogenic enteric virus genomes, thereby broadening our surveillance efforts. Moreover, we detected a sapovirus that was not detected by our validated gastroenteritis RT-qPCR panel, which exemplifies the strength of metagenomics. Our study shows that metagenomics can be used for public health gastroenteritis surveillance, the generation of reference-sets for molecular epidemiology, and how it compares to current surveillance strategies.

Elevated risk of infection with SARS-CoV-2 Beta, Gamma, and Delta variants compared with Alpha variant in vaccinated individuals.

The extent to which severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) break through infection- or vaccine-induced immunity is not well understood. We analyzed 28,578 sequenced SARS-CoV-2 samples from individuals with known immune status obtained through national community testing in the Netherlands from March to August 2021. We found evidence of an increased risk of infection by the Beta (B.1.351), Gamma (P.1), or Delta (B.1.617.2) variants compared with the Alpha (B.1.1.7) variant after vaccination. No clear differences were found between vaccines. However, the effect was larger in the first 14 to 59 days after complete vaccination compared with ≥60 days. In contrast to vaccine-induced immunity, there was no increased risk for reinfection with Beta, Gamma, or Delta variants relative to the Alpha variant in individuals with infection-induced immunity.

Molecular Epidemiology and Evolutionary Trajectory of Emerging Echovirus 30, Europe.

In 2018, an upsurge in echovirus 30 (E30) infections was reported in Europe. We conducted a large-scale epidemiologic and evolutionary study of 1,329 E30 strains collected in 22 countries in Europe during 2016-2018. Most E30 cases affected persons 0-4 years of age (29%) and 25-34 years of age (27%). Sequences were divided into 6 genetic clades (G1-G6). Most (53%) sequences belonged to G1, followed by G6 (23%), G2 (17%), G4 (4%), G3 (0.3%), and G5 (0.2%). Each clade encompassed unique individual recombinant forms; G1 and G4 displayed >2 unique recombinant forms. Rapid turnover of new clades and recombinant forms occurred over time. Clades G1 and G6 dominated in 2018, suggesting the E30 upsurge was caused by emergence of 2 distinct clades circulating in Europe. Investigation into the mechanisms behind the rapid turnover of E30 is crucial for clarifying the epidemiology and evolution of these enterovirus infections.

Recommendations for the introduction of metagenomic next-generation sequencing in clinical virology, part II: bioinformatic analysis and reporting.

Metagenomic next-generation sequencing (mNGS) is an untargeted technique for determination of microbial DNA/RNA sequences in a variety of sample types from patients with infectious syndromes. mNGS is still in its early stages of broader translation into clinical applications. To further support the development, implementation, optimization and standardization of mNGS procedures for virus diagnostics, the European Society for Clinical Virology (ESCV) Network on Next-Generation Sequencing (ENNGS) has been established. The aim of ENNGS is to bring together professionals involved in mNGS for viral diagnostics to share methodologies and experiences, and to develop application guidelines. Following the ENNGS publication Recommendations for the introduction of mNGS in clinical virology, part I: wet lab procedure in this journal, the current manuscript aims to provide practical recommendations for the bioinformatic analysis of mNGS data and reporting of results to clinicians.

Molecular epidemiology of mumps viruses in the Netherlands, 2017-2019.

Mumps cases continue to occur, also in countries with a relatively high vaccination rate. The last major outbreaks of mumps in the Netherlands were in 2009-2012 and thereafter, only small clusters and single cases were reported. Molecular epidemiology can provide insights in the circulation of mumps viruses. The aims of the present study were to analyze the molecular epidemiology of mumps viruses in the Netherlands in 2017-2019 and to compare the phylogenetic trees built from sequence data of near complete mumps virus genomes or from the SH gene and non-coding regions (SH+NCRs). To this end, Sanger sequence data from SH+NCRs were analyzed from 82 mumps genotype G viruses. In addition, near complete genomes were obtained from 10 mumps virus isolates using next-generation sequencing. Analysis of SH+NCRs sequences of mumps genotype G viruses revealed the presence of two major genetic lineages in the Netherlands, which was confirmed by analysis of near complete genomes. Comparison of phylogenetic trees built with SH+NCRs or near complete genomes indicated that the topology was similar, while somewhat longer branches were present in the phylogenetic tree with near complete genomes. These results confirm that analysis of SH + NCRs sequence data is a useful approach for molecular surveillance. Furthermore, data from recent mumps genotype G viruses might indicate (intermittent) circulation of mumps genotype G viruses in the Netherlands in 2017-2019.

The COMPARE Data Hubs.

Data sharing enables research communities to exchange findings and build upon the knowledge that arises from their discoveries. Areas of public and animal health as well as food safety would benefit from rapid data sharing when it comes to emergencies. However, ethical, regulatory and institutional challenges, as well as lack of suitable platforms which provide an infrastructure for data sharing in structured formats, often lead to data not being shared or at most shared in form of supplementary materials in journal publications. Here, we describe an informatics platform that includes workflows for structured data storage, managing and pre-publication sharing of pathogen sequencing data and its analysis interpretations with relevant stakeholders.

Proficiency Testing of Virus Diagnostics Based on Bioinformatics Analysis of Simulated In Silico High-Throughput Sequencing Data Sets.

Quality management and independent assessment of high-throughput sequencing-based virus diagnostics have not yet been established as a mandatory approach for ensuring comparable results. The sensitivity and specificity of viral high-throughput sequence data analysis are highly affected by bioinformatics processing using publicly available and custom tools and databases and thus differ widely between individuals and institutions. Here we present the results of the COMPARE [Collaborative Management Platform for Detection and Analyses of (Re-)emerging and Foodborne Outbreaks in Europe] in silico virus proficiency test. An artificial, simulated in silico data set of Illumina HiSeq sequences was provided to 13 different European institutes for bioinformatics analysis to identify viral pathogens in high-throughput sequence data. Comparison of the participants' analyses shows that the use of different tools, programs, and databases for bioinformatics analyses can impact the correct identification of viral sequences from a simple data set. The identification of slightly mutated and highly divergent virus genomes has been shown to be most challenging. Furthermore, the interpretation of the results, together with a fictitious case report, by the participants showed that in addition to the bioinformatics analysis, the virological evaluation of the results can be important in clinical settings. External quality assessment and proficiency testing should become an important part of validating high-throughput sequencing-based virus diagnostics and could improve the harmonization, comparability, and reproducibility of results. There is a need for the establishment of international proficiency testing, like that established for conventional laboratory tests such as PCR, for bioinformatics pipelines and the interpretation of such results.

Overview of Virus Metagenomic Classification Methods and Their Biological Applications.

Metagenomics poses opportunities for clinical and public health virology applications by offering a way to assess complete taxonomic composition of a clinical sample in an unbiased way. However, the techniques required are complicated and analysis standards have yet to develop. This, together with the wealth of different tools and workflows that have been proposed, poses a barrier for new users. We evaluated 49 published computational classification workflows for virus metagenomics in a literature review. To this end, we described the methods of existing workflows by breaking them up into five general steps and assessed their ease-of-use and validation experiments. Performance scores of previous benchmarks were summarized and correlations between methods and performance were investigated. We indicate the potential suitability of the different workflows for (1) time-constrained diagnostics, (2) surveillance and outbreak source tracing, (3) detection of remote homologies (discovery), and (4) biodiversity studies. We provide two decision trees for virologists to help select a workflow for medical or biodiversity studies, as well as directions for future developments in clinical viral metagenomics.

The correlation between severity of paraparesis and reduced density of resident antigen-presenting cells implicates an unknown role for the spinal perivascular macrophages in experimental autoimmune encephalomyelitis in rats.

To study alterations in the morphology of spinal perivascular macrophages (SPM) during experimental allergic encephalomyelitis (EAE), we labelled SPM by intracerebroventricular (i.c.v.) injection of horseradish peroxidase (HRP). As earlier electron microscopical analysis had shown severely damaged SPM, we suspected that each inflammatory process is accompanied by the death of SPM. To prove this hypothesis, we compared the numerical density of resident SPM (i.c.v. labelled in red by Fluoro-Ruby) with that of monocytes/macrophages recruited to the perivascular space (i.c.v. labelled in green by Fluoro-Emerald). At the peak of paraparesis, the density of resident SPM was reduced by 33%. Since this reduction contrasted sharply with earlier data indicating a massive increase in the density of SPM during EAE, we checked our findings after general or selective suppression of the immune response to myelin autoantigens with the drugs dexamethasone and copaxone, respectively. Dexamethasone treatment commenced after evident paraparesis accelerated recovery, but did not influence SPM density. Immunisation with copaxone completely prevented the occurrence of EAE (monitored by video-based motion analysis of tail motility); the subsequent histological analysis revealed no reduction in SPM density. Based on this inverse correlation between the severity of EAE and the density of resident macrophages, we conclude that SPM plays an important role in the pathogenesis of EAE.