Search

ABSTRACT Metagenomic shotgun sequencing (mNGS) can serve as a generic molecular diagnostic tool. An mNGS proficiency test (PT) was performed in six European veterinary and public health laboratories to detect porcine astroviruses in fecal material and the extracted RNA. While different mNGS workflows for the generation of mNGS data were used in the different laboratories, the bioinformatic analysis was standardized using a metagenomic read classifier as well as read mapping to selected astroviral reference genomes to assess the semiquantitative representation of astrovirus species mixtures. All participants successfully identified and classified most of the viral reads to the two dominant species. The normalized read counts obtained by aligning reads to astrovirus reference genomes by Bowtie2 were in line with Kraken read classification counts. Moreover, participants performed well in terms of repeatability when the fecal sample was tested in duplicate. However, the normalized read counts per detected astrovirus species differed substantially between participants, which was related to the different laboratory methods used for data generation. Further modeling of the mNGS data indicated the importance of selecting appropriate reference data for mNGS read classification. As virus- or sample-specific biases may apply, caution is needed when extrapolating this swine feces-based PT for the detection of other RNA viruses or using different sample types. The suitability of experimental design to a given pathogen/sample matrix combination, quality assurance, interpretation, and follow-up investigation remain critical factors for the diagnostic interpretation of mNGS results.IMPORTANCEMetagenomic shotgun sequencing (mNGS) is a generic molecular diagnostic method, involving laboratory preparation of samples, sequencing, bioinformatic analysis of millions of short sequences, and interpretation of the results. In this paper, we investigated the performance of mNGS on the detection of porcine astroviruses, a model for RNA viruses in a pig fecal material, among six European veterinary and public health laboratories. We showed that different methods for data generation affect mNGS performance among participants and that the selection of reference genomes is crucial for read classification. Follow-up investigation remains a critical factor for the diagnostic interpretation of mNGS results. The paper contributes to potential improvements of mNGS as a diagnostic tool in clinical settings.

Usutu virus (USUV) is an emerging flavivirus that is maintained in an enzootic cycle with mosquitoes as vectors and birds as amplifying hosts. In Europe, the virus has caused mass mortality of wild birds, mainly among Common Blackbird (Turdus merula) populations. While mosquitoes are the primary vectors for USUV, Common Blackbirds and other avian species are exposed to other arthropod ectoparasites, such as ticks. It is unknown, however, if ticks can maintain and transmit USUV. We addressed this question using in vitro and in vivo experiments and field collected data. USUV replicated in IRE/CTVM19 Ixodes ricinus tick cells and in injected ticks. Moreover, I. ricinus nymphs acquired the virus via artificial membrane blood-feeding and maintained the virus for at least 70 days. Transstadial transmission of USUV from nymphs to adults was confirmed in 4.9% of the ticks. USUV disseminated from the midgut to the haemocoel, and was transmitted via the saliva of the tick during artificial membrane blood-feeding. We further explored the role of ticks by monitoring USUV in questing ticks and in ticks feeding on wild birds in the Netherlands between 2016 and 2019. In total, 622 wild birds and the Ixodes ticks they carried were tested for USUV RNA. Of these birds, 48 (7.7%) carried USUV-positive ticks. The presence of negative-sense USUV RNA in ticks, as confirmed via small RNA-sequencing, showed active virus replication. In contrast, we did not detect USUV in 15,381 questing ticks collected in 2017 and 2019. We conclude that I. ricinus can be infected with USUV and can transstadially and horizontally transmit USUV. However, in comparison to mosquito-borne transmission, the role of I. ricinus ticks in the epidemiology of USUV is expected to be minor.

Abstract Systematic monitoring of SARS-CoV-2 co-infections between different lineages and assessing the risk of intra-host recombinant emergence are crucial for forecasting viral evolution. Here we present a comprehensive analysis of more than 2 million SARS-CoV-2 raw read datasets submitted to the European COVID-19 Data Portal to identify co-infections and intra-host recombination. Co-infection was observed in 0.35% of the investigated cases. Two independent procedures were implemented to detect intra-host recombination. We show that sensitivity is predominantly determined by the density of lineage-defining mutations along the genome, thus we used an expanded list of mutually exclusive defining mutations of specific variant combinations to increase statistical power. We call attention to multiple challenges rendering recombinant detection difficult and provide guidelines for the reduction of false positives arising from chimeric sequences produced during PCR amplification. Additionally, we identify three recombination hotspots of Delta – Omicron BA.1 intra-host recombinants.

Abstract Background Healthcare facilities have been challenged by the risk of SARS-CoV-2 transmission between healthcare workers (HCW) and patients. During the first wave of the COVID-19 pandemic, infections among HCW were observed, questioning infection prevention and control (IPC) measures implemented at that time. Aim This study aimed to identify nosocomial transmission routes of SARS-CoV-2 between HCW and patients in a tertiary care hospital. Methods All SARS-CoV-2 PCR positive HCW and patients identified between 1 March and 19 May 2020, were included in the analysis. Epidemiological data were collected from patient files and HCW contact tracing interviews. Whole genome sequences of SARS-CoV-2 were generated using Nanopore sequencing (WGS). Epidemiological clusters were identified, whereafter WGS and epidemiological data were combined for re-evaluation of epidemiological clusters and identification of potential transmission clusters. HCW infections were further classified into categories based on the likelihood that the infection was acquired via nosocomial transmission. Secondary cases were defined as COVID-19 cases in our hospital, part of a transmission cluster, of which the index case was either a patient or HCW from our hospital. Findings The study population consisted of 293 HCW and 245 patients. Epidemiological data revealed 36 potential epidemiological clusters, with an estimated 222 (75.7%) HCW as secondary cases. WGS results were available for 195 HCW (88.2%) and 20 patients (12.8%) who belonged to an epidemiological cluster. Re-evaluation of the epidemiological clusters, with the available WGS data identified 31 transmission clusters with 65 (29.4%) HCW as secondary cases. Transmission clusters were all part of 18 (50.0%) previously determined epidemiological clusters, demonstrating that several larger outbreaks actually consisted, of several smaller transmission clusters. A total of 21 (7.2%) HCW infections were classified as from confirmed nosocomial, of which 18 were acquired from another HCW and 3 from a patient. Conclusion The majority of SARS-CoV-2 infections among HCW could be attributed to community-acquired infection. Infections among HCW that could be classified as due to nosocomial transmission, were mainly caused by HCW-to-HCW transmission rather than patient-to-HCW transmission. It is important to recognize the uncertainties of cluster analyses based solely on epidemiological data.

In August 2021, a large-scale mortality event affected harbor porpoises (Phocoena phocoena) in the Netherlands. Pathology and ancillary testing of 22 animals indicated that the most likely cause of death was Erysipelothrix rhusiopathiae infection. This zoonotic agent poses a health hazard for cetaceans and possibly for persons handling cetacean carcasses.

Inleiding: Een grote SARS-CoV-2 uitbraak, resulterend in een hoge morbiditeit en sterfte onder bewoners, vond plaats in een Rotterdams verpleeghuis in de periode voordat vaccins beschikbaar waren. Wij beschrijven de uitbraak en reflecteren over geleerde lessen voor toekomstige uitbraakbestrijding nu SARS-CoV-2 evolueert. Methoden: Diagnosedatum en demografische gegevens (leeftijd, geslacht, etage) van positieve bewoners werden verzameld, evenals bijgewoonde groepsbijeenkomsten. Neus-keelswabs van bewoners/personeel werden getest op SARS-CoV-2 door middel van een polymerase chain reactie (PCR), waarbij een deel werd getypeerd middels Whole Genome Sequencing (WGS). Ook werd het ventilatiesysteem onderzocht. Resultaten: Tussen 28/11/2020 en 23/12/2020 testten 106/114 verpleeghuisbewoners (93%) SARS-CoV-2 positief (64% vrouwen, mediane leeftijd 89,5 jaar), evenals 76/164 (46%) personeelsleden. WGS toonde een groot cluster aan, maar er werd geen ‘super preader event’ geïdentificeerd, noch problemen met de ventilatie. De sterfte onder besmette bewoners was 42%. Discussie: Deze uitbraak ontwikkelde zich heel snel, mogelijk doordat personeel in de hectiek niet altijd consequent de persoonlijke beschermingsmiddelen gebruikte. De hoge sterfte is deels te verklaren uit de hoge mediane leeftijd. Aangezien SARS-CoV-2-uitbraken in de toekomst opnieuw zullen optreden, dienen risico’s in verpleeghuizen beperkt te worden, o.a. door het handhaven van een hoge vaccinatiegraad, met oog voor kwaliteit van leven van bewoners en personeel.

Abstract Background In fall 2020 when schools in the Netherlands operated under a limited set of COVID-19 measures, we conducted outbreaks studies in four secondary schools to gain insight in the level of school transmission and the role of SARS-CoV-2 transmission via air and surfaces. Methods Outbreak studies were performed between 11 November and 15 December 2020 when the wild-type variant of SARS-CoV-2 was dominant. Clusters of SARS-CoV-2 infections within schools were identified through a prospective school surveillance study. All school contacts of cluster cases, irrespective of symptoms, were invited for PCR testing twice within 48 h and 4–7 days later. Combined NTS and saliva samples were collected at each time point along with data on recent exposure and symptoms. Surface and active air samples were collected in the school environment. All samples were PCR-tested and sequenced when possible. Results Out of 263 sampled school contacts, 24 tested SARS-CoV-2 positive (secondary attack rate 9.1%), of which 62% remained asymptomatic and 42% had a weakly positive test result. Phylogenetic analysis on 12 subjects from 2 schools indicated a cluster of 8 and 2 secondary cases, respectively, but also other distinct strains within outbreaks. Of 51 collected air and 53 surface samples, none were SARS-CoV-2 positive. Conclusion Our study confirmed within school SARS-CoV-2 transmission and substantial silent circulation, but also multiple introductions in some cases. Absence of air or surface contamination suggests environmental contamination is not widespread during school outbreaks.

We detected a highly divergent SARS-CoV-2 Alpha variant in an immunocompromised person several months after the latest detection of the Alpha variant in the Netherlands. The patient was infected for 42 weeks despite several treatment regimens and disappearance of most clinical symptoms. We identified several potential immune escape mutations in the spike protein.

Abstract Background To understand the dynamics of infectious diseases, genomic epidemiology is increasingly advocated, with a need for rapid generation of genetic sequences during outbreaks for public health decision making. Here, we explore the use of metagenomic sequencing compared to specific amplicon- and capture-based sequencing, both on the Nanopore and the Illumina platform for generation of whole genomes of Usutu virus, Zika virus, West Nile virus, and Yellow Fever virus. Results We show that amplicon-based Nanopore sequencing can be used to rapidly obtain whole genome sequences in samples with a viral load up to Ct 33 and capture-based Illumina is the most sensitive method for initial virus determination. Conclusions The choice of sequencing approach and platform is important for laboratories wishing to start whole genome sequencing. Depending on the purpose of genome sequencing the best choice can differ. The insights presented in this work and the shown differences in data characteristics can guide labs to make a well informed choice.

SARS-CoV-2 was detected in mink farms in the Netherlands in the first wave of the pandemic with evidence of human-to-mink and mink-to-human transmission. Here, the authors investigate this outbreak using phylodynamic analysis and show that personnel links and spatial proximity are predictors of transmission between farms.

Wild birds are reservoirs of several zoonotic arboviruses including West Nile virus (WNV) and Usutu virus (USUV), and are often monitored as indicators for virus introduction and spread. To optimize the bird surveillance for arboviruses in the Netherlands and to explore the possibilities for citizen science in surveillance, we investigated the suitability of using alternative sample types from live and dead birds. The sensitivity of molecular detection via RT-PCR of viral RNA in feather, heart, lung, throat and cloaca swabs from dead birds, and serum, dried blood spots (DBS) and throat and cloaca swabs from live birds were compared. IgY antibody detection was also assessed from DBS relative to serum on protein-microarray and virus neutralization test. Feathers showed a high detection sensitivity for USUV RNA in both live and dead birds, and no significant decrease was observed in the RNA loads in the feathers after being stored dry at room temperature for 43 days. Additionally, viral RNAs extracted from feathers of day 0 and 43 were successfully sequenced. The results indicated no statistical significant difference in sensitivity and viral loads detection in heart, spleen, and lung relative to corresponding brain samples in dead birds. In live birds, viral RNA loads did not differ between throat and cloaca swabs. This study identified less-invasive sample types that allows involvement of citizens in collecting samples from wild birds for arbovirus surveillance. Sensitivity and specificity of DBS-based antibody detections were significantly lower and therefore need optimization.

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.

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.

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

Background: SARS-CoV-2 has spread globally. Currently, literature of SARS-CoV-2 in neonates is scarce. We present a case of a neonate with a high viral load and prolonged virus shedding. Methods: Epidemiology, clinical characteristics, treatment, laboratory data and follow-up information and the treatment of a neonate with COVID-19 were recorded. Results: A 7-day-old boy was admitted to the hospital with fever, lethargy and apnoea. He was found SARS-CoV-2 RNA positive with an exceptionally high viral load in nasopharyngeal swab and stool. The father and two maternity nurses at home had detectable SARS-CoV-2 RNA as well. Sequencing showed all strains belonged to the same cluster. The father was asymptomatic and the maternity nurses developed symptoms after visiting. In the mother, no SARS-CoV-2 RNA could be found. Six days after admission, the neonate was discharged after clinical improvement with oral antibiotics because of a possible pyelonephritis. Monitoring the course of this infection showed that SARS-CoV-2 RNA was detectable in the nasopharynx until day 19 and in stool until day 42 after symptom onset. Conclusions: This case shows that neonates can have a high viral load of SARS-CoV-2 and can shed the virus for over one month in stool. Despite the high viral load in the neonate, the mother and a sibling did not get infected.

Identifying the origin of the rabies virus (RABV) infection may have significant implications for control measures. Here, we identified the source of a RABV infection of two Nepalese migrants in Qatar by comparing their RABV genomes with RABV genomes isolated from the brains of a RABV infected camel and fox from Qatar.

Abstract Foodborne viral infections rank among the top 5 causes of disease, with noroviruses and hepatitis A causing the greatest burden globally. Contamination of foods by infected food handlers or through environmental pollution are the main sources of foodborne illness, with a lesser role for consumption of products from infected animals. Viral partial genomic sequencing has been used for more than two decades to track foodborne outbreaks and whole genome or metagenomics next-generation-sequencing (NGS) are new additions to the toolbox of food microbiology laboratories. We discuss developments in the field of targeted and metagenomic NGS, with an emphasis on application in food virology, the challenges and possible solutions towards future routine application.

An outbreak of SARS-CoV-2 Alpha variant (Pango lineage B.1.1.7) was detected at a primary school (School X) in Lansingerland, the Netherlands, in December 2020. The outbreak was studied retrospectively, and population-based screening was used to assess the extent of virus circulation and decelerate transmission. Cases were SARS-CoV-2 laboratory confirmed and were residents of Lansingerland (November 16th 2020 until February 22th 2021), or had an epidemiological link with School X or neighbouring schools. The SARS-CoV-2 variant was determined using variant PCR or whole genome sequencing. A questionnaire primarily assessed clinical symptoms. A total of 77 Alpha variant cases were found with an epidemiological link to School X, 16 Alpha variant cases linked to the neighbouring schools, and 146 Alpha variant cases among residents of Lansingerland without a link to the schools. The mean number of self-reported symptoms was not significantly different among Alpha variant infected individuals compared to non-Alpha infected individuals. The secondary attack rate (SAR) among Alpha variant exposed individuals in households was 52% higher compared to non-Alpha variant exposed individuals (p = 0.010), with the mean household age, and mean number of children and adults per household as confounders. Sequence analysis of 60 Alpha variant sequences obtained from cases confirmed virus transmission between School X and neighbouring schools, and showed that multiple introductions of the Alpha variant had already taken place in Lansingerland at the time of the study. The alpha variant caused a large outbreak at both locations of School X, and subsequently spread to neighbouring schools, and households. Population-based screening (together with other public health measures) nearly stopped transmission of the outbreak strain, but did not prevent variant replacement in the Lansingerland municipality.

Rabies virus (RABV) has a broad host range and infects multiple cell types throughout the infection cycle. Next-generation sequencing (NGS) and minor variant analysis are powerful tools for studying virus populations within specific hosts and tissues, leading to novel insights into the mechanisms of host-switching and key factors for infecting specific cell types. In this study we investigated RABV populations and minor variants in both original (non-passaged) samples and in vitro-passaged isolates of various CNS regions (hippocampus, medulla oblongata and spinal cord) of a fatal human rabies case, and of multiple CNS and non-CNS tissues of experimentally infected mice. No differences in virus populations were detected between the human CNS regions, and only one non-synonymous single nucleotide polymorphism (SNP) was detected in the fifth in vitro passage of virus isolated from the spinal cord. However, the appearance of this SNP shows the importance of sequencing newly passaged virus stocks before further use. Similarly, we did not detect apparent differences in virus populations isolated from different CNS and non-CNS tissues of experimentally infected mice. Sequencing of viruses obtained from pharyngeal swab and salivary gland proved difficult, and we propose methods for improving sampling.

SARS-CoV-2 outbreaks on 69 Dutch mink farms in 2020 were studied to identify risk factors for virus introduction and transmission and to improve surveillance and containment measures. Clinical signs, laboratory test results, and epidemiological aspects were investigated, such as the date and reason of suspicion, housing, farm size and distances, human contact structure, biosecurity measures, and presence of wildlife, pets, pests, and manure management. On seven farms, extensive random sampling was performed, and age, coat color, sex, and clinical signs were recorded. Mild to severe respiratory signs and general diseases such as apathy, reduced feed intake, and increased mortality were detected on 62/69 farms. Throat swabs were more likely to result in virus detection than rectal swabs. Clinical signs differed between virus clusters and were more severe for dark-colored mink, males, and animals infected later during the year. Geographical clustering was found for one virus cluster. Shared personnel could explain some cases, but other transmission routes explaining farm-to-farm spread were not elucidated. An early warning surveillance system, strict biosecurity measures, and a (temporary) ban on mink farming and vaccinating animals and humans can contribute to reducing the risks of the virus spreading and acquisition of potential mutations relevant to human and animal health.

We tested samples collected from camels, camel workers, and other animals in Sudan and Qatar in 2015 and 2017 for evidence of Middle East respiratory syndrome coronavirus (MERS-CoV) infection. MERS-CoV antibodies were abundant in Sudan camels, but we found no evidence of MERS-CoV infection in camel workers, other livestock, or bats.

Market surveillance showed continuing circulation of avian influenza A(H5N6) virus in live poultry markets in Guangdong Province in 2017, despite compulsory vaccination for avian influenza A(H5Nx) and A(H7N9). We analyzed H5N6 viruses from 2014–2018 from Guangdong Province, revealing antigenic drift and decreased antibody response against the vaccine strain in vaccinated chickens.

The Usutu virus (USUV) is a mosquito-borne zoonotic flavivirus. Despite its continuous circulation in Europe, knowledge on the pathology, cellular and tissue tropism and pathogenetic potential of different circulating viral lineages is still fragmentary. Here, macroscopic and microscopic evaluations are performed in association with the study of cell and tissue tropism and comparison of lesion severity of two circulating virus lineages (Europe 3; Africa 3) in 160 Eurasian blackbirds (Turdus merula) in the Netherlands. Results confirm hepatosplenomegaly, coagulative necrosis and lymphoplasmacytic inflammation as major patterns of lesions and, for the first time, vasculitis as a novel virus-associated lesion. A USUV and Plasmodium spp. co-infection was commonly identified. The virus was associated with lesions by immunohistochemistry and was reported most commonly in endothelial cells and blood circulating and tissue mononucleated cells, suggesting them as a major route of entry and spread. A tropism for mononuclear phagocytes cells was further supported by viral labeling in multinucleated giant cells. The involvement of ganglionic neurons and epithelial cells of the gastrointestinal tract suggests a possible role of oral transmission, while the involvement of feather follicle shafts and bulbs suggests their use as a diagnostic sample for live bird testing. Finally, results suggest similar pathogenicity for the two circulating lineages.

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.

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.

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.

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.