The virome of German bats: comparing virus discovery approaches.

Bats are known to be reservoirs of several highly pathogenic viruses. Hence, the interest in bat virus discovery has been increasing rapidly over the last decade. So far, most studies have focused on a single type of virus detection method, either PCR, virus isolation or virome sequencing. Here we present a comprehensive approach in virus discovery, using all three discovery methods on samples from the same bats. By family-specific PCR screening we found sequences of paramyxoviruses, adenoviruses, herpesviruses and one coronavirus. By cell culture we isolated a novel bat adenovirus and bat orthoreovirus. Virome sequencing revealed viral sequences of ten different virus families and orders: three bat nairoviruses, three phenuiviruses, one orbivirus, one rotavirus, one orthoreovirus, one mononegavirus, five parvoviruses, seven picornaviruses, three retroviruses, one totivirus and two thymoviruses were discovered. Of all viruses identified by family-specific PCR in the original samples, none was found by metagenomic sequencing. Vice versa, none of the viruses found by the metagenomic virome approach was detected by family-specific PCRs targeting the same family. The discrepancy of detected viruses by different detection approaches suggests that a combined approach using different detection methods is necessary for virus discovery studies.

First detection of bat-borne Issyk-Kul virus in Europe.

Bats have been gaining attention as potential reservoir hosts of numerous viruses pathogenic to animals and man. Issyk-Kul virus, a member of the family Nairoviridae, was first isolated in the 1970s from vespertilionid bats in Central Asia. Issyk-Kul virus has been described as human-pathogenic virus, causing febrile outbreaks in humans with headaches, myalgia and nausea. Here we describe the detection of a novel strain of Issyk-Kul virus from Eptesicus nilssonii in Germany. This finding indicates for the first time the prevalence of these zoonotic viruses in Europe.

Zwiesel bat banyangvirus, a potentially zoonotic Huaiyangshan banyangvirus (Formerly known as SFTS)-like banyangvirus in Northern bats from Germany.

Bats are reservoir hosts for several emerging and re-emerging viral pathogens causing morbidity and mortality in wildlife, animal stocks and humans. Various viruses within the family Phenuiviridae have been detected in bats, including the highly pathogenic Rift Valley fever virus and Malsoor virus, a novel Banyangvirus with close genetic relation to Huaiyangshan banyangvirus (BHAV)(former known as Severe fever with thrombocytopenia syndrome virus, SFTSV) and Heartland virus (HRTV), both of which have caused severe disease with fatal casualties in humans. In this study we present the whole genome of a novel Banyangvirus, named Zwiesel bat banyangvirus, revealed through deep sequencing of the Eptesicus nilssonii bat virome. The detection of the novel bat banyangvirus, which is in close phylogenetic relationship with the pathogenic HRTV and BHAV, underlines the possible impact of emerging phenuiviruses on public health.

A novel Coltivirus-related virus isolated from free-tailed bats from Côte d'Ivoire is able to infect human cells in vitro.

BACKGROUND: Zoonotic transmission events play a major role in the emergence of novel diseases. While such events are virtually impossible to predict, wildlife screening for potential emerging pathogens can be a first step. Driven by recent disease epidemics like severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS), and Ebola, bats have gained special interest as reservoirs of emerging viruses. METHODS: As part of a bigger study investigating pathogens in African bats we screened animals for the presence of known and unknown viruses. RESULTS: We isolated and characterised a novel reovirus from blood of free-tailed bats (Chaereophon aloysiisabaudiae) captured in 2006 in Côte d'Ivoire. The virus showed closest relationship with two human pathogenic viruses, Colorado tick fever virus and Eyach virus, and was able to infect various human cell lines in vitro. CONCLUSION: The study shows the presence of a coltivirus-related virus in bats from Sub-Sahara Africa. Serological studies could help to assess its impact on humans or wildlife health.

Comparing Viral Metagenomic Extraction Methods.

A crucial step in the molecular detection of viruses in clinical specimens is the efficient extraction of viral nucleic acids. The total yield of viral nucleic acid from a clinical specimen is dependent on the specimen's volume, the initial virus concentration and the effectiveness provided by the extraction method. Recent Next Generation Sequencing (NGS)-based diagnostic approaches (i.e. metagenomics) provide a molecular 'open view' into the sample, as they theoretically generate sequence reads of any nucleic acid present in a specimen in a statistically representative manner. However, since a higher virus-related read output promises better sensitivity in the subsequent bioinformatic analysis, the extraction method selected determines the reliability of diagnostic NGS. In this study nine commercially available kits for nucleic acid extraction were compared regarding the simultaneous isolation of DNA and RNA by real-time PCR,four of which were selected for subsequent comparison by NGS (QIAamp Viral RNA Mini Kit, QIAamp DNA Blood Mini Kit, QIAamp cador Pathogen Mini Kit and QIAamp MinElute Virus Spin Kit). The nucleic acid yields and the sequence read output were compared for four different model viruses comprising Reovirus, Orthomyxovirus, Orthopoxvirus and Paramyxovirus, each at defined but varying concentrations in the same sample. The total amount of nucleic acid was processed to sequence the RNA (as cDNA) and the DNA with quantification by Qubit and virus-specific quantitative real-time PCRs. NGS libraries were prepared for sequencing on the Illumina HiSeq 1500 system. Finally, the percentage of reads assignable to each virus was determined via mapping. Evaluation of different commercial nucleic acid extraction kits with four different viruses indicates little variation in the read numbers obtained for transcribed RNA or DNA by NGS. Since NGSis increasingly being used as a tool in diagnostics of infectious diseases, the individual steps of the complete process have to be validated carefully. Here we could show that for virus identification in liquid clinical specimens, any nucleic acid extraction kit that is performing well for PCR diagnostics can be used for NGS diagnostics as well and that the selection of the kit has only a minor impact on the yield of viral reads.

Comparison of 454 Ultra-Deep Sequencing and Allele-Specific Real-Time PCR with Regard to the Detection of Emerging Drug-Resistant Minor HIV-1 Variants after Antiretroviral Prophylaxis for Vertical Transmission.

BACKGROUND: Pregnant HIV-infected women were screened for the development of HIV-1 drug resistance after implementation of a triple-antiretroviral transmission prophylaxis as recommended by the WHO in 2006. The study offered the opportunity to compare amplicon-based 454 ultra-deep sequencing (UDS) and allele-specific real-time PCR (ASPCR) for the detection of drug-resistant minor variants in the HIV-1 reverse transcriptase (RT). METHODS: Plasma samples from 34 Tanzanian women were previously analysed by ASPCR for key resistance mutations in the viral RT selected by AZT, 3TC, and NVP (K70R, K103N, Y181C, M184V, T215Y/F). In this study, the RT region of the same samples was investigated by amplicon-based UDS for resistance mutations using the 454 GS FLX System. RESULTS: Drug-resistant HIV-variants were identified in 69% (20/29) of women by UDS and in 45% (13/29) by ASPCR. The absolute number of resistance mutations identified by UDS was twice that identified by ASPCR (45 vs 24). By UDS 14 of 24 ASPCR-detected resistance mutations were identified at the same position. The overall concordance between UDS and ASPCR was 61.0% (25/41). The proportions of variants quantified by UDS were approximately 2-3 times lower than by ASPCR. Amplicon generation from samples with viral loads below 20,000 copies/ml failed more frequently by UDS compared to ASPCR (limit of detection = 650 copies/ml), resulting in missing or insufficient sequence coverage. CONCLUSIONS: Both methods can provide useful information about drug-resistant minor HIV-1 variants. ASPCR has a higher sensitivity than UDS, but is restricted to single resistance mutations. In contrast, UDS is limited by its requirement for high viral loads to achieve sufficient sequence coverage, but the sequence information reveals the complete resistance patterns within the genomic region analysed. Improvements to the UDS limit of detection are in progress, and UDS could then facilitate monitoring of drug-resistant minor variants in the HIV-1 quasispecies.

SuRankCo: supervised ranking of contigs in de novo assemblies.

BACKGROUND: Evaluating the quality and reliability of a de novo assembly and of single contigs in particular is challenging since commonly a ground truth is not readily available and numerous factors may influence results. Currently available procedures provide assembly scores but lack a comparative quality ranking of contigs within an assembly. RESULTS: We present SuRankCo, which relies on a machine learning approach to predict quality scores for contigs and to enable the ranking of contigs within an assembly. The result is a sorted contig set which allows selective contig usage in downstream analysis. Benchmarking on datasets with known ground truth shows promising sensitivity and specificity and favorable comparison to existing methodology. CONCLUSIONS: SuRankCo analyzes the reliability of de novo assemblies on the contig level and thereby allows quality control and ranking prior to further downstream and validation experiments.

What caused the outbreak of ESBL-producing Klebsiella pneumoniae in a neonatal intensive care unit, Germany 2009 to 2012? Reconstructing transmission with epidemiological analysis and whole-genome sequencing.

OBJECTIVE: We aimed to retrospectively reconstruct the timing of transmission events and pathways in order to understand why extensive preventive measures and investigations were not sufficient to prevent new cases. METHODS: We extracted available information from patient charts to describe cases and to compare them to the normal population of the ward. We conducted a cohort study to identify risk factors for pathogen acquisition. We sequenced the available isolates to determine the phylogenetic relatedness of Klebsiella pneumoniae isolates on the basis of their genome sequences. RESULTS: The investigation comprises 37 cases and the 10 cases with ESBL (extended-spectrum beta-lactamase)-producing K. pneumoniae bloodstream infection. Descriptive epidemiology indicated that a continuous transmission from person to person was most likely. Results from the cohort study showed that 'frequent manipulation' (a proxy for increased exposure to medical procedures) was significantly associated with being a case (RR 1.44, 95% CI 1.02 to 2.19). Genome sequences revealed that all 48 bacterial isolates available for sequencing from 31 cases were closely related (maximum genetic distance, 12 single nucleotide polymorphisms). Based on our calculation of evolutionary rate and sequence diversity, we estimate that the outbreak strain was endemic since 2008. CONCLUSIONS: Epidemiological and phylogenetic analyses consistently indicated that there were additional, undiscovered cases prior to the onset of microbiological screening and that the spread of the pathogen remained undetected over several years, driven predominantly by person-to-person transmission. Whole-genome sequencing provided valuable information on the onset, course and size of the outbreak, and on possible ways of transmission.

Pipasic: similarity and expression correction for strain-level identification and quantification in metaproteomics.

MOTIVATION: Metaproteomic analysis allows studying the interplay of organisms or functional groups and has become increasingly popular also for diagnostic purposes. However, difficulties arise owing to the high sequence similarity between related organisms. Further, the state of conservation of proteins between species can be correlated with their expression level, which can lead to significant bias in results and interpretation. These challenges are similar but not identical to the challenges arising in the analysis of metagenomic samples and require specific solutions. RESULTS: We introduce Pipasic (peptide intensity-weighted proteome abundance similarity correction) as a tool that corrects identification and spectral counting-based quantification results using peptide similarity estimation and expression level weighting within a non-negative lasso framework. Pipasic has distinct advantages over approaches only regarding unique peptides or aggregating results to the lowest common ancestor, as demonstrated on examples of viral diagnostics and an acid mine drainage dataset. AVAILABILITY AND IMPLEMENTATION: Pipasic source code is freely available from https://sourceforge.net/projects/pipasic/. CONTACT: RenardB@rki.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Fatal monkeypox in wild-living sooty mangabey, Côte d'Ivoire, 2012.

We isolated a monkeypox virus from a wild-living monkey, a sooty mangabey, found dead in Taï National Park, Côte d'Ivoire, in March 2012. The whole-genome sequence obtained from this isolate and directly from clinical specimens showed its close relationship to monkeypox viruses from Western Africa.

Genome-wide comparison of cowpox viruses reveals a new clade related to Variola virus.

Zoonotic infections caused by several orthopoxviruses (OPV) like monkeypox virus or vaccinia virus have a significant impact on human health. In Europe, the number of diagnosed infections with cowpox viruses (CPXV) is increasing in animals as well as in humans. CPXV used to be enzootic in cattle; however, such infections were not being diagnosed over the last decades. Instead, individual cases of cowpox are being found in cats or exotic zoo animals that transmit the infection to humans. Both animals and humans reveal local exanthema on arms and legs or on the face. Although cowpox is generally regarded as a self-limiting disease, immunosuppressed patients can develop a lethal systemic disease resembling smallpox. To date, only limited information on the complex and, compared to other OPV, sparsely conserved CPXV genomes is available. Since CPXV displays the widest host range of all OPV known, it seems important to comprehend the genetic repertoire of CPXV which in turn may help elucidate specific mechanisms of CPXV pathogenesis and origin. Therefore, 22 genomes of independent CPXV strains from clinical cases, involving ten humans, four rats, two cats, two jaguarundis, one beaver, one elephant, one marah and one mongoose, were sequenced by using massive parallel pyrosequencing. The extensive phylogenetic analysis showed that the CPXV strains sequenced clearly cluster into several distinct clades, some of which are closely related to Vaccinia viruses while others represent different clades in a CPXV cluster. Particularly one CPXV clade is more closely related to Camelpox virus, Taterapox virus and Variola virus than to any other known OPV. These results support and extend recent data from other groups who postulate that CPXV does not form a monophyletic clade and should be divided into multiple lineages.

MultiPSQ: a software solution for the analysis of diagnostic n-plexed pyrosequencing reactions.

BACKGROUND: Pyrosequencing can be applied for Single-Nucleotide-Polymorphism (SNP)-based pathogen typing or for providing sequence information of short DNA stretches. However, for some pathogens molecular typing cannot be performed relying on a single SNP or short sequence stretch, necessitating the consideration of several genomic regions. A promising rapid approach is the simultaneous application of multiple sequencing primers, called multiplex pyrosequencing. These primers generate a fingerprint-pyrogram which is constituted by the sum of all individual pyrograms originating from each primer used. METHODS: To improve pyrosequencing-based pathogen typing, we have developed the software tool MultiPSQ that expedites the analysis and evaluation of multiplex-pyrograms. As a proof of concept, a multiplex pyrosequencing assay for the typing of orthopoxviruses was developed to analyse clinical samples diagnosed in the German Consultant Laboratory for Poxviruses. RESULTS: The software tool MultiPSQ enabled the analysis of multiplex-pyrograms originating from various pyrosequencing primers. Thus several target regions can be used for pathogen typing based on pyrosequencing. As shown with a proof of concept assay, SNPs present in different orthopoxvirus strains could be identified correctly with two primers by MultiPSQ. CONCLUSIONS: Software currently available is restricted to a fixed number of SNPs and sequencing primers, severely limiting the usefulness of this technique. In contrast, our novel software MultiPSQ allows analysis of data from multiplex pyrosequencing assays that contain any number of sequencing primers covering any number of polymorphisms.

Comparison of host cell gene expression in cowpox, monkeypox or vaccinia virus-infected cells reveals virus-specific regulation of immune response genes.

BACKGROUND: Animal-borne orthopoxviruses, like monkeypox, vaccinia and the closely related cowpox virus, are all capable of causing zoonotic infections in humans, representing a potential threat to human health. The disease caused by each virus differs in terms of symptoms and severity, but little is yet know about the reasons for these varying phenotypes. They may be explained by the unique repertoire of immune and host cell modulating factors encoded by each virus. In this study, we analysed the specific modulation of the host cell's gene expression profile by cowpox, monkeypox and vaccinia virus infection. We aimed to identify mechanisms that are either common to orthopoxvirus infection or specific to certain orthopoxvirus species, allowing a more detailed description of differences in virus-host cell interactions between individual orthopoxviruses. To this end, we analysed changes in host cell gene expression of HeLa cells in response to infection with cowpox, monkeypox and vaccinia virus, using whole-genome gene expression microarrays, and compared these to each other and to non-infected cells. RESULTS: Despite a dominating non-responsiveness of cellular transcription towards orthopoxvirus infection, we could identify several clusters of infection-modulated genes. These clusters are either commonly regulated by orthopoxvirus infection or are uniquely regulated by infection with a specific orthopoxvirus, with major differences being observed in immune response genes. Most noticeable was an induction of genes involved in leukocyte migration and activation in cowpox and monkeypox virus-infected cells, which was not observed following vaccinia virus infection. CONCLUSION: Despite their close genetic relationship, the expression profiles induced by infection with different orthopoxviruses vary significantly. It may be speculated that these differences at the cellular level contribute to the individual characteristics of cowpox, monkeypox and vaccinia virus infections in certain host species.

mPSQed: a software for the design of multiplex pyrosequencing assays.

Molecular-based diagnostic assays are the gold standard for infectious diseases today, since they allow a rapid and sensitive identification and typing of various pathogens. While PCR can be designed to be specific for a certain pathogen, a subsequent sequence analysis is frequently required for confirmation or typing. The design of appropriate PCR-based assays is a complex task, especially when conserved discriminating polymorphisms are rare or if the number of types which need to be differentiated is high. One extremely useful but underused method for this purpose is the multiplex pyrosequencing technique. Unfortunately there is no software available to aid researchers in designing multiplex pyrosequencing assays. Here, we present mPSQed (Multiplex PyroSeQuencing EDitor), a program targeted at closing this gap. We also present the design of an exemplarily theoretical assay for the differentiation of human adenovirus types A-F using two pyrosequencing primers on two distinct PCR products, designed quickly and easily using our software.

Genome analysis of bat adenovirus 2: indications of interspecies transmission.

The genome of bat adenovirus 2 was sequenced and analyzed. It is similar in size (31,616 bp) to the genomes of bat adenovirus 3 and canine adenoviruses 1 and 2. These four viruses are monophyletic and share an identical genome organization, with one E3 gene and four E4 genes unique to this group among the mastadenoviruses. These findings suggest that canine adenoviruses may have originated by interspecies transfer of a vespertilionid bat adenovirus.