Optimization and validation of sample preparation for metagenomic sequencing of viruses in clinical samples

Comparison of the performance of two targeted metagenomic virus capture probe-based methods using reference control materials and clinical samples

Viral enrichment by probe hybridization has been reported to significantly increase the sensitivity of viral metagenomics. This study compares the analytical performance of two targeted metagenomic virus capture probe-based methods: (i) SeqCap EZ HyperCap by Roche (ViroCap) and (ii) Twist Comprehensive Viral Research Panel workflow, for diagnostic use. Sensitivity, specificity, and limit of detection were analyzed using 25 synthetic viral sequences spiked in increasing proportions of human background DNA, eight clinical samples, and American Type Culture Collection (ATCC) Virome Virus Mix. Sensitivity and specificity were 95% and higher for both methods using the synthetic and reference controls as gold standard. Combining thresholds for viral sequence read counts and genome coverage [respectively 500 reads per million (RPM) and 10% coverage] resulted in optimal prediction of true positive results. Limits of detection were approximately 50-500 copies/mL for both methods as determined by ddPCR. Increasing proportions of spike-in cell-free human background sequences up to 99.999% (50 ng/mL) did not negatively affect viral detection, suggesting effective capture of viral sequences. These data show analytical performances in ranges applicable to clinical samples, for both probe hybridization metagenomic approaches. This study supports further steps toward more widespread use of viral metagenomics for pathogen detection, in clinical and surveillance settings using low biomass samples.

IMPORTANCE

Viral metagenomics has been gradually applied for broad-spectrum pathogen detection of infectious diseases, surveillance of emerging diseases, and pathogen discovery. Viral enrichment by probe hybridization methods has been reported to significantly increase the sensitivity of viral metagenomics. During the past years, a specific hybridization panel distributed by Roche has been adopted in a broad range of different clinical and zoonotic settings. Recently, Twist Bioscience has released a new hybridization panel targeting human and animal viruses. This is the first report comparing the performance of viral metagenomic hybridization panels.

Retrospective Genotyping of Enteroviruses Using a Diagnostic Nanopore Sequencing Workflow

Enteroviruses are among the most common viruses pathogenic to humans. They are associated with various forms of disease, ranging from mild respiratory illness to severe neurological diseases. In recent years, an increasing number of isolated cases of children developing meningitis or encephalitis as a result of enterovirus infection have been reported, as well as discrete enterovirus D68 outbreaks in North America in 2014 and 2016. We developed an assay to rapidly genotype enteroviruses by sequencing a region within the VP1 gene using nanopore Flongles. We retrospectively analyzed enterovirus-/rhinovirus-positive clinical samples from the Zurich, Switzerland area mainly collected during two seasons in 2019/2020 and 2021/2022. Respiratory, cerebrospinal fluid, and stool samples were analyzed. Whole-genome sequencing was performed on samples with ambiguous genotyping results and enterovirus D68-positive samples. Out of 255 isolates, a total of 95 different genotypes were found. A difference in the prevalence of enterovirus and rhinovirus infections was observed for both sample type and age group. In particular, children aged 0-4 years showed a higher frequency of enterovirus infections. Comparing the respiratory seasons, a higher prevalence was found, especially for enterovirus A and rhinovirus A after the SARS-CoV-2 pandemic. The enterovirus genotyping workflow provides a rapid diagnostic tool for individual analysis and continuous enterovirus surveillance.

Next-Generation Sequencing for the Detection of Microbial Agents in Avian Clinical Samples

Direct-targeted next-generation sequencing (tNGS), with its undoubtedly superior diagnostic capacity over real-time PCR (RT-PCR), and direct-non-targeted NGS (ntNGS), with its higher capacity to identify and characterize multiple agents, are both likely to become diagnostic methods of choice in the future. tNGS is a rapid and sensitive method for precise characterization of suspected agents. ntNGS, also known as agnostic diagnosis, does not require a hypothesis and has been used to identify unsuspected infections in clinical samples. Implemented in the form of multiplexed total DNA metagenomics or as total RNA sequencing, the approach produces comprehensive and actionable reports that allow semi-quantitative identification of most of the agents present in respiratory, cloacal, and tissue samples. The diagnostic benefits of the use of direct tNGS and ntNGS are high specificity, compatibility with different types of clinical samples (fresh, frozen, FTA cards, and paraffin-embedded), production of nearly complete infection profiles (viruses, bacteria, fungus, and parasites), production of "semi-quantitative" information, direct agent genotyping, and infectious agent mutational information. The achievements of NGS in terms of diagnosing poultry problems are described here, along with future applications. Multiplexing, development of standard operating procedures, robotics, sequencing kits, automated bioinformatics, cloud computing, and artificial intelligence (AI) are disciplines converging toward the use of this technology for active surveillance in poultry farms. Other advances in human and veterinary NGS sequencing are likely to be adaptable to avian species in the future.

Rapid and sensitive single-sample viral metagenomics using Nanopore Flongle sequencing

The ability of viral metagenomic Next-Generation Sequencing (mNGS) to unbiasedly detect nucleic acids in a clinical sample is a powerful tool for advanced diagnosis of viral infections. When clinical symptoms do not provide a clear differential diagnosis, extensive laboratory testing with virus-specific PCR and serology can be replaced by a single viral mNGS analysis. However, widespread diagnostic use of viral mNGS is thus far limited by long sample-to-result times, as most protocols rely on Illumina sequencing, which provides high and accurate sequencing output but is time-consuming and expensive. Here, we describe the development of an mNGS protocol based on the more cost-effective Nanopore Flongle sequencing with decreased turnaround time and lower, yet sufficient sequencing output to provide sensitive virus detection. Sample preparation (6 h) and sequencing (2 h) times are substantially reduced compared to Illumina mNGS and allow detection of DNA/RNA viruses at low input (up to 33-38 cycle threshold of specific qPCR). Although Flongles yield lower sequencing output, direct comparison with Illumina mNGS on diverse clinical samples showed similar results. Collectively, the novel Nanopore mNGS approach is specifically tailored for use in clinical diagnostics and provides a rapid and cost-effective mNGS strategy for individual testing of severe cases.

The International Virus Bioinformatics Meeting 2023

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.

Evaluation of extraction and enrichment methods for recovery of respiratory RNA viruses in a metagenomics approach

Viral metagenomics is increasingly applied in viral detection and virome characterization. Different extraction and enrichment techniques may be adopted, however, reports on their effective influence on viral recovery is often conflicting. Using a three step enrichment steps, the effect of three extraction kits and the influence of DNase treatment with or without rRNA removal for respiratory RNA virus recovery from nasopharyngeal swab samples was evaluated. The viral cocktail containing six different RNA viruses pooled in equal volume were subjected to the different extraction and enrichment methods, sequenced using the Illumina MiSeq, and analysed using Genome Detective. The PureLink® Viral RNA/DNA Mini Kit (PureLink) was highly efficient with better recovery of all the viral agents in the cocktail. The use of rRNA treatment resulted in increased viral recovery with PureLink and QIAamp® Viral RNA Mini kit, while having comparable recovery rate as DNase only with the QIAamp® MinElute Virus Spin Kit. The observed low reads and genome coverage of some of the viruses could be attributed to their low abundance. Depending on sample matrix, extraction choice and enrichment strategy may influence recovery of respiratory RNA virus in metagenomics studies, therefore individual evaluation and adoption may be necessary for a robust result.

Evolution and emergence of mosquito-borne viruses of medical importance: towards a routine metagenomic surveillance approach

During the last two decades, the world has witnessed the emergence and re-emergence of arthropod-borne viruses, better known as arboviruses. The close contact between sylvatic, rural and peri-urban vector species and humans has been mainly determined by the environment-modifying human activity. The resulting interactions have led to multiple dead-end host infections and have allowed sylvatic arboviruses to eventually adapt to new vectors and hosts, contributing to the establishment of urban transmission cycles of some viruses with enormous epidemiologic impact. The metagenomic next-generation sequencing (NGS) approach has allowed obtaining unbiased sequence information of millions of DNA and RNA molecules from clinical and environmental samples. Robust bioinformatics tools have enabled the assembly of individual sequence reads into contigs and scaffolds partially or completely representing the genomes of the microorganisms and viruses being present in biological samples of clinical relevance. In this review, we describe the different ecological scenarios for the emergence of viral diseases, the virus adaptation process required for the establishment of a new transmission cycle and the usefulness of NGS and computational methods for the discovery and routine genomic surveillance of mosquito-borne viruses in their ecosystems.

Critical Assessment of Whole Genome and Viral Enrichment Shotgun Metagenome on the Characterization of Stool Total Virome in Hepatocellular Carcinoma Patients

Viruses are the most abundant form of life on earth and play important roles in a broad range of ecosystems. Currently, two methods, whole genome shotgun metagenome (WGSM) and viral-like particle enriched metagenome (VLPM) sequencing, are widely applied to compare viruses in various environments. However, there is no critical assessment of their performance in recovering viruses and biological interpretation in comparative viral metagenomic studies. To fill this gap, we applied the two methods to investigate the stool virome in hepatocellular carcinoma (HCC) patients and healthy controls. Both WGSM and VLPM methods can capture the major diversity patterns of alpha and beta diversities and identify the altered viral profiles in the HCC stool samples compared with healthy controls. Viral signatures identified by both methods showed reductions of Faecalibacterium virus Taranis in HCC patients' stool. Ultra-deep sequencing recovered more viruses in both methods, however, generally, 3 or 5 Gb were sufficient to capture the non-fragmented long viral contigs. More lytic viruses were detected than lysogenetic viruses in both methods, and the VLPM can detect the RNA viruses. Using both methods would identify shared and specific viral signatures and would capture different parts of the total virome.

Next Generation Sequencing Approaches to Characterize the Respiratory Tract Virome

The COVID-19 pandemic and heightened perception of the risk of emerging viral infections have boosted the efforts to better understand the virome or complete repertoire of viruses in health and disease, with a focus on infectious respiratory diseases. Next-generation sequencing (NGS) is widely used to study microorganisms, allowing the elucidation of bacteria and viruses inhabiting different body systems and identifying new pathogens. However, NGS studies suffer from a lack of standardization, in particular, due to various methodological approaches and no single format for processing the results. Here, we review the main methodological approaches and key stages for studies of the human virome, with an emphasis on virome changes during acute respiratory viral infection, with applications for clinical diagnostics and epidemiologic analyses.

Development and implementation of a customised rapid syndromic diagnostic test for severe pneumonia

Background: The diagnosis of pneumonia has been hampered by a reliance on bacterial cultures which take several days to return a result, and are frequently negative. In critically ill patients this leads to the use of empiric, broad-spectrum antimicrobials and compromises good antimicrobial stewardship. The objective of this study was to establish the performance of a syndromic molecular diagnostic approach, using a custom TaqMan array card (TAC) covering 52 respiratory pathogens, and assess its impact on antimicrobial prescribing. Methods: The TAC was validated against a retrospective multi-centre cohort of broncho-alveolar lavage samples. The TAC was assessed prospectively in patients undergoing investigation for suspected pneumonia, with a comparator cohort formed of patients investigated when the TAC laboratory team were unavailable. Co-primary outcomes were sensitivity compared to conventional microbiology and, for the prospective study, time to result. Metagenomic sequencing was performed to validate findings in prospective samples. Antibiotic free days (AFD) were compared between the study cohort and comparator group. Results: 128 stored samples were tested, with sensitivity of 97% (95% confidence interval (CI) 88-100%). Prospectively, 95 patients were tested by TAC, with 71 forming the comparator group. TAC returned results 51 hours (interquartile range 41-69 hours) faster than culture and with sensitivity of 92% (95% CI 83-98%) compared to conventional microbiology. 94% of organisms identified by sequencing were detected by TAC. There was a significant difference in the distribution of AFDs with more AFDs in the TAC group (p=0.02). TAC group were more likely to experience antimicrobial de-escalation (odds ratio 2.9 (95%1.5-5.5)). Conclusions: Implementation of a syndromic molecular diagnostic approach to pneumonia led to faster results, with high sensitivity and impact on antibiotic prescribing.

Development and implementation of a customised rapid syndromic diagnostic test for severe pneumonia

Background: The diagnosis of pneumonia has been hampered by a reliance on bacterial cultures which take several days to return a result, and are frequently negative. In critically ill patients this leads to the use of empiric, broad-spectrum antimicrobials and compromises good antimicrobial stewardship. The objective of this study was to establish the performance of a syndromic molecular diagnostic approach, using a custom TaqMan array card (TAC) covering 52 respiratory pathogens, and assess its impact on antimicrobial prescribing. Methods: The TAC was validated against a retrospective multi-centre cohort of broncho-alveolar lavage samples. The TAC was assessed prospectively in patients undergoing investigation for suspected pneumonia, with a comparator cohort formed of patients investigated when the TAC laboratory team were unavailable. Co-primary outcomes were sensitivity compared to conventional microbiology and, for the prospective study, time to result. Metagenomic sequencing was performed to validate findings in prospective samples. Antibiotic free days (AFD) were compared between the study cohort and comparator group. Results: 128 stored samples were tested, with sensitivity of 97% (95% confidence interval (CI) 88-100%). Prospectively, 95 patients were tested by TAC, with 71 forming the comparator group. TAC returned results 51 hours (interquartile range 41-69 hours) faster than culture and with sensitivity of 92% (95% CI 83-98%) compared to conventional microbiology. 94% of organisms identified by sequencing were detected by TAC. There was a significant difference in the distribution of AFDs with more AFDs in the TAC group (p=0.02). TAC group were more likely to experience antimicrobial de-escalation (odds ratio 2.9 (95%1.5-5.5)). Conclusions: Implementation of a syndromic molecular diagnostic approach to pneumonia led to faster results, with high sensitivity and impact on antibiotic prescribing.

Metagenomic sequencing reveals time, host, and body compartment-specific viral dynamics after lung transplantation

BACKGROUND

The virome of lung transplant recipients (LTRs) under immunosuppressive therapy is dominated by non-pathogenic Anelloviridae and further includes several pathogenic viruses such as Herpesviruses or respiratory viruses. It is unclear whether the donor-derived virome in the transplanted lung influences recipient virome dynamics in other body compartments and if so, to which degree. Likewise, it is unknown whether dependencies exist among virus populations that mutually shape viral loads and kinetics.

RESULTS

To address these questions, we characterized viral communities in airways and plasma of 49 LTRs and analyzed their abundance patterns in a data modeling approach. We found distinct viral clusters that were specific for body compartments and displayed independent dynamics. These clusters robustly gathered specific viral species across the patient cohort. In the lung, viral cluster abundance associated with time after transplantation and we detected mutual exclusion of viral species within the same human host. In plasma, viral cluster dynamics were associated with the indication for transplantation lacking significant short-time changes. Interestingly, pathogenic viruses in the plasma co-occurred specifically with Alpha torque virus genogroup 4 and Gamma torque virus strains suggesting shared functional or ecological requirements.

CONCLUSIONS

In summary, the detailed analysis of virome dynamics after lung transplantation revealed host, body compartment, and time-specific dependency patterns among viruses. Furthermore, our results suggested genetic adaptation to the host microenvironment at the level of the virome and support the hypothesis of functional complementarity between Anellovirus groups and other persistent viruses. Video abstract.

Detection of Low-Copy Human Virus DNA upon Prolonged Formalin Fixation

Formalin fixation, albeit an outstanding method for morphological and molecular preservation, induces DNA damage and cross-linking, which can hinder nucleic acid screening. This is of particular concern in the detection of low-abundance targets, such as persistent DNA viruses. In the present study, we evaluated the analytical sensitivity of viral detection in lung, liver, and kidney specimens from four deceased individuals. The samples were either frozen or incubated in formalin (±paraffin embedding) for up to 10 days. We tested two DNA extraction protocols for the control of efficient yields and viral detections. We used short-amplicon qPCRs (63–159 nucleotides) to detect 11 DNA viruses, as well as hybridization capture of these plus 27 additional ones, followed by deep sequencing. We observed marginally higher ratios of amplifiable DNA and scantly higher viral genoprevalences in the samples extracted with the FFPE dedicated protocol. Based on the findings in the frozen samples, most viruses were detected regardless of the extended fixation times. False-negative calls, particularly by qPCR, correlated with low levels of viral DNA (<250 copies/million cells) and longer PCR amplicons (>150 base pairs). Our data suggest that low-copy viral DNAs can be satisfactorily investigated from FFPE specimens, and encourages further examination of historical materials.

Development and implementation of a customised rapid syndromic diagnostic test for severe pneumonia

Background: The diagnosis of pneumonia has been hampered by a reliance on bacterial cultures which take several days to return a result, and are frequently negative. In critically ill patients this leads to the use of empiric, broad-spectrum antimicrobials and compromises good antimicrobial stewardship. The objective of this study was to establish the performance of a syndromic molecular diagnostic approach, using a custom TaqMan array card (TAC) covering 52 respiratory pathogens, and assess its impact on antimicrobial prescribing. Methods: The TAC was validated against a retrospective multi-centre cohort of broncho-alveolar lavage samples. The TAC was assessed prospectively in patients undergoing investigation for suspected pneumonia, with a comparator cohort formed of patients investigated when the TAC laboratory team were unavailable. Co-primary outcomes were sensitivity compared to conventional microbiology and, for the prospective study, time to result. Metagenomic sequencing was performed to validate findings in prospective samples. Antibiotic free days (AFD) were compared between the study cohort and comparator group. Results: 128 stored samples were tested, with sensitivity of 97% (95% confidence interval (CI) 88-100%). Prospectively, 95 patients were tested by TAC, with 71 forming the comparator group. TAC returned results 51 hours (interquartile range 41-69 hours) faster than culture and with sensitivity of 92% (95% CI 83-98%) compared to conventional microbiology. 94% of organisms identified by sequencing were detected by TAC. There was a significant difference in the distribution of AFDs with more AFDs in the TAC group (p=0.02). TAC group were more likely to experience antimicrobial de-escalation (odds ratio 2.9 (95%1.5-5.5)). Conclusions: Implementation of a syndromic molecular diagnostic approach to pneumonia led to faster results, with high sensitivity and impact on antibiotic prescribing.

Enrichment approach for unbiased sequencing of respiratory syncytial virus directly from clinical samples

Background: Nasopharyngeal samples contain higher quantities of bacterial and host nucleic acids relative to viruses; presenting challenges during virus metagenomics sequencing, which underpins agnostic sequencing protocols. We aimed to develop a viral enrichment protocol for unbiased whole-genome sequencing of respiratory syncytial virus (RSV) from nasopharyngeal samples using the Oxford Nanopore Technology (ONT) MinION platform. Methods: We assessed two protocols using RSV positive samples. Protocol 1 involved physical pre-treatment of samples by centrifugal processing before RNA extraction, while Protocol 2 entailed direct RNA extraction without prior enrichment. Concentrates from Protocol 1 and RNA extracts from Protocol 2 were each divided into two fractions; one was DNase treated while the other was not. RNA was then extracted from both concentrate fractions per sample and RNA from both protocols converted to cDNA, which was then amplified using the tagged Endoh primers through Sequence-Independent Single-Primer Amplification (SISPA) approach, a library prepared, and sequencing done. Statistical significance during analysis was tested using the Wilcoxon signed-rank test. Results: DNase-treated fractions from both protocols recorded significantly reduced host and bacterial contamination unlike the untreated fractions (in each protocol p<0.01). Additionally, DNase treatment after RNA extraction (Protocol 2) enhanced host and bacterial read reduction compared to when done before (Protocol 1). However, neither protocol yielded whole RSV genomes. Sequenced reads mapped to parts of the nucleoprotein (N gene) and polymerase complex (L gene) from Protocol 1 and 2, respectively. Conclusions: DNase treatment was most effective in reducing host and bacterial contamination, but its effectiveness improved if done after RNA extraction than before. We attribute the incomplete genome segments to amplification biases resulting from the use of short length random sequence (6 bases) in tagged Endoh primers. Increasing the length of the random nucleotides from six hexamers to nine or 12 in future studies may reduce the coverage biases.

Deep Impact of Random Amplification and Library Construction Methods on Viral Metagenomics Results

Clinical metagenomics is a broad-range agnostic detection method of pathogens, including novel microorganisms. A major limit is the low pathogen load compared to the high background of host nucleic acids. To overcome this issue, several solutions exist, such as applying a very high depth of sequencing, or performing a relative enrichment of viral genomes associated with capsids. At the end, the quantity of total nucleic acids is often below the concentrations recommended by the manufacturers of library kits, which necessitates to random amplify nucleic acids. Using a pool of 26 viruses representative of viral diversity, we observed a deep impact of the nature of sample (total nucleic acids versus RNA only), the reverse transcription, the random amplification and library construction method on virus recovery. We further optimized the two most promising methods and assessed their performance with fully characterized reference virus stocks. Good genome coverage and limit of detection lower than 100 or 1000 genome copies per mL of plasma, depending on the genome viral type, were obtained from a three million reads dataset. Our study reveals that optimized random amplification is a technique of choice when insufficient amounts of nucleic acid are available for direct libraries constructions.

Biomimetic Virus-Like Particles as Severe Acute Respiratory Syndrome Coronavirus 2 Diagnostic Tools

Coronavirus disease 2019 (COVID-19) is a highly transmissible disease that has affected more than 90% of the countries worldwide. At least 17 million individuals have been infected, and some countries are still battling first or second waves of the pandemic. Nucleic acid tests, especially reverse transcription polymerase chain reaction (RT-PCR), have become the workhorse for early detection of COVID-19 infection. Positive controls for the molecular assays have been developed to validate each test and to provide high accuracy. However, most available positive controls require cold-chain distribution and cannot serve as full-process control. To overcome these shortcomings, we report the production of biomimetic virus-like particles (VLPs) as SARS-CoV-2 positive controls. A SARS-CoV-2 detection module for RT-PCR was encapsidated into VLPs from a bacteriophage and a plant virus. The chimeric VLPs were obtained either by in vivo reconstitution and coexpression of the target detection module and coat proteins or by in vitro assembly of purified detection module RNA sequences and coat proteins. These VLP-based positive controls mimic SARS-CoV-2 packaged ribonucleic acid (RNA) while being noninfectious. Most importantly, we demonstrated that the positive controls are scalable, stable, and can serve broadly as controls, from RNA extraction to PCR in clinical settings.

Implementation of next-generation sequencing for virus identification in veterinary diagnostic laboratories

The value of next-generation sequencing (NGS)-based applications for testing purposes in human medicine is widely recognized. Although NGS-based metagenomic screening may be of interest in veterinary medicine, in particular for intensively farmed livestock species such as pigs, there is a lack of protocols tailored to veterinary requirements, likely because of the high diversity of species and samples. Therefore, we developed an NGS-based protocol for use in veterinary virology and present here different applications in porcine medicine. To develop the protocol, each step of sample preparation was optimized using porcine samples spiked with various RNA and DNA viruses. The resulting protocol was tested with clinical samples previously confirmed to be positive for specific viruses by a diagnostic laboratory. Additionally, we validated the protocol in an NGS viral metagenomics ring trial and tested the protocol on viral multiplex reference material (NIBSC, U.K.). We applied our ViroScreen protocol successfully for 1) virus identification, 2) virus characterization, and 3) herd screening. We identified torque teno sus virus and atypical porcine pestivirus in a neurologic case, determined the full-length genome sequence of swine influenza A virus in field samples, and screened pigs using pen floor fecal samples and chewing rope liquid.

Treatments of porcine fecal samples affect high-throughput virome sequencing results

The porcine enteric microbiota is currently extensively studied, taking advantage of developments in high-throughput sequencing technologies. However, the viral part of the microbiota, the virome, is being lightly explored, and the impact of the pretreatments used before sequencing the viruses is barely considered. In this study, the impacts of filtration, RNase and DNase treatments on virus reads recovery and diversity after sequencing on a MiSeq platform were assessed on fecal samples individually taken at <3, 5, 12 and 20 weeks from two piglets. None of the four pretreatment series affected the virus read averages or influenced diversity, but the samples with the higher proportion of reads corresponding to an entry in the "nt" database were those receiving the least number of pretreatments. The enzymatic pretreatments affected the detection of the single-stranded RNA viruses of Aichivirus C, porcine astrovirus, Sapovirus and posavirus, which is worrisome, as these viruses can be involved in swine diarrhea. If enzymatic pretreatments are used when sequencing using a high-throughput method, it may impact single-stranded RNA virus recovery, but not the overall virome diversity. Therefore, filtrated samples may be the better option, reducing the amount of bacterial genetic material while preserving the virus reads.

Metagenomics-Based Proficiency Test of Smoked Salmon Spiked with a Mock Community

An inter-laboratory proficiency test was organized to assess the ability of participants to perform shotgun metagenomic sequencing of cold smoked salmon, experimentally spiked with a mock community composed of six bacteria, one parasite, one yeast, one DNA, and two RNA viruses. Each participant applied its in-house wet-lab workflow(s) to obtain the metagenomic dataset(s), which were then collected and analyzed using MG-RAST. A total of 27 datasets were analyzed. Sample pre-processing, DNA extraction protocol, library preparation kit, and sequencing platform, influenced the abundance of specific microorganisms of the mock community. Our results highlight that despite differences in wet-lab protocols, the reads corresponding to the mock community members spiked in the cold smoked salmon, were both detected and quantified in terms of relative abundance, in the metagenomic datasets, proving the suitability of shotgun metagenomic sequencing as a genomic tool to detect microorganisms belonging to different domains in the same food matrix. The implementation of standardized wet-lab protocols would highly facilitate the comparability of shotgun metagenomic sequencing dataset across laboratories and sectors. Moreover, there is a need for clearly defining a sequencing reads threshold, to consider pathogens as detected or undetected in a food sample.

Comparison of Nucleic Acid Extraction Methods for a Viral Metagenomics Analysis of Respiratory Viruses

Viral metagenomics next-generation sequencing (mNGS) is increasingly being used to characterize the human virome. The impact of viral nucleic extraction on virome profiling has been poorly studied. Here, we aimed to compare the sensitivity and sample and reagent contamination of three extraction methods used for viral mNGS: two automated platforms (eMAG; MagNA Pure 24, MP24) and the manual QIAamp Viral RNA Mini Kit (QIAamp). Clinical respiratory samples (positive for Respiratory Syncytial Virus or Herpes Simplex Virus), one mock sample (including five viruses isolated from respiratory samples), and a no-template control (NTC) were extracted and processed through an mNGS workflow. QIAamp yielded a lower proportion of viral reads for both clinical and mock samples. The sample cross-contamination was higher when using MP24, with up to 36.09% of the viral reads mapping to mock viruses in the NTC (vs. 1.53% and 1.45% for eMAG and QIAamp, respectively). The highest number of viral reads mapping to bacteriophages in the NTC was found with QIAamp, suggesting reagent contamination. Our results highlight the importance of the extraction method choice for accurate virome characterization.

Next-generation sequencing in the diagnosis of viral encephalitis: sensitivity and clinical limitations

Identification of pathogens causing viral encephalitis remains challenging, and in over 50% of cases the etiologic factor remains undetermined. Next-generation sequencing (NGS) based metagenomics has been successfully used to detect novel and rare infections, but its value for routine diagnosis of encephalitis remains unclear. The aim of the present study was to determine the sensitivity of shotgun metagenomic sequencing protocols, which include preamplification, and testing it against cerebrospinal fluid (CSF) samples from encephalitis patients. For sensitivity testing HIV and HBV positive sera were serially diluted in CSF from an uninfected patient. NGS repeatedly detected HIV and HBV sequences present at concentrations from 10⁵ to 10² and from 10⁵ to 10 viral copies/reaction, respectively. However, when the same protocols were applied to RT-PCR/PCR positive CSF samples from 6 patients with enteroviral encephalitis (median viral load 47 copies/ml) and 15 patients with HSV, CMV or VZV encephalitis (median viral load 148 copies/ml), only 7 (28.6%) were identified as positive. In conclusions, while NGS has the advantage of being able to identify a wide range of potential pathogens it seems to be less sensitive compared to the standard amplification-based assays in the diagnosis of encephalitis, where low viral loads are common.

Metagenomic Virome Sequencing in Living Donor and Recipient Kidney Transplant Pairs Revealed JC Polyomavirus Transmission

Background

Before kidney transplantation, donors and recipients are routinely screened for viral pathogens using specific tests. Little is known about unrecognized viruses of the urinary tract that potentially result in transmission. Using an open metagenomic approach, we aimed to comprehensively assess virus transmission in living-donor kidney transplantation.

Methods

Living kidney donors and their corresponding recipients were enrolled at the time of transplantation. Follow-up study visits for recipients were scheduled 4–6 weeks and 1 year thereafter. At each visit, plasma and urine samples were collected and transplant recipients were evaluated for signs of infection or other transplant-related complications. For metagenomic analysis, samples were enriched for viruses, amplified by anchored random polymerase chain reaction (PCR), and sequenced using high-throughput metagenomic sequencing. Viruses detected by sequencing were confirmed using real-time PCR.

Results

We analyzed a total of 30 living kidney donor and recipient pairs, with a follow-up of at least 1 year. In addition to viruses commonly detected during routine post-transplant virus monitoring, metagenomic sequencing detected JC polyomavirus (JCPyV) in the urine of 7 donors and their corresponding recipients. Phylogenetic analysis confirmed infection with the donor strain in 6 cases, suggesting transmission from the transplant donor to the recipient, despite recipient seropositivity for JCPyV at the time of transplantation.

Conclusions

Metagenomic sequencing identified frequent transmission of JCPyV from kidney transplant donors to recipients. Considering the high incidence rate, future studies within larger cohorts are needed to define the relevance of JCPyV infection and the donor’s virome for transplant outcomes.

An Optimized Metagenomic Approach for Virome Detection of Clinical Pharyngeal Samples With Respiratory Infection

Respiratory virus infections are one of the major causes of acute respiratory disease or exacerbation of chronic obstructive pulmonary disease (COPD). However, next-generation sequencing has not been used for routine viral detection in clinical respiratory samples owing to its sophisticated technology. Here, several pharyngeal samples with COPD were collected to enrich viral particles using an optimized method (M3), which involved M1 with centrifugation, filtration, and concentration, M2 (magnetic beads) combined with mixed nuclease digestion, and M4 with no pretreatment as a control. Metagenomic sequencing and bioinformatics analyses showed that the M3 method for viral enrichment was superior in both viral sequencing composition and viral taxa when compared to M1, M2, and M4. M3 acquired the most viral reads and more complete sequences within 15-h performance, indicating that it might be feasible for viral detection in multiple respiratory samples in clinical practice. Based on sequence similarity analysis, 12 human viruses, including nine Anelloviruses and three coronaviruses, were characterized. Coronavirus OC43 with the largest number of viral reads accounted for nearly complete (99.8%) genome sequences, indicating that it may be a major viral pathogen involved in exacerbation of COPD.

Co-occurrence of Campylobacter Species in Children From Eastern Ethiopia, and Their Association With Environmental Enteric Dysfunction, Diarrhea, and Host Microbiome

High Campylobacter prevalence during early childhood has been associated with stunting and environmental enteric dysfunction (EED), especially in low resource settings. This study assessed the prevalence, diversity, abundance, and co-occurrence of Campylobacter spp. in stools from children in a rural area of eastern Ethiopia and their association with microbiome, diarrhea, and EED in children. Stool samples (n = 100) were collected from randomly selected children (age range: 360-498 days) in five kebeles in Haramaya District, Ethiopia. Diarrhea, compromised gut permeability, and gut inflammation were observed in 48, 45, and 57% of children, respectively. Campylobacter prevalence and species diversity were assessed using PCR and meta-total RNA sequencing (MeTRS). The prevalence of Campylobacter spp. in the children's stools was 50% (41-60%) by PCR and 88% (80-93.6%) by MeTRS (P < 0.01). Further, seven Campylobacter species (Campylobacter jejuni, Campylobacter upsaliensis, Campylobacter hyointestinalis, Campylobacter coli, Campylobacter sp. RM6137, uncultured Campylobacter sp., and Campylobacter sp. RM12175) were detected by MeTRS in at least 40% of children stools in high abundance (>1.76-log read per million per positive stool sample). Four clusters of Campylobacter species (5-12 species per cluster) co-occurred in the stool samples, suggesting that Campylobacter colonization of children may have occurred through multiple reservoirs or from a reservoir in which several Campylobacter species may co-inhabit. No associations between Campylobacter spp., EED, and diarrhea were detected in this cross-sectional study; however, characteristic microbiome profiles were identified based on the prevalence of Campylobacter spp., EED severity, and diarrhea. Forty-seven bacterial species were correlated with Campylobacter, and 13 of them also correlated with gut permeability, gut inflammation and/or EED severity. Forty-nine species not correlated with Campylobacter were correlated with gut permeability, gut inflammation, EED severity and/or diarrhea. This study demonstrated that (1) in addition to C. jejuni and C. coli, multiple non-thermophilic Campylobacter spp. (i.e., Campylobacter hyointestinalis, Campylobacter fetus, and Campylobacter concisus) were frequently detected in the children's stools and (2) the Campylobacter, gut permeability, gut inflammation, EED severity, and diarrhea were associated with characteristic microbiome composition. Additional spatial and longitudinal studies are needed to identify environmental reservoirs and sources of infection of children with disparate Campylobacter species and to better define their associations with EED in low-income countries.

Portable sequencer in the fight against infectious disease

Infectious disease is still a major threat in the world today. Five decades ago, it was considered soon to be eradicated, but the adaptation of pathogens to environmental pressure, such as antimicrobials, encouraged the emergence and reemergence of infectious disease. The fight with infectious disease starts with prevention, diagnosis, and treatment. Diagnosis can be upheld by observing the cause of disease under the microscope or detecting the presence of nucleic acid and proteins of the pathogens. The molecular techniques span from classical polymerase chain reaction (PCR) to sequencing the nucleic acid composition. Here, we are reviewing the works have been undertaken to utilize a portable sequencer, MinION, in various aspects of infectious disease management.

Application of shotgun metagenomics to smoked salmon experimentally spiked: Comparison between sequencing and microbiological data using different bioinformatic approaches

The aims of this study were i) to evaluate the possibility to detect and possibly quantify microorganisms belonging to different domains experimentally spiked in smoked salmon at known concentrations using shotgun metagenomics; ii) to compare the sequencing results using four bioinformatic tools. The salmon was spiked with six species of bacteria, including potential foodborne pathogens, as well as Cryptosporidium parvum, Saccharomyces cerevisiae and Bovine alphaherpesvirus 1. After spiking, the salmon was kept refrigerated before DNA extraction, library preparation and sequencing at 7 Gbp in paired ends at 150 bp. The bioinformatic tools named MG-RAST, OneCodex, CosmosID and MgMapper were used for the sequence analysis and the data provided were compared using STAMP. All bacteria spiked in the salmon were identified using all bioinformatic tools. Such tools were also able to assign the higher abundances to the species Propionibacterium freudenreichii spiked at the highest concentration in comparison to the other bacteria. Nevertheless, different abundances were quantified for bacteria spiked in the salmon at the same cell concentration. Cryptosporidium parvum was detected by all bioinformatics tools, while Saccharomyces cerevisiae by MG-RAST only. Finally, the DNA virus was detected by CosmosID and OneCodex only. Overall, the results of this study showed that shotgun metagenomics can be applied to detect microorganisms belonging to different domains in the same food sample. Nevertheless, a direct correlation between cell concentration of each spiked microorganism and number of corresponding reads cannot be established yet.

Evaluation of Direct Metagenomics and Target Enriched Approaches for High-throughput Sequencing of Field Rabies Viruses

INTRODUCTION

High-throughput sequencing (HTS) identifies random viral fragments in environmental samples metagenomically. High reliability gains it broad application in virus evolution, host-virus interaction, and pathogenicity studies. Deep sequencing of field samples with content of host genetic material and bacteria often produces insufficient data for metagenomics and must be preceded by target enrichment. The main goal of the study was the evaluation of HTS for complete genome sequencing of field-case rabies viruses (RABVs).

MATERIAL AND METHODS

The material was 23 RABVs isolated mainly from red foxes and one European bat lyssavirus-1 isolate propagated in neuroblastoma cells. Three methods of RNA isolation were tested for the direct metagenomics and RABV-enriched approaches. Deep sequencing was performed with a MiSeq sequencer (Illumina) and reagent v3 kit. Bioinformatics data were evaluated by Kraken and Centrifuge software and de novo assembly was done with metaSPAdes.

RESULTS

Testing RNA extraction procedures revealed the deep sequencing scope superiority of the combined TRIzol/column method. This HTS methodology made it possible to obtain complete genomes of all the RABV isolates collected in the field. Significantly greater rates of RABV genome coverages (over 5,900) were obtained with RABV enrichment. Direct metagenomic studies sequenced the full length of 6 out of 16 RABV isolates with a medium coverage between 1 and 71.

CONCLUSION

Direct metagenomics gives the most realistic illustration of the field sample microbiome, but with low coverage. For deep characterisation of viruses, e.g. for spatial and temporal phylogeography during outbreaks, target enrichment is recommended as it covers sequences much more completely.

Filovirus-reactive antibodies in humans and bats in Northeast India imply zoonotic spillover

Bats are reservoirs for several zoonotic pathogens, including filoviruses. Recent work highlights the diversity of bat borne filoviruses in Asia. High risk activities at the bat-human interface pose the threat of zoonotic virus transmission. We present evidence for prior exposure of bat harvesters and two resident fruit bat species to filovirus surface glycoproteins by screening sera in a multiplexed serological assay. Antibodies reactive to two antigenically distinct filoviruses were detected in human sera and to three individual filoviruses in bats in remote Northeast India. Sera obtained from Eonycteris spelaea bats showed similar patterns of cross-reactivity as human samples, suggesting them as the species responsible for the spillover. In contrast, sera from Rousettus leschenaultii bats reacted to two different virus glycoproteins. Our results indicate circulation of several filoviruses in bats and the possibility for filovirus transmission from bats to humans.

Focused virus surveillance at human-wildlife interfaces enables proactive detection of potentially epidemic pathogens. Filoviruses, including ebolaviruses and marburgviruses, are pathogens with epidemic potential. They were previously detected in bats and have caused disease outbreaks in humans with a high case fatality rate. Here, we tested sera obtained from bats and humans at a high-risk interface for the presence of filovirus reactive antibodies. Human participants were engaged in annual bat hunts, possibly exposing them to bat-borne viruses. We report the exposure of humans to filoviruses that were likely derived from the two sampled bat species. The bats contain antibodies raised to presumably three distinct filoviruses. Our findings suggest bats in South Asia act as a reservoir host of a diverse range of filoviruses and filovirus spillover occurs through human exposure to these bats.

Two Years of Viral Metagenomics in a Tertiary Diagnostics Unit: Evaluation of the First 105 Cases

Metagenomic next-generation sequencing (mNGS) can capture the full spectrum of viral pathogens in a specimen and has the potential to become an all-in-one solution for virus diagnostics. To date, clinical application is still in an early phase and limitations remain. Here, we evaluated the impact of viral mNGS for cases analyzed over two years in a tertiary diagnostics unit. High throughput mNGS was performed upon request by the treating clinician in cases where the etiology of infection remained unknown or the initial differential diagnosis was very broad. The results were compared to conventional routine testing regarding outcome and workload. In total, 163 specimens from 105 patients were sequenced. The main sample types were cerebrospinal fluid (34%), blood (33%) and throat swabs (10%). In the majority of the cases, viral encephalitis/meningitis or respiratory infection was suspected. In parallel, conventional virus diagnostic tests were performed (mean 18.5 individually probed targets/patients). mNGS detected viruses in 34 cases (32%). While often confirmatory, in multiple cases, the identified viruses were not included in the selected routine diagnostic tests. Two years of mNGS in a tertiary diagnostics unit demonstrated the advantages of a single, untargeted approach for comprehensive, rapid and efficient virus diagnostics, confirming the utility of mNGS in complementing current routine tests.

The use of next generation sequencing for improving food safety: Translation into practice

Next Generation Sequencing (NGS) combined with powerful bioinformatic approaches are revolutionising food microbiology. Whole genome sequencing (WGS) of single isolates allows the most detailed comparison possible hitherto of individual strains. The two principle approaches for strain discrimination, single nucleotide polymorphism (SNP) analysis and genomic multi-locus sequence typing (MLST) are showing concordant results for phylogenetic clustering and are complementary to each other. Metabarcoding and metagenomics, applied to total DNA isolated from either food materials or the production environment, allows the identification of complete microbial populations. Metagenomics identifies the entire gene content and when coupled to transcriptomics or proteomics, allows the identification of functional capacity and biochemical activity of microbial populations. The focus of this review is on the recent use and future potential of NGS in food microbiology and on current challenges. Guidance is provided for new users, such as public health departments and the food industry, on the implementation of NGS and how to critically interpret results and place them in a broader context. The review aims to promote the broader application of NGS technologies within the food industry as well as highlight knowledge gaps and novel applications of NGS with the aim of driving future research and increasing food safety outputs from its wider use.

Operational models and criteria for incorporating microbial whole genome sequencing in hospital microbiology - A systematic literature review

BACKGROUND

Microbial whole genome sequencing (WGS) has many advantages over standard microbiological methods. However, it is not yet widely implemented in routine hospital diagnostics due to notable challenges.

OBJECTIVES

The aim was to extract managerial, financial and clinical criteria supporting the decision to implement WGS in routine diagnostic microbiology, across different operational models of implementation in the hospital setting.

METHODS

This was a systematic review of literature identified through PubMed and Web of Science. English literature studies discussing the applications of microbial WGS without limitation on publication date were eligible. A narrative approach for categorization and synthesis of the sources identified was adopted.

RESULTS

A total of 98 sources were included. Four main alternative operational models for incorporating WGS in clinical microbiology laboratories were identified: full in-house sequencing and analysis, full outsourcing of sequencing and analysis and two hybrid models combining in-house/outsourcing of the sequencing and analysis components. Six main criteria (and multiple related sub-criteria) for WGS implementation emerged from our review and included cost (e.g. the availability of resources for capital and operational investment); manpower (e.g. the ability to provide training programmes or recruit trained personnel), laboratory infrastructure (e.g. the availability of supplies and consumables or sequencing platforms), bioinformatics requirements (e.g. the availability of valid analysis tools); computational infrastructure (e.g. the availability of storage space or data safety arrangements); and quality control (e.g. the existence of standardized procedures).

CONCLUSIONS

The decision to incorporate WGS in routine diagnostics involves multiple, sometimes competing, criteria and sub-criteria. Mapping these criteria systematically is an essential stage in developing policies for adoption of this technology, e.g. using a multicriteria decision tool. Future research that will prioritize criteria and sub-criteria that were identified in our review in the context of operational models will inform decision-making at clinical and managerial levels with respect to effective implementation of WGS for routine use. Beyond WGS, similar decision-making challenges are expected with respect to future integration of clinical metagenomics.

virMine: automated detection of viral sequences from complex metagenomic samples

Metagenomics has enabled sequencing of viral communities from a myriad of different environments. Viral metagenomic studies routinely uncover sequences with no recognizable homology to known coding regions or genomes. Nevertheless, complete viral genomes have been constructed directly from complex community metagenomes, often through tedious manual curation. To address this, we developed the software tool virMine to identify viral genomes from raw reads representative of viral or mixed (viral and bacterial) communities. virMine automates sequence read quality control, assembly, and annotation. Researchers can easily refine their search for a specific study system and/or feature(s) of interest. In contrast to other viral genome detection tools that often rely on the recognition of viral signature sequences, virMine is not restricted by the insufficient representation of viral diversity in public data repositories. Rather, viral genomes are identified through an iterative approach, first omitting non-viral sequences. Thus, both relatives of previously characterized viruses and novel species can be detected, including both eukaryotic viruses and bacteriophages. Here we present virMine and its analysis of synthetic communities as well as metagenomic data sets from three distinctly different environments: the gut microbiota, the urinary microbiota, and freshwater viromes. Several new viral genomes were identified and annotated, thus contributing to our understanding of viral genetic diversity in these three environments.

Quality control implementation for universal characterization of DNA and RNA viruses in clinical respiratory samples using single metagenomic next-generation sequencing workflow

Background

In recent years, metagenomic Next-Generation Sequencing (mNGS) has increasingly been used for an accurate assumption-free virological diagnosis. However, the systematic workflow evaluation on clinical respiratory samples and implementation of quality controls (QCs) is still lacking.

Methods

A total of 3 QCs were implemented and processed through the whole mNGS workflow: a no-template-control to evaluate contamination issues during the process; an internal and an external QC to check the integrity of the reagents, equipment, the presence of inhibitors, and to allow the validation of results for each sample. The workflow was then evaluated on 37 clinical respiratory samples from patients with acute respiratory infections previously tested for a broad panel of viruses using semi-quantitative real-time PCR assays (28 positive samples including 6 multiple viral infections; 9 negative samples). Selected specimens included nasopharyngeal swabs (n = 20), aspirates (n = 10), or sputums (n = 7).

Results

The optimal spiking level of the internal QC was first determined in order to be sufficiently detected without overconsumption of sequencing reads. According to QC validation criteria, mNGS results were validated for 34/37 selected samples. For valid samples, viral genotypes were accurately determined for 36/36 viruses detected with PCR (viral genome coverage ranged from 0.6 to 100%, median = 67.7%). This mNGS workflow allowed the detection of DNA and RNA viruses up to a semi-quantitative PCR Ct value of 36. The six multiple viral infections involving 2 to 4 viruses were also fully characterized. A strong correlation between results of mNGS and real-time PCR was obtained for each type of viral genome (R² ranged from 0.72 for linear single-stranded (ss) RNA viruses to 0.98 for linear ssDNA viruses).

Conclusions

Although the potential of mNGS technology is very promising, further evaluation studies are urgently needed for its routine clinical use within a reasonable timeframe. The approach described herein is crucial to bring standardization and to ensure the quality of the generated sequences in clinical setting. We provide an easy-to-use single protocol successfully evaluated for the characterization of a broad and representative panel of DNA and RNA respiratory viruses in various types of clinical samples.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3446-5) contains supplementary material, which is available to authorized users.

Metagenomic analysis of viral nucleic acid extraction methods in respiratory clinical samples

BACKGROUND

Numerous protocols for viral enrichment and genome amplification have been created. However, the direct identification of viral genomes from clinical specimens using next-generation sequencing (NGS) still has its challenges. As a selected viral nucleic acid extraction method may determine the sensitivity and reliability of NGS, it is still valuable to evaluate the extraction efficiency of different extraction kits using clinical specimens directly.

RESULTS

In this study, we performed qRT-PCR and viral metagenomic analysis of the extraction efficiency of four commonly used Qiagen extraction kits: QIAamp Viral RNA Mini Kit (VRMK), QIAamp MinElute Virus Spin Kit (MVSK), RNeasy Mini Kit (RMK), and RNeasy Plus Micro Kit (RPMK), using a mixed respiratory clinical sample without any pre-treatment. This sample contained an adenovirus (ADV), influenza virus A (Flu A), human parainfluenza virus 3 (PIV3), human coronavirus OC43 (OC43), and human metapneumovirus (HMPV). The quantity and quality of the viral extracts were significantly different among these kits. The highest threshold cycle(Ct)values for ADV and OC43 were obtained by using the RPMK. The MVSK had the lowest Ct values for ADV and PIV3. The RMK revealed the lowest detectability for HMPV and PIV3. The most effective rate of NGS data at 67.47% was observed with the RPMK. The other three kits ranged between 12.1-26.79% effectiveness rates for the NGS data. Most importantly, compared to the other three kits the highest proportion of non-host reads was obtained by the RPMK. The MVSK performed best with the lowest Ct value of 20.5 in the extraction of ADV, while the RMK revealed the best extraction efficiency by NGS analysis.

CONCLUSIONS

The evaluation of viral nucleic acid extraction efficiency is different between NGS and qRT-PCR analysis. The RPMK was most applicable for the metagenomic analysis of viral RNA and enabled more sensitive identification of the RNA virus genome in respiratory clinical samples. In addition, viral RNA extraction kits were also applicable for metagenomic analysis of the DNA virus. Our results highlighted the importance of nucleic acid extraction kit selection, which has a major impact on the yield and number of viral reads by NGS analysis. Therefore, the choice of extraction method for a given viral pathogen needs to be carefully considered.

Temporal dynamics of the lung and plasma viromes in lung transplant recipients

The human virome plays an important role for the clinical outcome of lung transplant recipients (LTRs). While pathogenic viruses may cause severe infections, non-pathogenic viruses may serve as potential markers for the level of immunosuppression. However, neither the complexity of the virome in different compartments nor the dynamics of the virus populations posttransplantation are yet understood. Therefore, in this study the virome was analyzed by metagenomic sequencing in simultaneously withdrawn bronchoalveolar lavage (BAL) and plasma samples of 15 LTRs. In seven patients, also follow-up samples were investigated for abundance and dynamics of virus populations posttransplantation. Five eukaryotic and two prokaryotic virus families were identified in BAL, and nine eukaryotic and two prokaryotic families in plasma. Anelloviruses were the most abundant in both compartments, followed by Herpes- and Coronaviruses. Virus abundance was significantly higher in LTRs than in healthy controls (Kruskal-Wallis test, p<0.001). Up to 48 different anellovirus strains were identified within a single LTR. Analyses in the follow-up patients revealed for the first time a highly complex and unique dynamics of individual anellovirus strains in the posttransplantation period. The abundance of anelloviruses in plasma was inversely correlated with that of other eukaryotic viruses (Pearson correlation coefficient r = -0.605; p<0.05). A broad spectrum of virus strains co-exists in BAL and plasma of LTRs. Especially for the anelloviruses, a high degree of co-infections and a highly individual and complex dynamics after transplantation was observed. The biological impact of these findings and their relation to clinical variables remain to be elucidated by future analyses.

Evaluation of bias induced by viral enrichment and random amplification protocols in metagenomic surveys of saliva DNA viruses

BACKGROUND

Viruses are key players regulating microbial ecosystems. Exploration of viral assemblages is now possible thanks to the development of metagenomics, the most powerful tool available for studying viral ecology and discovering new viruses. Unfortunately, several sources of bias lead to the misrepresentation of certain viruses within metagenomics workflows, hindering the shift from merely descriptive studies towards quantitative comparisons of communities. Therefore, benchmark studies on virus enrichment and random amplification protocols are required to better understand the sources of bias.

RESULTS

We assessed the bias introduced by viral enrichment on mock assemblages composed of seven DNA viruses, and the bias from random amplification methods on human saliva DNA viromes, using qPCR and deep sequencing, respectively. While iodixanol cushions and 0.45 μm filtration preserved the original composition of nuclease-protected viral genomes, low-force centrifugation and 0.22 μm filtration removed large viruses. Comparison of unamplified and randomly amplified saliva viromes revealed that multiple displacement amplification (MDA) induced stochastic bias from picograms of DNA template. However, the type of bias shifted to systematic using 1 ng, with only a marginal influence by amplification time. Systematic bias consisted of over-amplification of small circular genomes, and under-amplification of those with extreme GC content, a negative bias that was shared with the PCR-based sequence-independent, single-primer amplification (SISPA) method. MDA based on random priming provided by a DNA primase activity slightly outperformed those based on random hexamers and SISPA, which may reflect differences in ability to handle sequences with extreme GC content. SISPA viromes showed uneven coverage profiles, with high coverage peaks in regions with low linguistic sequence complexity. Despite misrepresentation of certain viruses after random amplification, ordination plots based on dissimilarities among contig profiles showed perfect overlapping of related amplified and unamplified saliva viromes and strong separation from unrelated saliva viromes. This result suggests that random amplification bias has a minor impact on beta diversity studies.

CONCLUSIONS

Benchmark analyses of mock and natural communities of viruses improve understanding and mitigate bias in metagenomics surveys. Bias induced by random amplification methods has only a minor impact on beta diversity studies of human saliva viromes.

Protocol for Generating Infectious RNA Viromes from Complex Biological Samples

This chapter proposes a simple, standardized protocol for generating RNA viromes from complex host-associated biological samples of various origins. Compared to other existing protocols to generate RNA viromes, this protocol preserves the infectivity of viral particles and allows for downstream applications such as viral characterization and isolation tests.