Search

Objectives: To obtain insight into SARS-CoV-2 clustering and transmission routes during outbreaks in the predominantly migrant workforce of the fruit and vegetable packaging industry of South Holland, the Netherlands, May to July 2020. Design: This mixed-methods study applied direct observation and interviews, epidemiologic investigation, source and contact data analysis and whole-genome sequencing. Results: We detected 46 SARS-CoV-2 cases and 4 outbreaks with a proportional representation of labour migrant and native workers in 6 unrelated facilities. Complete viral genome sequences revealed at least 3 clusters of native workers and labour migrants, 2 within and 1 between facilities. On-site inspections found adequate implementation of preventative measures to which both native workers and labour migrants showed suboptimal adherence. Being a labour migrant was associated with living in shared housing, but not with more contacts or different sources. Conclusions: The fruit and vegetable packaging industry gave the impression of sufficient preparedness and control. Suboptimal adherence to the facilities’ preventative guidelines could have facilitated work floor transmission. Community and household transmission are likely to have contributed to outbreaks. We encourage further research into risk factors for transmission in labour migrants and application of these insights into targeted public health policy.

After the first wave of SARS-CoV-2 infections in spring 2020, Europe experienced a resurgence of the virus starting in late summer 2020 that was deadlier and more difficult to contain1. Relaxed intervention measures and summer travel have been implicated as drivers of the second wave2. Here we build a phylogeographical model to evaluate how newly introduced lineages, as opposed to the rekindling of persistent lineages, contributed to the resurgence of COVID-19 in Europe. We inform this model using genomic, mobility and epidemiological data from 10 European countries and estimate that in many countries more than half of the lineages circulating in late summer resulted from new introductions since 15 June 2020. The success in onward transmission of newly introduced lineages was negatively associated with the local incidence of COVID-19 during this period. The pervasive spread of variants in summer 2020 highlights the threat of viral dissemination when restrictions are lifted, and this needs to be carefully considered in strategies to control the current spread of variants that are more transmissible and/or evade immunity. Our findings indicate that more effective and coordinated measures are required to contain the spread through cross-border travel even as vaccination is reducing disease burden. ispartof: NATURE vol:595 issue:7869 pages:713-+ ispartof: location:England status: published

The emergence of SARS-CoV-2 variants harboring mutations in the spike (S) protein has raised concern about potential immune escape. Here, we studied humoral and cellular immune responses to wild type SARS-CoV-2 and the B.1.1.7 and B.1.351 variants of concern in a cohort of 121 BNT162b2 mRNA-vaccinated health care workers (HCW). Twenty-three HCW recovered from mild COVID-19 disease and exhibited a recall response with high levels of SARS-CoV-2-specific functional antibodies and virus-specific T cells after a single vaccination. Specific immune responses were also detected in seronegative HCW after one vaccination, but a second dose was required to reach high levels of functional antibodies and cellular immune responses in all individuals. Vaccination-induced antibodies cross-neutralized the variants B.1.1.7 and B.1.351, but the neutralizing capacity and Fc-mediated functionality against B.1.351 was consistently 2-to 4-fold lower than to the homologous virus. In addition, peripheral blood mononuclear cells were stimulated with peptide pools spanning the mutated S regions of B.1.1.7 and B.1.351 to detect cross-reactivity of SARS-CoV-2-specific T cells with variants. Importantly, we observed no differences in CD4+ T-cell activation in response to variant antigens, indicating that the B.1.1.7 and B.1.351 S proteins do not escape T-cell-mediated immunity elicited by the wild type S protein. In conclusion, this study shows that some variants can partially escape humoral immunity induced by SARS-CoV-2 infection or BNT162b2 vaccination, but S-specific CD4+ T-cell activation is not affected by the mutations in the B.1.1.7 and B.1.351 variants.

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501YV2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected., Wellcome Open Research, 6

Metagenomic sequencing is increasingly being used in clinical settings for difficult to diagnose cases. The performance of viral metagenomic protocols relies to a large extent on the bioinformatic analysis. In this study, the European Society for Clinical Virology (ESCV) Network on NGS (ENNGS) initiated a benchmark of metagenomic pipelines currently used in clinical virological laboratories.MethodsMetagenomic datasets from 13 clinical samples from patients with encephalitis or viral respiratory infections characterized by PCR were selected. The datasets were analysed with 13 different pipelines currently used in virological diagnostic laboratories of participating ENNGS members. The pipelines and classification tools were: Centrifuge, DAMIAN, DIAMOND, DNASTAR, FEVIR, Genome Detective, Jovian, MetaMIC, MetaMix, One Codex, RIEMS, VirMet, and Taxonomer. Performance, characteristics, clinical use, and user-friendliness of these pipelines were analysed.ResultsOverall, viral pathogens with high loads were detected by all the evaluated metagenomic pipelines. In contrast, lower abundance pathogens and mixed infections were only detected by 3/13 pipelines, namely DNASTAR, FEVIR, and MetaMix. Overall sensitivity ranged from 80% (10/13) to 100% (13/13 datasets). Overall positive predictive value ranged from 71-100%. The majority of the pipelines classified sequences based on nucleotide similarity (8/13), only a minority used amino acid similarity, and 6 of the 13 pipelines assembled sequences de novo. No clear differences in performance were detected that correlated with these classification approaches. Read counts of target viruses varied between the pipelines over a range of 2-3 log, indicating differences in limit of detection.ConclusionA wide variety of viral metagenomic pipelines is currently used in the participating clinical diagnostic laboratories. Detection of low abundant viral pathogens and mixed infections remains a challenge, implicating the need for standardization and validation of metagenomic analysis for clinical diagnostic use. Future studies should address the selective effects due to the choice of different reference viral databases.

Experiments in which complex virome sequencing data is generated remain difficult to explore and unpack for scientists without a background in data science. The processing of raw sequencing data by high throughput sequencing workflows usually results in contigs in FASTA format coupled to an annotation file linking the contigs to a reference sequence or taxonomic identifier. The next step is to compare the virome of different samples based on the metadata of the experimental setup and extract sequences of interest that can be used in subsequent analyses. The viromeBrowser is an application written in the opensource R shiny framework that was developed in collaboration with end-users and is focused on three common data analysis steps. First, the application allows interactive filtering of annotations by default or custom quality thresholds. Next, multiple samples can be visualized to facilitate comparison of contig annotations based on sample specific metadata values. Last, the application makes it easy for users to extract sequences of interest in FASTA format. With the interactive features in the viromeBrowser we aim to enable scientists without a data science background to compare and extract annotation data and sequences from virome sequencing analysis results.

Following the first wave of SARS-CoV-2 infections in spring 2020, Europe experienced a resurgence of the virus starting late summer that was deadlier and more difficult to contain. Relaxed intervention measures and summer travel have been implicated as drivers of the second wave. Here, we build a phylogeographic model to evaluate how newly introduced lineages, as opposed to the rekindling of persistent lineages, contributed to the COVID-19 resurgence in Europe. We inform this model using genomic, mobility and epidemiological data from 10 West European countries and estimate that in many countries more than 50% of the lineages circulating in late summer resulted from new introductions since June 15th. The success in onwards transmission of these lineages is predicted by SARS-CoV-2 incidence during this period. Relatively early introductions from Spain into the United Kingdom contributed to the successful spread of the 20A.EU1/B.1.177 variant. The pervasive spread of variants that have not been associated with an advantage in transmissibility highlights the threat of novel variants of concern that emerged more recently and have been disseminated by holiday travel. Our findings indicate that more effective and coordinated measures are required to contain spread through cross-border travel.

The current SARS-CoV-2 pandemic has rapidly become a major global health problem for which public health surveillance is crucial to monitor virus spread. Given the presence of viral RNA in feces in around 40% of infected persons, wastewater-based epidemiology has been proposed as an addition to disease-based surveillance to assess the spread of the virus at the community level. Here we have explored the possibility of using next-generation sequencing (NGS) of sewage samples to evaluate the diversity of SARS-CoV-2 at the community level from routine wastewater testing, and compared these results with the virus diversity in patients from the Netherlands and Belgium. Phylogenetic analysis revealed the presence of viruses belonging to the most prevalent clades (19A, 20A and 20B) in both countries. Clades 19B and 20C were not identified, while they were present in clinical samples during the same period. Low frequency variant (LFV) analysis showed that some known LFVs can be associated with particular clusters within a clade, different to those of their consensus sequences, suggesting the presence of at least 2 clades within a single sewage sample. Additionally, combining genome consensus and LFV analyses we found a total of 57 unique mutations in the SARS-CoV-2 genome which have not been described before. In conclusion, this work illustrates how NGS analysis of wastewater can be used to approximate the diversity of SARS-CoV-2 viruses circulating in a community.

INTRODUCTION: Metagenomic sequencing is increasingly being used in clinical settings for difficult to diagnose cases. The performance of viral metagenomic protocols relies to a large extent on the bioinformatic analysis. In this study, the European Society for Clinical Virology (ESCV) Network on NGS (ENNGS) initiated a benchmark of metagenomic pipelines currently used in clinical virological laboratories. METHODS: Metagenomic datasets from 13 clinical samples from patients with encephalitis or viral respiratory infections characterized by PCR were selected. The datasets were analyzed with 13 different pipelines currently used in virological diagnostic laboratories of participating ENNGS members. The pipelines and classification tools were: Centrifuge, DAMIAN, DIAMOND, DNASTAR, FEVIR, Genome Detective, Jovian, MetaMIC, MetaMix, One Codex, RIEMS, VirMet, and Taxonomer. Performance, characteristics, clinical use, and user-friendliness of these pipelines were analyzed. RESULTS: Overall, viral pathogens with high loads were detected by all the evaluated metagenomic pipelines. In contrast, lower abundance pathogens and mixed infections were only detected by 3/13 pipelines, namely DNASTAR, FEVIR, and MetaMix. Overall sensitivity ranged from 80% (10/13) to 100% (13/13 datasets). Overall positive predictive value ranged from 71-100%. The majority of the pipelines classified sequences based on nucleotide similarity (8/13), only a minority used amino acid similarity, and 6 of the 13 pipelines assembled sequences de novo. No clear differences in performance were detected that correlated with these classification approaches. Read counts of target viruses varied between the pipelines over a range of 2-3 log, indicating differences in limit of detection. CONCLUSION: A wide variety of viral metagenomic pipelines is currently used in the participating clinical diagnostic laboratories. Detection of low abundant viral pathogens and mixed infections remains a challenge, implicating the need for standardization and validation of metagenomic analysis for clinical diagnostic use. Future studies should address the selective effects due to the choice of different reference viral databases. ispartof: JOURNAL OF CLINICAL VIROLOGY vol:141 ispartof: location:Netherlands status: published

Background: 10 days after the first reported case of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in the Netherlands (on Feb 27, 2020), 55 (4%) of 1497 health-care workers in nine hospitals located in the south of the Netherlands had tested positive for SARS-CoV-2 RNA. We aimed to gain insight in possible sources of infection in health-care workers. Methods: We did a cross-sectional study at three of the nine hospitals located in the south of the Netherlands. We screened health-care workers at the participating hospitals for SARS-CoV-2 infection, based on clinical symptoms (fever or mild respiratory symptoms) in the 10 days before screening. We obtained epidemiological data through structured interviews with health-care workers and combined this information with data from whole-genome sequencing of SARS-CoV-2 in clinical samples taken from health-care workers and patients. We did an in-depth analysis of sources and modes of transmission of SARS-CoV-2 in health-care workers and patients. Findings: Between March 2 and March 12, 2020, 1796 (15%) of 12 022 health-care workers were screened, of whom 96 (5%) tested positive for SARS-CoV-2. We obtained complete and near-complete genome sequences from 50 health-care workers and ten patients. Most sequences were grouped in three clusters, with two clusters showing local circulation within the region. The noted patterns were consistent with multiple introductions into the hospitals through community-acquired infections and local amplification in the community. Interpretation: Although direct transmission in the hospitals cannot be ruled out, our data do not support widespread nosocomial transmission as the source of infection in patients or health-care workers. Funding: EU Horizon 2020 (RECoVer, VEO, and the European Joint Programme One Health METASTAVA), and the National Institute of Allergy and Infectious Diseases, National Institutes of Health.

Here we describe a SARS-CoV-2 variant with diminished amplification of the ORF ORF1ab target in the Cobas® dual-target SARS-CoV-2 assay resulting in a discrepancy of Ct-values (Ct-value 20.7 for the E-gene and Ct-value 30.2 for ORF1ab). Five unique nucleotide mutations were identified in ORF1ab: C11450A (nsp10) C14178T (RdRp), G15006T (RdRp), G18394T (Hel), and G20995T (Hel). This case highlights the importance of surveillance of genomic regions used in molecular diagnostics and the importance of the public release of target regions used to update commercial and in-house developed SARS-CoV-2 PCR tests. This work underpins the importance of using dual-targets in molecular diagnostic assays to limit the change of false-negative results due to primer and/or probe mismatches.

The current SARS-CoV-2 pandemic has rapidly become a major global health problem for which public health surveillance is crucial to monitor virus spread. Given the presence of viral RNA in feces in around 40% of infected persons, wastewater-based epidemiology has been proposed as an addition to disease-based surveillance to assess the spread of the virus at the community level. Here we have explored the possibility of using next-generation sequencing (NGS) of sewage samples to evaluate the diversity of SARS-CoV-2 at the community level from routine wastewater testing, and compared these results with the virus diversity in patients from the Netherlands and Belgium. Phylogenetic analysis revealed the presence of viruses belonging to the most prevalent clades (19A, 20A and 20B) in both countries. Clades 19B and 20C were not identified, while they were present in clinical samples during the same period. Low frequency variant (LFV) analysis showed that some known LFVs can be associated with particular clusters within a clade, different to those of their consensus sequences, suggesting the presence of at least 2 clades within a single sewage sample. Additionally, combining genome consensus and LFV analyses we found a total of 57 unique mutations in the SARS-CoV-2 genome which have not been described before. In conclusion, this work illustrates how NGS analysis of wastewater can be used to approximate the diversity of SARS-CoV-2 viruses circulating in a community.

BackgroundThis study aimed to assess the contribution of asymptomatic and presymptomatic residents and staff in SARS-CoV-2 transmission during a large outbreak in a Dutch nursing home.MethodsObservational study in a 185-bed nursing home with two consecutive testing strategies: testing of symptomatic cases only, and weekly facility-wide testing of staff and residents regardless of symptoms. Nasopharyngeal and oropharyngeal testing with RT-PCR for SARs-CoV-2 was conducted with a standardized symptom assessment. Positive samples with a cycle threshold (CT) value below 32 were selected for sequencing.Results185 residents and 244 staff participated. Sequencing identified one cluster. In the symptom-based test strategy period 3/39 residents were presymptomatic versus 38/74 residents in the period of weekly facility-wide testing (p-valueConclusionsThe frequency of a/presymptomatic residents compared to staff suggests that a/presymptomatic residents could be unrecognized symptomatic cases. However, symptomatic and presymptomatic/unrecognized symptomatic residents both have the same potential for viral shedding. The high prevalence symptomatic staff found in facility-wide testing suggests that staff has difficulty attributing their symptoms to possible SARS-CoV-2 infection. Weekly testing was an effective strategy for early identification of SARS-Cov-2 cases, resulting in fast isolation and mitigation of this outbreak.

The zoonotic origin of the SARS-CoV-2 pandemic is still unknown. Animal experiments have shown that non-human primates, cats, ferrets, hamsters, rabbits and bats can be infected by SARS-CoV-2. In addition, SARS-CoV-2 RNA has been detected in felids, mink and dogs in the field. Here, we describe an in-depth investigation of outbreaks on 16 mink farms and humans living or working on these farms, using whole genome sequencing. We conclude that the virus was initially introduced from humans and has evolved, most likely reflecting widespread circulation among mink in the beginning of the infection period several weeks prior to detection. At the moment, despite enhanced biosecurity, early warning surveillance and immediate culling of infected farms, there is ongoing transmission between mink farms with three big transmission clusters with unknown modes of transmission. We also describe the first animal to human transmissions of SARS-CoV-2 in mink farms.One sentence summarySARS-CoV-2 transmission on mink farms.

Background The Dutch province of Limburg borders the German district of Heinsberg, which had a large cluster of COVID-19 cases linked to local carnival activities before any cases were reported in the Netherlands. However, Heinsberg was not included as an area reporting local or community transmission per the national case definition at the time. In early March, two residents from a long-term care facility (LTCF) in Sittard, a Dutch town located in close vicinity to the district of Heinsberg, tested positive for COVID-19. In this study we aimed to determine whether cross-border introduction of the virus took place by analysing the LTCF outbreak in Sittard, both epidemiologically and microbiologically. Methods Surveys and semi-structured oral interviews were conducted with all present LTCF residents by health care workers during regular points of care for information on new or unusual signs and symptoms of disease. Both throat and nasopharyngeal swabs were taken from residents suspect of COVID-19, based on regional criteria, for the detection of SARS-CoV-2 by Real-time Polymerase Chain Reaction. Additionally, whole genome sequencing was performed using a SARS-CoV-2 specific amplicon-based Nanopore sequencing approach. Moreover, twelve random residents were sampled for possible asymptomatic infections. Results Out of 99 residents, 46 got tested for COVID-19. Out of the 46 tested residents, nineteen (41%) tested positive for COVID-19, including 3 asymptomatic residents. CT-values for asymptomatic residents seemed higher compared to symptomatic residents. Eleven samples were sequenced, along with three random samples from COVID-19 patients hospitalized in the regional hospital at the time of the LTCF outbreak. All samples were linked to COVID-19 cases from the cross-border region of Heinsberg, Germany. Conclusions Sequencing combined with epidemiological data was able to virtually prove cross-border transmission at the start of the Dutch COVID-19 epidemic. Our results highlight the need for cross-border collaboration and adjustment of national policy to emerging region-specific needs along borders in order to establish coordinated implementation of infection control measures to limit the spread of COVID-19.

In April 2020, respiratory disease and increased mortality were observed in farmed mink on two farms in the Netherlands. In both farms, at least one worker had been found positive for SARS-CoV-2. Necropsies of the mink revealed interstitial pneumonia, and organ and swab samples tested positive for SARS-CoV-2 RNA by qPCR. Variations in viral genomes point at between-mink transmission on the farms and lack of infection link between the farms. Inhalable dust in the mink houses contained viral RNA, indicating possible exposure of workers.

In April 2019, a suspect cluster of enterovirus cases was reported in a neonatology department in Guangdong, China, resulting in five deaths. We aimed to investigate the pathogen profiles in fatal cases, the circulation and transmission pattern of the viruses by combining metatranscriptomic, phylogenetic, and epidemiological analyses. Metatranscriptomic sequencing was used to characterize the enteroviruses. Clinical and environmental surveillance in the local population was performed to understand the prevalence and genetic diversity of the viruses in the local population. The possible source(s), evolution, transmission, and recombination of the viruses were investigated by incorporating genomes from the current outbreak, from local retrospective surveillance, and from public databases. Metatranscriptomic analysis identified Echovirus 11 (E11) in three fatal cases. Seroprevalence of neutralization antibody to E11 was 35 to 44 per cent in 3–15 age groups of general population, and the viruses were associated with various clinical symptoms. From the viral phylogeny, nosocomial transmissions were identified and all E11 2019 outbreak strains were closely related with E11 strains circulating in local population 2017–19. Frequent recombination occurred among the 2019 Guangdong E11 outbreak strains and various genotypes in enterovirus B species. This study provides an example of combining advanced genetic technology and epidemiological surveillance in pathogen diagnosis, source(s), and transmission tracing during an infectious disease outbreak. The result highlights the hidden E11 circulation and the risk of viral transmission and infection in the young age population in China. Frequent recombination between Guangdong-like strains and other enterovirus genotypes also implies the prevalence of these emerging E11 strains.

Whole genome sequencing can be used to characterize and to trace viral outbreaks. Nanopore-based whole genome sequencing protocols have been described for several different viruses. These approaches utilize an overlapping amplicon-based approach which can be used to target a specific virus or group of genetically related viruses. In addition to confirmation of the virus presence, sequencing can be used for genomic epidemiology studies, to track viruses and unravel origins, reservoirs and modes of transmission. For such applications, it is crucial to understand possible effects of the error rate associated with the platform used. Routine application in clinical and public health settings require that this is documented with every important change in the protocol. Previously, a protocol for whole genome Usutu virus sequencing on the nanopore sequencing platform was validated (R9.4 flowcell) by direct comparison to Illumina sequencing. Here, we describe the method used to determine the required read coverage, using the comparison between the R10 flow cell and Illumina sequencing as an example.

On 22 August, a common whitethroat in the Netherlands tested positive for West Nile virus lineage 2. The same bird had tested negative in spring. Subsequent testing of Culex mosquitoes collected in August and early September in the same location generated two of 44 positive mosquito pools, providing first evidence for enzootic transmission in the Netherlands. Sequences generated from the positive mosquito pools clustered with sequences that originate from Germany, Austria and the Czech Republic. © 2020 European Centre for Disease Prevention and Control (ECDC). All rights reserved.

The first highly pathogenic (HP) influenza A/H7N9 was reported in Guangdong in January 2017. To investigate the emergence and spread of HP A/H7N9 in Guangdong province, we sequenced 297 viruses (58 HP A/H7N9, 19 low pathogenic (LP) A/H7N9, and 220 A/H9N2) during 2016–2017. Our analysis showed that during the fifth wave, three A/H7N9 lineages were co-circulating in Guangdong: the local LP Pearl River Delta (PRD) lineage (13%), the newly imported LP Yangtze River Delta (YRD) lineage (23%), and the HP YRD lineage (64%). Previously circulating YRD-lineage LP during the third wave evolved to the YRD-lineage HP A/H7N9 in Guangdong. All YRD-lineage LP detected during the fifth wave most likely originated from newly imported viruses into Guangdong. Genotype comparison of HP A/H7N9 suggests limited outward spread of HP A/H7N9 to other provinces. The distribution of HP A/H7N9 cleavage site variants on live poultry markets differed from that found in humans, suggesting a V1-type cleavage site may facilitate human infections.

The current coronavirus disease 2019 (COVID-19) pandemic is the first to apply whole-genome sequencing near to real time, with over 2 million severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) whole-genome sequences generated and shared through the GISAID platform. This genomic resource informed public health decision-making throughout the pandemic; it also allowed detection of mutations that might affect virulence, pathogenesis, host range or immune escape as well as the effectiveness of SARS-CoV-2 diagnostics and therapeutics. However, genotype-to-phenotype predictions cannot be performed at the rapid pace of genomic sequencing. To prepare for the next phase of the pandemic, a systematic approach is needed to link global genomic surveillance and timely assessment of the phenotypic characteristics of novel variants, which will support the development and updating of diagnostics, vaccines, therapeutics and nonpharmaceutical interventions. This Review summarizes the current knowledge on key viral mutations and variants and looks to the next phase of surveillance of the evolving pandemic.

Sewage samples have a high potential benefit for surveillance of circulating pathogens because they are easy to obtain and reflect population-wide circulation of pathogens. These type of samples typically contain a great diversity of viruses. Therefore, one of the main challenges of metagenomic sequencing of sewage for surveillance is sequence annotation and interpretation. Especially for high-threat viruses, false positive signals can trigger unnecessary alerts, but true positives should not be missed. Annotation thus requires high sensitivity and specificity. To better interpret annotated reads for high-threat viruses, we attempt to determine how classifiable they are in a background of reads of closely related low-threat viruses. As an example, we attempted to distinguish poliovirus reads, a virus of high public health importance, from other enterovirus reads. A sequence-based deep learning algorithm was used to classify reads as either polio or non-polio enterovirus. Short reads were generated from 500 polio and 2,000 non-polio enterovirus genomes as a training set. By training the algorithm on this dataset we try to determine, on a single read level, which short reads can reliably be labeled as poliovirus and which cannot. After training the deep learning algorithm on the generated reads we were able to calculate the probability with which a read can be assigned to a poliovirus genome or a non-poliovirus genome. We show that the algorithm succeeds in classifying the reads with high accuracy. The probability of assigning the read to the correct class was related to the location in the genome to which the read mapped, which conformed with our expectations since some regions of the genome are more conserved than others. Classifying short reads of high-threat viral pathogens seems to be a promising application of sequence-based deep learning algorithms. Also, recent developments in software and hardware have facilitated the development and training of deep learning algorithms. Further plans of this work are to characterize the hard-to-classify regions of the poliovirus genome, build larger training databases, and expand on the current approach to other viruses.

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., 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.

Highly pathogenic avian influenza (HPAI) H5 clade 2.3.4.4 viruses were first introduced into Europe in late 2014 and re-introduced in late 2016, following detections in Asia and Russia. In contrast to the 2014–15 H5N8 wave, there was substantial local virus amplification in wild birds in Europe in 2016–17 and associated wild bird mortality, with evidence for occasional gene exchange with low pathogenic avian influenza (LPAI) viruses. Since December 2017, several European countries have again reported events or outbreaks with HPAI H5N6 reassortant viruses in both wild birds and poultry, respectively. Previous phylogenetic studies have shown that the two earliest incursions of HPAI H5N8 viruses originated in Southeast Asia and subsequently spread to Europe. In contrast, this study indicates that recent HPAI H5N6 viruses evolved from the H5N8 2016–17 viruses during 2017 by reassortment of a European HPAI H5N8 virus and wild host reservoir LPAI viruses. The genetic and phenotypic differences between these outbreaks and the continuing detections of HPAI viruses in Europe are a cause of concern for both animal and human health. The current co-circulation of potentially zoonotic HPAI and LPAI virus strains in Asia warrants the determination of drivers responsible for the global spread of Asian lineage viruses and the potential threat they pose to public health.

Recently, four new viruses belonging to an unassigned family within the order Picornavirales were identified in excrements of healthy carp (fisavirus) and pigs (posavirus 1, 2 and 3). We report the detection and characterization of a fifth virus present in human faeces. The virus, named human stool-associated RNA virus (husavirus), contains a single ORF encoding a putative 2993 AA polyprotein, with a Hel-Pro-Pol replication block, typical for the Picornavirales. Phylogenetic analysis revealed that the closest relative to husavirus is posavirus 1, and together they cluster with fisavirus, posavirus 2 and 3 and a roundworm (Ascaris suum) derived virus. Husavirus was detected in eight human stool samples collected in 1984 (n = 2), 1985 (n = 4), 1995 (n = 1) and 2014 (n = 1). From three strains of husavirus from 1984 and 1985 the full genome sequence was determined, showing less than 5 % intraspecies variation in the nucleotide composition. The host of this virus remains to be determined.

Lemurs are highly endangered mammals inhabiting the forests of Madagascar. In this study, we performed virus discovery on serum samples collected from 84 wild lemurs and identified viral sequence fragments from 4 novel viruses within the family Flaviviridae, including members of the genera Hepacivirus and Pegivirus. The sifaka hepacivirus (SifHV, two genotypes) and pegivirus (SifPgV, two genotypes) were discovered in the diademed sifaka (Propithecus diadema), while other pegiviral fragments were detected in samples from the indri (Indri indri, IndPgV) and the weasel sportive lemur (Lepilemur mustelinus, LepPgV). Although data are preliminary, each viral species appeared host species-specific and frequent infection was detected (18 of 84 individuals were positive for at least one virus). The complete coding sequence and partial 5’ and 3’ untranslated regions (UTRs) were obtained for SifHV and its genomic organization was consistent with that of other hepaciviruses, with one unique polyprotein and highly structured UTRs. Phylogenetic analyses showed the SifHV belonged to a clade that includes several viral species identified in rodents from Asia and North America, while SifPgV and IndPgV were more closely related to pegiviral species A and C, that include viruses found in humans as well as New- and Old-World monkeys. Our results support the current proposed model of virus-host co-divergence with frequent occurrence of cross-species transmission for these genera and highlight how the discovery of more members of the Flaviviridae can help clarify the ecology and evolutionary history of these viruses. Furthermore, this knowledge is important for conservation and captive management of lemurs.

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.

Viral infections during the prenatal or early childhood periods are one of the environmental factors which might play an etiological role in psychoses. Several studies report higher antibody levels against viruses during pregnancy in blood of mothers of offspring with psychotic disorders, but the presence of such viruses has never been demonstrated. The goal of this study was to investigate the potential association between viral infections during pregnancy and progeny with psychotic disorders and, for this purpose, we performed a nested case control study involving pregnant mothers of offspring with schizophrenia or bipolar disorder with psychotic features (cases, N = 43) and pregnant women with healthy offspring (controls, N = 95). Since several potential viral candidates have been suggested in prior work, a broad-spectrum virus detection system was necessary. A metagenomic analysis performed with the virus discovery method VIDISCA-454 revealed only common blood-associated viruses in all cohorts. However, a significantly lower viral prevalence was detected in the group of cases and in the sub-population of pregnant mothers of offspring with schizophrenia (p < 0.05). Consistent with the existing inverse correlation between the level of these viruses and the immunocompetence of an individual, we hypothesized the presence of a higher immune activity during pregnancy in mothers whose offspring later develop a psychotic disorder as compared to controls. Combining our results with previously available literature data on antibody levels during the gestation period suggests that a more prominent maternal immune activity can be considered a risk factor for developing psychosis. (C) 2015 Elsevier B.V. All rights reserved.

Causative agents for more than 30 percent of respiratory infections remain unidentified, suggesting that unknown respiratory pathogens might be involved. In this study, antibody capture VIDISCA-454 (virus discovery cDNA-AFLP combined with Roche 454 high-throughput sequencing) resulted in the discovery of a novel type of rhinovirus C (RV-C). The virus has an RNA genome of at least 7054 nt and carries the characteristics of rhinovirus C species. The gene encoding viral protein 1, which is used for typing, has only 81% nucleotide sequence identity with the closest known RV-C type, and, therefore, the virus represents the first member of a novel type, named RV-C54.

Porcine stool-associated RNA virus (posavirus), and Human stool-associated RNA virus (husavirus) are viruses in the order Picornavirales recently described in porcine and human fecal samples. The tentative group (Posa and Posa-like viruses: PPLVs) also includes fish stool-associated RNA virus (fisavirus) as well as members detected in insects (Drosophila subobscura and Anopheles sinensis) and parasites (Ascaris suum). As part of an agnostic deep sequencing survey of animal and human viruses in Vietnam, we detected three husaviruses in human fecal samples, two of which share 97–98% amino acid identity to Dutch husavirus strains and one highly divergent husavirus with only 25% amino acid identity to known husaviruses. In addition, the current study found forty-seven complete posavirus genomes from pigs, ten novel rat stool-associated RNA virus genomes (tentatively named rasavirus), and sixteen novel bat stool-associated RNA virus genomes (tentatively named basavirus). The five expected Picornavirales protein domains (helicase, 3C-protease, RNA-dependent RNA polymerase, and two Picornavirus capsid domain) were found to be encoded by all PPLV genomes. In addition, a nucleotide composition analysis revealed that the PPLVs shared compositional properties with arthropod viruses and predicted non-mammalian hosts for all PPLV lineages. The study adds seventy-six genomes to the twenty-nine PPLV genomes currently available and greatly extends our sequence knowledge of this group of viruses within the Picornavirales order.

Several novel clades of astroviruses have recently been identified in human faecal samples. Here, we describe a novel astrovirus-like RNA virus detected in human stools, which we have tentatively named bastrovirus. The genome of this novel virus consists of 6,300 nucleotides organized in three open reading frames. Several sequence divergent strains were detected sharing 67-93 per cent nucleotide identity. Bastrovirus encodes a putative structural protein that is homologous to the capsid protein found in members of the Astroviridae family (45% amino acid identity). The virus also encodes a putative non-structural protein that is genetically distant from astroviruses but shares some homology to the non-structural protein encoded by members of the Hepeviridae family (28% amino acid identity). This novel bastrovirus is present in 8.7 per cent (35/400) of faecal samples collected from 300 HIV-1-positive and 100 HIV-1-negative individuals suggesting common occurrence of the virus. However, whether the source of the virus is infected human cells or other, for example, dietary, remains to be determined.

Coordinated and synchronous virological surveillance for zoonotic viruses in both human clinical cases and animal reservoirs provides an opportunity to identify interspecies virus movement. Rotavirus is an important cause of viral gastroenteritis in humans and animals. We have documented the rotavirus diversity within co-located humans and animals sampled from the Mekong delta region of Vietnam using a primer-independent, agnostic, deep sequencing approach. A total of 296 stool samples (146 from diarrhoeal human patients and 150 from pigs living in the same geographical region) were directly sequenced, generating the genomic sequences of 60 human rotaviruses (all group A) and 31 porcine rotaviruses (13 group A, 7 group B, 6 group C and 5 group H). Phylogenetic analyses showed the co-circulation of multiple distinct rotavirus group A (RVA) genotypes/strains, many of which were divergent from the strain components of licensed RVA vaccines, as well as considerable virus diversity in pigs including full genomes of rotaviruses in groups B, C and H, none of which have been previously reported in Vietnam. Furthermore the detection of an atypical RVA genotype constellation (G4-P[6]-I1-R1-C1-M1-A8-N1-T7-E1-H1) in a human patient and a pig from the same region provides some evidence for a zoonotic event

The list of recently discovered gastrointestinal viruses is expanding rapidly. Whether these agents are actually involved in a disease such as diarrhea is the essential question, yet difficult to answer. In this review a summary of all viruses found in diarrhea is presented, together with the current knowledge about their connection to disease.

Background Gastrointestinal symptoms, in particular diarrhoea, are common in non-treated HIV-1 infected individuals. Although various enteric pathogens have been implicated, the aetiology of diarrhoea remains unexplained in a large proportion of HIV-1 infected patients. Our aim is to identify the cause of diarrhoea for patients that remain negative in routine diagnostics. Methods In this study stool samples of 196 HIV-1 infected persons, including 29 persons with diarrhoea, were examined for enteropathogens and HIV-1. A search for unknown and unexpected viruses was performed using virus discovery cDNA-AFLP combined with Roche-454 sequencing (VIDISCA-454). Results HIV-1 RNA was detected in stool of 19 patients with diarrhoea (66%) compared to 75 patients (45%) without diarrhoea. In 19 of the 29 diarrhoea cases a known enteropathogen could be identified (66%). Next to these known causative agents, a range of recently identified viruses was identified via VIDISCA-454: cosavirus, Aichi virus, human gyrovirus, and non-A non-B hepatitis virus. Moreover, a novel virus was detected which was named immunodeficiency-associated stool virus (IASvirus). However, PCR based screening for these viruses showed that none of these novel viruses was associated with diarrhoea. Notably, among the 34% enteropathogen-negative cases, HIV-1 RNA shedding in stool was more frequently observed (80%) compared to enteropathogen-positive cases (47%), indicating that HIV-1 itself is the most likely candidate to be involved in diarrhoea. Conclusion Unexplained diarrhoea in HIV-1 infected patients is probably not caused by recently described or previously unknown pathogens, but it is more likely that HIV-1 itself plays a role in intestinal mucosal abnormalities which leads to diarrhoea.

Since 2009 pandemic, international health authorities recommended monitoring severe and complicated cases of respiratory disease, that is, severe acute respiratory infection (SARI) and acute respiratory distress syndrome (ARDS). We evaluated the proportion of SARI/ARDS cases and deaths due to influenza A(H1N1)pdm09 infection and the impact of other respiratory viruses during pandemic and postpandemic period (2009-2011) in northern Italy; additionally we searched for unknown viruses in those cases for which diagnosis remained negative. 206 respiratory samples were collected from SARI/ARDS cases and analyzed by real-time RT-PCR/PCR to investigate influenza viruses and other common respiratory pathogens; also, a virus discovery technique (VIDISCA-454) was applied on those samples tested negative to all pathogens. Influenza A(H1N1)pdm09 virus was detected in 58.3% of specimens, with a case fatality rate of 11.3%. The impact of other respiratory viruses was 19.4%, and the most commonly detected viruses were human rhinovirus/enterovirus and influenza A(H3N2). VIDISCA-454 enabled the identification of one previously undiagnosed measles infection. Nearly 22% of SARI/ARDS cases did not obtain a definite diagnosis. In clinical practice, great efforts should be dedicated to improving the diagnosis of severe respiratory disease; the introduction of innovative molecular technologies, as VIDISCA-454, will certainly help in reducing such "diagnostic gap." © 2014 Elena Pariani et al

Discovery of new viruses has been boosted by novel deep sequencing technologies. Currently, many viruses can be identified by sequencing without knowledge of the pathogenicity of the virus. However, attributing the presence of a virus in patient material to a disease in the patient can be a challenge. One approach to meet this challenge is identification of viral sequences based on enrichment by autologous patient antibody capture. This method facilitates identification of viruses that have provoked an immune response within the patient and may increase the sensitivity of the current virus discovery techniques. To demonstrate the utility of this method, virus discovery deep sequencing (VIDISCA-454) was performed on clinical samples from 19 patients: 13 with a known respiratory viral infection and 6 with a known gastrointestinal viral infection. Patient sera was collected from one to several months after the acute infection phase. Input and antibody capture material was sequenced and enrichment was assessed. In 18 of the 19 patients, viral reads from immunogenic viruses were enriched by antibody capture (ranging between 1.5x to 343x in respiratory material, and 1.4x to 53x in stool). Enriched reads were also determined in an identity independent manner by using a novel algorithm Xcompare. In 16 of the 19 patients, 21% to 100% of the enriched reads were derived from infecting viruses. In conclusion, the technique provides a novel approach to specifically identify immunogenic viral sequences among the bulk of sequences which are usually encountered during virus discovery metagenomics.

Virus discovery combining sequence unbiased amplification with next generation sequencing is now state-of-the-art. We have previously determined that the performance of the unbiased amplification technique which is operational at our institute, VIDISCA-454, is efficient when respiratory samples are used as input. The performance of the assay is, however, not known for other clinical materials like blood or stool samples. Here, we investigated the sensitivity of VIDISCA-454 with feces-suspensions and serum samples that are positive and that have been quantified for norovirus and human immunodeficiency virus type 1, respectively. The performance of VIDISCA-454 in serum samples was equal to its performance in respiratory material, with an estimated lower threshold of 1,000 viral genome copies. The estimated threshold in feces-suspension is around 200,000 viral genome copies. The decreased sensitivity in feces suspension is mainly due to sequences that share no recognizable identity with known sequences. Most likely these sequences originate from bacteria and phages which are not completely sequenced.

We have developed a full genome virus detection process that combines sensitive nucleic acid preparation optimised for virus identification in fecal material with Illumina MiSeq sequencing and a novel post-sequencing virus identification algorithm. Enriched viral nucleic acid was converted to double-stranded DNA and subjected to Illumina MiSeq sequencing. The resulting short reads were processed with a novel iterative Python algorithm SLIM for the identification of sequences with homology to known viruses. De novo assembly was then used to generate full viral genomes. The sensitivity of this process was demonstrated with a set of fecal samples from HIV-1 infected patients. A quantitative assessment of the mammalian, plant, and bacterial virus content of this compartment was generated and the deep sequencing data were sufficient to assembly 12 complete viral genomes from 6 virus families. The method detected high levels of enteropathic viruses that are normally controlled in healthy adults, but may be involved in the pathogenesis of HIV-1 infection and will provide a powerful tool for virus detection and for analyzing changes in the fecal virome associated with HIV-1 progression and pathogenesis.

Spread and emergence of West Nile virus (WNV) in Europe have been very different from those observed in North America. Here, we describe key drivers by combining viral genome sequences with epidemiological data and possible factors of spread into phylodynamic models. WNV in Europe has greater lineage diversity than other regions of the world, suggesting repeated introductions and local amplification. Among the six lineages found in Europe, WNV-2a is predominant, has spread to at least 14 countries and evolved into two major co-circulating clusters (A and B). Both of these seem to originate from regions of Central Europe. Viruses of Cluster A emerged earlier and have spread towards the west of Europe with higher genetic diversity. Amongst multiple drivers, high agriculture activities were associated with both spread direction and velocity. Our study suggests future surveillance activities should be strengthened in Central Europe and Southeast European countries, and enhanced monitoring should be targeted to areas with high agriculture activities.