Genome analysis of the novel putative rotavirus species K.

Rotaviruses are causative agents of diarrhea in humans and animals. Currently, the species rotavirus A-J (RVA-RVJ) and the putative species RVK and RVL are defined, mainly based on their genome sequence identities. RVK strains were first identified in 2019 in common shrews (Sorex aranaeus) in Germany; however, only short sequence fragments were available so far. Here, we analyzed the complete coding regions of strain RVK/shrew-wt/GER/KS14-0241/2013, which showed highest sequence identities with RVC. The amino acid sequence identity of VP6, which is used for rotavirus species definition, reached only 51% with other rotavirus reference strains thus confirming classification of RVK as a separate species. Phylogenetic analyses for the deduced amino acid sequences of all 11 virus proteins showed, that for most of them RVK and RVC formed a common branch within the RVA-like phylogenetic clade. Only the tree for the highly variable NSP4 showed a different branching; however, with very low bootstrap support. Comparison of partial nucleotide sequences of other RVK strains from common shrews of different regions in Germany indicated a high degree of sequence variability (61-97% identity) within the putative species. These RVK strains clustered separately from RVC genotype reference strains in phylogenetic trees indicating diversification of RVK independent from RVC. The results indicate that RVK represents a novel rotavirus species, which is most closely related to RVC.

Clonal relation between Salmonella enterica subspecies enterica serovar Dublin strains of bovine and food origin in Germany.

Salmonella enterica subspecies enterica serovar Dublin (S. Dublin) is a host-adapted serovar causing enteritis and/or systemic diseases in cattle. As the serovar is not host-restricted, it may cause infections in other animals, including humans with severe illness and higher mortality rates than other non-typhoidal serovars. As human infections are mainly caused by contaminated milk, milk products and beef, information on the genetic relationship of S. Dublin strains from cattle and food should be evaluated. Whole-genome sequencing (WGS) of 144 S. Dublin strains from cattle and 30 strains from food origin was performed. Multilocus sequence typing (MLST) revealed mostly sequence type ST-10 from both, cattle and food isolates. In total, 14 of 30 strains from food origin were clonally related to at least one strain from cattle, as detected by core-genome single nucleotide polymorphisms typing as well as core-genome MLST. The remaining 16 foodborne strains fit into the genome structure of S. Dublin in Germany without outliers. WGS proved to be a powerful tool not only to gain information on the epidemiology of Salmonella strains but also to detect clonal relations between organisms isolated from different stages of production. This study has shown a high genetic correlation between S. Dublin strains from cattle and food and, therefore, the potential to cause human infections. S. Dublin strains of both origins share an almost identical set of virulence factors, emphasizing their potential to cause severe clinical manifestations in animals, but also in humans and thus the need for effective control of S. Dublin in a farm-to-fork strategy.

Genome sequence analysis of a novel rotavirus strain indicates a broad genetic diversity of rotavirus A in shrews.

Rotavirus A (RVA) is an etiologic agent of diarrhea in humans and animals. It shows a high degree of genetic heterogeneity. Although distinct associations of RVA genotypes with certain host species are common, interspecies-transmission has also been described. Recently, RVA strains, which are genetically distinct and cluster basally to all other RVA strains in phylogenetic trees, have been identified in common shrews (Sorex araneus). Here, the genome sequence analysis of another RVA strain (RVA/Common Shrew-wt/GER/KS11-0893/2010/G42P[58]) from a common shrew from Germany is described. Generally, the strain shows low sequence identities to established strains, which is reflected by the assessment of the novel genotypes G42-P[58]-I32-R28-C24-M24-A39-N28-T28-E32-H28 to its genome segments. Specifically, the strain is phylogenetically distant from previously described RVA strains of common shrews, whereas it is more closely related to other avian and mammalian RVA strains including those from Asian house shrews (Suncus murinus). The results indicate that a broad variety of diverse RVA strains can be found in shrews suggesting a significant role of these animals in rotavirus evolution.

PathoLive-Real-Time Pathogen Identification from Metagenomic Illumina Datasets.

Over the past years, NGS has become a crucial workhorse for open-view pathogen diagnostics. Yet, long turnaround times result from using massively parallel high-throughput technologies as the analysis can only be performed after sequencing has finished. The interpretation of results can further be challenged by contaminations, clinically irrelevant sequences, and the sheer amount and complexity of the data. We implemented PathoLive, a real-time diagnostics pipeline for the detection of pathogens from clinical samples hours before sequencing has finished. Based on real-time alignment with HiLive2, mappings are scored with respect to common contaminations, low-entropy areas, and sequences of widespread, non-pathogenic organisms. The results are visualized using an interactive taxonomic tree that provides an easily interpretable overview of the relevance of hits. For a human plasma sample that was spiked in vitro with six pathogenic viruses, all agents were clearly detected after only 40 of 200 sequencing cycles. For a real-world sample from Sudan, the results correctly indicated the presence of Crimean-Congo hemorrhagic fever virus. In a second real-world dataset from the 2019 SARS-CoV-2 outbreak in Wuhan, we found the presence of a SARS coronavirus as the most relevant hit without the novel virus reference genome being included in the database. For all samples, clinically irrelevant hits were correctly de-emphasized. Our approach is valuable to obtain fast and accurate NGS-based pathogen identifications and correctly prioritize and visualize them based on their clinical significance: PathoLive is open source and available on GitLab and BioConda.

Whole Genome Sequence Analysis of a Prototype Strain of the Novel Putative Rotavirus Species L.

Rotaviruses infect humans and animals and are a main cause of diarrhea. They are non-enveloped viruses with a genome of 11 double-stranded RNA segments. Based on genome analysis and amino acid sequence identities of the capsid protein VP6, the rotavirus species A to J (RVA-RVJ) have been defined so far. In addition, rotaviruses putatively assigned to the novel rotavirus species K (RVK) and L (RVL) have been recently identified in common shrews (Sorex araneus), based on partial genome sequences. Here, the complete genome sequence of strain KS14/0241, a prototype strain of RVL, is presented. The deduced amino acid sequence for VP6 of this strain shows only up to 47% identity to that of RVA to RVJ reference strains. Phylogenetic analyses indicate a clustering separated from the established rotavirus species for all 11 genome segments of RVL, with the closest relationship to RVH and RVJ within the phylogenetic RVB-like clade. The non-coding genome segment termini of RVL showed conserved sequences at the 5'-end (positive-sense RNA strand), which are common to all rotaviruses, and those conserved among the RVB-like clade at the 3'-end. The results are consistent with a classification of the virus into a novel rotavirus species L.

Genetic and biological characteristics of species A rotaviruses detected in common shrews suggest a distinct evolutionary trajectory.

Species A rotaviruses (RVAs) are important aetiological agents of severe diarrhoea in young children. They are also widely distributed in mammals and birds, and increasing evidence indicates the possibility of zoonotic transmission of RVA strains between animals and humans. Moreover, reassortment of the eleven segments of the RVA genome can result in rapid biological changes and may influence pathogenic properties. Here, the nearly complete genome of an RVA strain from a common shrew (Sorex araneus) was sequenced, which showed high nucleotide sequence similarity to additionally determined partial sequences from common shrew RVAs but only very low identity (below 68 per cent) to RVAs from other animal species and humans. New genotypes were assigned to most genome segments of the novel common shrew RVA strain KS14/269, resulting in the genome constellation G39-P[55]-I27-R26-C22-M22-A37-N26-T26-E30-H26. Phylogenetic analyses clustered the common shrew RVAs as ancestral branches of other mammalian and avian RVAs for most of the genome segments, which is in contrast to the phylogeny of the hosts. Nevertheless, conserved sequences typical for all RVAs were identified at the 5'- and 3'- non-coding segment termini. To explore whether the common shrew RVA can exchange genetic material with other mammalian RVAs by reassortment, a reverse genetics system based on the simian RVA strain SA11 was used. However, no viable reassortants could be rescued by exchanging the VP4-, VP6-, or VP7-encoding genome segment alone or in combinations. It can be concluded that highly divergent RVAs are present in common shrews, indicating an evolution of these viruses largely separated from other mammalian and avian RVAs. The zoonotic potential of the virus seems to be low but needs to be further analysed in future.

Future-proofing and maximizing the utility of metadata: The PHA4GE SARS-CoV-2 contextual data specification package.

BACKGROUND: The Public Health Alliance for Genomic Epidemiology (PHA4GE) (https://pha4ge.org) is a global coalition that is actively working to establish consensus standards, document and share best practices, improve the availability of critical bioinformatics tools and resources, and advocate for greater openness, interoperability, accessibility, and reproducibility in public health microbial bioinformatics. In the face of the current pandemic, PHA4GE has identified a need for a fit-for-purpose, open-source SARS-CoV-2 contextual data standard. RESULTS: As such, we have developed a SARS-CoV-2 contextual data specification package based on harmonizable, publicly available community standards. The specification can be implemented via a collection template, as well as an array of protocols and tools to support both the harmonization and submission of sequence data and contextual information to public biorepositories. CONCLUSIONS: Well-structured, rich contextual data add value, promote reuse, and enable aggregation and integration of disparate datasets. Adoption of the proposed standard and practices will better enable interoperability between datasets and systems, improve the consistency and utility of generated data, and ultimately facilitate novel insights and discoveries in SARS-CoV-2 and COVID-19. The package is now supported by the NCBI's BioSample database.

What WGS Reveals about Salmonella enterica subsp. enterica in Wildlife in Germany.

The aim of this study was to gain an overview of the genetic diversity of Salmonella found in wildlife in Germany. We were particularly interested in exploring whether wildlife acts as a reservoir of certain serovars/subtypes or antimicrobial resistance (AMR) genes. Moreover, we wanted to explore the potential of Salmonella in spreading from wildlife to livestock and humans. To answer these questions, we sequenced 260 Salmonella enterica subsp. enterica isolates sampled between 2002 and 2020 from wildlife across Germany, using short-read whole genome sequencing. We found, consistent with previous findings, that some Salmonella sequence types are associated with certain animal species, such as S. Choleraesuis ST145 with wild boar and S. Enteritidis ST183 with hedgehogs. Antibiotic resistance was detected in 14.2% of all isolates, with resistance against important WATCH group antibiotics present in a small number of isolates. We further found that wildlife isolates do not form separate phylogenetic clusters distant to isolates from domestic animals and foodstuff, thus indicating frequent transmission events between these reservoirs. Overall, our study shows that Salmonella in German wildlife are diverse, with a low AMR burden and close links to Salmonella populations of farm and food-production environments.

Decentralized Investigation of Bacterial Outbreaks Based on Hashed cgMLST.

Whole-genome sequencing (WGS)-based outbreak investigation has proven to be a valuable method for the surveillance of bacterial pathogens. Its utility has been successfully demonstrated using both gene-by-gene (cgMLST or wgMLST) and single-nucleotide polymorphism (SNP)-based approaches. Among the obstacles of implementing a WGS-based routine surveillance is the need for an exchange of large volumes of sequencing data, as well as a widespread reluctance to share sequence and metadata in public repositories, together with a lacking standardization of suitable bioinformatic tools and workflows. To address these issues, we present chewieSnake, an intuitive and simple-to-use cgMLST workflow. ChewieSnake builds on the allele calling software chewBBACA and extends it by the concept of allele hashing. The resulting hashed allele profiles can be readily compared between laboratories without the need of a central allele nomenclature. The workflow fully automates the computation of the allele distance matrix, cluster membership, and phylogeny and summarizes all important findings in an interactive HTML report. Furthermore, chewieSnake can join allele profiles generated at different laboratories and identify shared clusters, including a stable and intercommunicable cluster nomenclature, thus facilitating a joint outbreak investigation. We demonstrate the feasibility of the proposed approach with a thorough method comparison using publically available sequencing data for Salmonella enterica. However, chewieSnake is readily applicable to all bacterial taxa, provided that a suitable cgMLST scheme is available. The workflow is freely available as an open-source tool and can be easily installed via conda or docker.

Species-Specific Quality Control, Assembly and Contamination Detection in Microbial Isolate Sequences with AQUAMIS.

Sequencing of whole microbial genomes has become a standard procedure for cluster detection, source tracking, outbreak investigation and surveillance of many microorganisms. An increasing number of laboratories are currently in a transition phase from classical methods towards next generation sequencing, generating unprecedented amounts of data. Since the precision of downstream analyses depends significantly on the quality of raw data generated on the sequencing instrument, a comprehensive, meaningful primary quality control is indispensable. Here, we present AQUAMIS, a Snakemake workflow for an extensive quality control and assembly of raw Illumina sequencing data, allowing laboratories to automatize the initial analysis of their microbial whole-genome sequencing data. AQUAMIS performs all steps of primary sequence analysis, consisting of read trimming, read quality control (QC), taxonomic classification, de-novo assembly, reference identification, assembly QC and contamination detection, both on the read and assembly level. The results are visualized in an interactive HTML report including species-specific QC thresholds, allowing non-bioinformaticians to assess the quality of sequencing experiments at a glance. All results are also available as a standard-compliant JSON file, facilitating easy downstream analyses and data exchange. We have applied AQUAMIS to analyze ~13,000 microbial isolates as well as ~1000 in-silico contaminated datasets, proving the workflow's ability to perform in high throughput routine sequencing environments and reliably predict contaminations. We found that intergenus and intragenus contaminations can be detected most accurately using a combination of different QC metrics available within AQUAMIS.

Outcome of Different Sequencing and Assembly Approaches on the Detection of Plasmids and Localization of Antimicrobial Resistance Genes in Commensal Escherichia coli.

Antimicrobial resistance (AMR) is a major threat to public health worldwide. Currently, AMR typing changes from phenotypic testing to whole-genome sequence (WGS)-based detection of resistance determinants for a better understanding of the isolate diversity and elements involved in gene transmission (e.g., plasmids, bacteriophages, transposons). However, the use of WGS data in monitoring purposes requires suitable techniques, standardized parameters and approved guidelines for reliable AMR gene detection and prediction of their association with mobile genetic elements (plasmids). In this study, different sequencing and assembly strategies were tested for their suitability in AMR monitoring in Escherichia coli in the routines of the German National Reference Laboratory for Antimicrobial Resistances. To assess the outcomes of the different approaches, results from in silico predictions were compared with conventional phenotypic- and genotypic-typing data. With the focus on (fluoro)quinolone-resistant E.coli, five qnrS-positive isolates with multiple extrachromosomal elements were subjected to WGS with NextSeq (Illumina), PacBio (Pacific BioSciences) and ONT (Oxford Nanopore) for in depth characterization of the qnrS1-carrying plasmids. Raw reads from short- and long-read sequencing were assembled individually by Unicycler or Flye or a combination of both (hybrid assembly). The generated contigs were subjected to bioinformatics analysis. Based on the generated data, assembly of long-read sequences are error prone and can yield in a loss of small plasmid genomes. In contrast, short-read sequencing was shown to be insufficient for the prediction of a linkage of AMR genes (e.g., qnrS1) to specific plasmid sequences. Furthermore, short-read sequencing failed to detect certain duplications and was unsuitable for genome finishing. Overall, the hybrid assembly led to the most comprehensive typing results, especially in predicting associations of AMR genes and mobile genetic elements. Thus, the use of different sequencing technologies and hybrid assemblies currently represents the best approach for reliable AMR typing and risk assessment.

Toward an Integrated Genome-Based Surveillance of Salmonella enterica in Germany.

Despite extensive monitoring programs and preventative measures, Salmonella spp. continue to cause tens of thousands human infections per year, as well as many regional and international food-borne outbreaks, that are of great importance for public health and cause significant socio-economic costs. In Germany, salmonellosis is the second most common cause of bacterial diarrhea in humans and is associated with high hospitalization rates. Whole-genome sequencing (WGS) combined with data analysis is a high throughput technology with an unprecedented discriminatory power, which is particularly well suited for targeted pathogen monitoring, rapid cluster detection and assignment of possible infection sources. However, an effective implementation of WGS methods for large-scale microbial pathogen detection and surveillance has been hampered by the lack of standardized methods, uniform quality criteria and strategies for data sharing, all of which are essential for a successful interpretation of sequencing data from different sources. To overcome these challenges, the national GenoSalmSurv project aims to establish a working model for an integrated genome-based surveillance system of Salmonella spp. in Germany, based on a decentralized data analysis. Backbone of the model is the harmonization of laboratory procedures and sequencing protocols, the implementation of open-source bioinformatics tools for data analysis at each institution and the establishment of routine practices for cross-sectoral data sharing for a uniform result interpretation. With this model, we present a working solution for cross-sector interpretation of sequencing data from different sources (such as human, veterinarian, food, feed and environmental) and outline how a decentralized data analysis can contribute to a uniform cluster detection and facilitate outbreak investigations.

Proficiency Testing of Metagenomics-Based Detection of Food-Borne Pathogens Using a Complex Artificial Sequencing Dataset.

Metagenomics-based high-throughput sequencing (HTS) enables comprehensive detection of all species comprised in a sample with a single assay and is becoming a standard method for outbreak investigation. However, unlike real-time PCR or serological assays, HTS datasets generated for pathogen detection do not easily provide yes/no answers. Rather, results of the taxonomic read assignment need to be assessed by trained personnel to gain information thereof. Proficiency tests are important instruments of validation, harmonization, and standardization. Within the European Union funded project COMPARE [COllaborative Management Platform for detection and Analyses of (Re-) emerging and foodborne outbreaks in Europe], we conducted a proficiency test to scrutinize the ability to assess diagnostic metagenomics data. An artificial dataset resembling shotgun sequencing of RNA from a sample of contaminated trout was provided to 12 participants with the request to provide a table with per-read taxonomic assignments at species level and a report with a summary and assessment of their findings, considering different categories like pathogen, background, or contaminations. Analysis of the read assignment tables showed that the software used reliably classified the reads taxonomically overall. However, usage of incomplete reference databases or inappropriate data pre-processing caused difficulties. From the combination of the participants' reports with their read assignments, we conclude that, although most species were detected, a number of important taxa were not or not correctly categorized. This implies that knowledge of and awareness for potentially dangerous species and contaminations need to be improved, hence, capacity building for the interpretation of diagnostic metagenomics datasets is necessary.

German-Wide Interlaboratory Study Compares Consistency, Accuracy and Reproducibility of Whole-Genome Short Read Sequencing.

We compared the consistency, accuracy and reproducibility of next-generation short read sequencing between ten laboratories involved in food safety (research institutes, state laboratories, universities and companies) from Germany and Austria. Participants were asked to sequence six DNA samples of three bacterial species (Campylobacter jejuni, Listeria monocytogenes and Salmonella enterica) in duplicate, according to their routine in-house sequencing protocol. Four different types of Illumina sequencing platforms (MiSeq, NextSeq, iSeq, NovaSeq) and one Ion Torrent sequencing instrument (S5) were involved in the study. Sequence quality parameters were determined for all data sets and centrally compared between laboratories. SNP and cgMLST calling were performed to assess the reproducibility of sequence data collected for individual samples. Overall, we found Illumina short read data to be more accurate (higher base calling accuracy, fewer miss-assemblies) and consistent (little variability between independent sequencing runs within a laboratory) than Ion Torrent sequence data, with little variation between the different Illumina instruments. Two laboratories with Illumina instruments submitted sequence data with lower quality, probably due to the use of a library preparation kit, which shows difficulty in sequencing low GC genome regions. Differences in data quality were more evident after assembling short reads into genome assemblies, with Ion Torrent assemblies featuring a great number of allele differences to Illumina assemblies. Clonality of samples was confirmed through SNP calling, which proved to be a more suitable method for an integrated data analysis of Illumina and Ion Torrent data sets in this study.

First Detection of GES-5-Producing Escherichia coli from Livestock-An Increasing Diversity of Carbapenemases Recognized from German Pig Production.

Resistance to carbapenems due to carbapenemase-producing Enterobacteriaceae (CPE) is an increasing threat to human health worldwide. In recent years, CPE could be found only sporadically from livestock, but concern rose that livestock might become a reservoir for CPE. In 2019, the first GES carbapenemase-producing Escherichia coli from livestock was detected within the German national monitoring on antimicrobial resistance. The isolate was obtained from pig feces and was phenotypically resistant to meropenem and ertapenem. The isolate harbored three successive blaGES genes encoding for GES-1, GES-5 and GES-5B in an incomplete class-I integron on a 12 kb plasmid (pEC19-AB02908; Acc. No. MT955355). The strain further encoded for virulence-associated genes typical for uropathogenic E. coli, which might hint at an increased pathogenic potential. The isolate produced the third carbapenemase detected from German livestock. The finding underlines the importance CPE monitoring and detailed characterization of new isolates.

Complete Genome Sequence of Salmonella enterica subsp. diarizonae Serovar 61:k:1,5,(7) Strain 14-SA00836-0, Isolated from Human Urine.

Salmonella enterica subsp. diarizonae serovar 61:k:1,5,(7) is commonly associated with sheep. Occasionally, the serovar has been found to also infect humans. Here, we report the complete genome sequence of strain 14-SA00836-0, isolated from human urine. To our knowledge, this is the first reported complete genome sequence of this serovar isolated from a human clinical sample.

Whole genome sequence analysis of cell culture-adapted rotavirus A strains from chicken.

Rotavirus A (RVA) is a major cause of gastroenteritis in humans and mammalian animals, and has also been abundantly detected in avian species. Avian RVA infection is associated with diarrhea, reduced growth and increased mortality, leading to economic losses in the poultry industry. Avian RVA forms a unique genetic clade within the whole RVA species. However, up to now, only a few full-length avian RVA genomes have been published and only a small number of avian RVA strains have been adapted to grow in cell culture for subsequent studies. Here, the four cell culture-adapted chicken RVA strains 02V0002G3, 04V0027G6, 05V0500F6 and 06V0661G1 were characterized in more detail. Transmission electron microscopy of the viruses derived from culture supernatant showed a typical triple-layered morphology of rotavirus particles; in addition, strain 06V0661G1 showed a high proportion of double-layered particles. The (nearly) complete genome sequences of the viruses were determined using next-generation sequencing (NGS). The resulting sequences were compared to full-length or partial sequences of the strains previously determined using Sanger sequencing; and a few nucleotide mismatches, some of them resulting in amino acid substitutions, were identified. The genomes of strains 02V0002G3, 04V0027G6 and 05V0500F6 were closely related to each other showing a G19-P[30]-I11-R6-C6-M7-A16-N6-T8-E10-H8 genotype constellation. Strain 06V0661G1 carries the VP4 genotype P[31] in the same genetic backbone like the other strains. However, further sequence analysis showed that the genes of this strain, especially that encoding NSP3, clustered more separately from the other strains in phylogenetic trees. The characterized cell culture-adapted chicken RVA strains may be useful for future studies investigating genetic diversity and replication of avian rotaviruses, as well as for the development of vaccines and diagnostic tools.

Establishment of a Plasmid-Based Reverse Genetics System for the Cell Culture-Adapted Hepatitis E Virus Genotype 3c Strain 47832c.

The hepatitis E virus (HEV) causes acute and chronic hepatitis in humans. Investigation of HEV replication is hampered by the lack of broadly applicable, efficient cell culture systems and tools for site-directed mutagenesis of HEV. The cell culture-adapted genotype 3c strain 47832c, which represents a typical genotype predominantly detected in Europe, has previously been used for several basic and applied research studies. Here, a plasmid-based reverse genetics system was developed for this strain, which efficiently rescued the infectious virus without the need for in vitro RNA transcription. The cotransfection of T7 RNA polymerase-expressing BSR/T7 cells with one plasmid encoding the full-length viral genome and two helper plasmids encoding vaccinia virus capping enzymes resulted in the production of infectious HEV, which could be serially passaged on A549/D3 cells. The parental and recombinant virus exhibited similar replication kinetics. A single point mutation creating an additional restriction enzyme site could be successfully introduced into the virus genome of progeny virus, indicating that the system is suitable for site-directed mutagenesis. This system is the first plasmid-based HEV reverse genetics system, as well as the first reverse genetics system for HEV genotype 3c, and should therefore be of broad use for basic and applied HEV research.