Multicenter benchmarking of short and long read wet lab protocols for clinical viral metagenomics.

Metagenomics is gradually being implemented for diagnosing infectious diseases. However, in-depth protocol comparisons for viral detection have been limited to individual sets of experimental workflows and laboratories. In this study, we present a benchmark of metagenomics protocols used in clinical diagnostic laboratories initiated by the European Society for Clinical Virology (ESCV) Network on NGS (ENNGS). A mock viral reference panel was designed to mimic low biomass clinical specimens. The panel was used to assess the performance of twelve metagenomic wet lab protocols currently in use in the diagnostic laboratories of participating ENNGS member institutions. Both Illumina and Nanopore, shotgun and targeted capture probe protocols were included. Performance metrics sensitivity, specificity, and quantitative potential were assessed using a central bioinformatics pipeline. Overall, viral pathogens with loads down to 104 copies/ml (corresponding to CT values of 31 in our PCR assays) were detected by all the evaluated metagenomic wet lab protocols. In contrast, lower abundant mixed viruses of CT values of 35 and higher were detected only by a minority of the protocols. Considering the reference panel as the gold standard, optimal thresholds to define a positive result were determined per protocol, based on the horizontal genome coverage. Implementing these thresholds, sensitivity and specificity of the protocols ranged from 67 to 100 % and 87 to 100 %, respectively. A variety of metagenomic protocols are currently in use in clinical diagnostic laboratories. Detection of low abundant viral pathogens and mixed infections remains a challenge, implying the need for standardization of metagenomic analysis for use in clinical settings.

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.

Characterization of SARS-CoV-2 replication in human H1299/ACE2 cells: A versatile and practical infection model for antiviral research and beyond.

A range of cell culture infection models have been used to study SARS-CoV-2 and perform antiviral drug research. Commonly used African green monkey Vero, human lung-derived Calu-3 and ACE2+TMPRSS2-expressing A549 cells, each have their limitations. Here, we describe human ACE2-expressing H1299 lung cells as a more efficient and robust model for SARS-CoV-2 research. These cells are as easy to handle as Vero cells, support SARS-CoV-2 replication to high titers, display a functional innate immune response and are suitable for plaque assays, microscopy, the production of (genetically stable) virus stocks and antiviral assays. H1299/ACE2-based (CPE reduction) assays can be performed without adding a P-gP drug efflux pump inhibitor, which is often required in Vero-based assays. Moreover, H1299/ACE2 cells allowed us to perform CPE reduction assays with omicron variants that did not work in Vero-based assays. In summary, H1299/ACE2 cells are a versatile infection model to study SARS-CoV-2 replication in the context of antiviral drug development and virus-host interaction studies.

A comparison of five Illumina, Ion Torrent, and nanopore sequencing technology-based approaches for whole genome sequencing of SARS-CoV-2.

Rapid identification of the rise and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern remains critical for monitoring of the efficacy of diagnostics, therapeutics, vaccines, and control strategies. A wide range of SARS-CoV-2 next-generation sequencing (NGS) methods have been developed over the last years, but cross-sequence technology benchmarking studies have been scarce. In the current study, 26 clinical samples were sequenced using five protocols: AmpliSeq SARS-CoV-2 (Illumina), EasySeq RC-PCR SARS-CoV-2 (Illumina/NimaGen), Ion AmpliSeq SARS-CoV-2 (Thermo Fisher), custom primer sets (Oxford Nanopore Technologies (ONT)), and capture probe-based viral metagenomics (Roche/Illumina). Studied parameters included genome coverage, depth of coverage, amplicon distribution, and variant calling. The median SARS-CoV-2 genome coverage of samples with cycle threshold (Ct) values of 30 and lower ranged from 81.6 to 99.8% for, respectively, the ONT protocol and Illumina AmpliSeq protocol. Correlation of coverage with PCR Ct values varied per protocol. Amplicon distribution signatures differed across the methods, with peak differences of up to 4 log10 at disbalanced positions in samples with high viral loads (Ct values ≤ 23). Phylogenetic analyses of consensus sequences showed clustering independent of the workflow used. The proportion of SARS-CoV-2 reads in relation to background sequences, as a (cost-)efficiency metric, was the highest for the EasySeq protocol. The hands-on time was the lowest when using EasySeq and ONT protocols, with the latter additionally having the shortest sequence runtime. In conclusion, the studied protocols differed on a variety of the studied metrics. This study provides data that assist laboratories when selecting protocols for their specific setting.

Non-Toxigenic Clostridioides difficile Strain E4 (NTCD-E4) Prevents Establishment of Primary C. difficile Infection by Epidemic PCR Ribotype 027 in an In Vitro Human Gut Model.

Clostridioides difficile infection (CDI) remains a significant healthcare burden. Non-toxigenic C. difficile (NTCD) strains have shown a benefit in preventing porcine enteritis and in human recurrent CDI. In this study, we evaluated the efficacy of metronidazole-resistant NTCD-E4 in preventing CDI facilitated by a range of antimicrobials in an in vitro human gut model. NTCD-E4 spores (at a dose of 107) were instilled 7 days before a clinical ribotype (RT) 027 (at the same dose) strain (210). In separate experiments, four different antimicrobials were used to perturb gut microbiotas; bacterial populations and cytotoxin production were determined using viable counting and Vero cell cytotoxicity, respectively. RT027 and NTCD-E4 proliferated in the in vitro model when inoculated singly, with RT027 demonstrating high-level cytotoxin (3-5-log10-relative units) production. In experiments where the gut model was pre-inoculated with NTCD-E4, RT027 was remained quiescent and failed to produce cytotoxins. NTCD-E4 showed mutations in hsmA and a gene homologous to CD196-1331, previously linked to medium-dependent metronidazole resistance, but lacked other metronidazole resistance determinants. This study showed that RT027 was unable to elicit simulated infection in the presence of NTCD-E4 following stimulation by four different antimicrobials. These data complement animal and clinical studies in suggesting NTCD offer prophylactic potential in the management of human CDI.

Engineering potent live attenuated coronavirus vaccines by targeted inactivation of the immune evasive viral deubiquitinase.

Coronaviruses express a papain-like protease (PLpro) that is required for replicase polyprotein maturation and also serves as a deubiquitinating enzyme (DUB). In this study, using a Middle East respiratory syndrome virus (MERS-CoV) PLpro modified virus in which the DUB is selectively inactivated, we show that the PLpro DUB is an important MERS-CoV interferon antagonist and virulence factor. Although the DUB-negative rMERS-CoVMA replicates robustly in the lungs of human dipeptidyl peptidase 4 knock-in (hDPP4 KI) mice, it does not cause clinical symptoms. Interestingly, a single intranasal vaccination with DUB-negative rMERS-CoVMA induces strong and sustained neutralizing antibody responses and sterilizing immunity after a lethal wt virus challenge. The survival of naïve animals also significantly increases when sera from animals vaccinated with the DUB-negative rMERS-CoVMA are passively transferred, prior to receiving a lethal virus dose. These data demonstrate that DUB-negative coronaviruses could be the basis of effective modified live attenuated vaccines.

Immune Determinants of Viral Clearance in Hospitalised COVID-19 Patients: Reduced Circulating Naïve CD4+ T Cell Counts Correspond with Delayed Viral Clearance.

Virus-specific cellular and humoral responses are major determinants for protection from critical illness after SARS-CoV-2 infection. However, the magnitude of the contribution of each of the components to viral clearance remains unclear. Here, we studied the timing of viral clearance in relation to 122 immune parameters in 102 hospitalised patients with moderate and severe COVID-19 in a longitudinal design. Delayed viral clearance was associated with more severe disease and was associated with higher levels of SARS-CoV-2-specific (neutralising) antibodies over time, increased numbers of neutrophils, monocytes, basophils, and a range of pro-inflammatory cyto-/chemokines illustrating ongoing, partially Th2 dominating, immune activation. In contrast, early viral clearance and less critical illness correlated with the peak of neutralising antibodies, higher levels of CD4 T cells, and in particular naïve CD4+ T cells, suggesting their role in early control of SARS-CoV-2 possibly by proving appropriate B cell help. Higher counts of naïve CD4+ T cells also correlated with lower levels of MIF, IL-9, and TNF-beta, suggesting an indirect role in averting prolonged virus-induced tissue damage. Collectively, our data show that naïve CD4+ T cell play a critical role in rapid viral T cell control, obviating aberrant antibody and cytokine profiles and disease deterioration. These data may help in guiding risk stratification for severe COVID-19.

Genome-wide diversity of Zika virus: Exploring spatio-temporal dynamics to guide a new nomenclature proposal.

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized 'African' and 'Asian' genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

Performance of Five Metagenomic Classifiers for Virus Pathogen Detection Using Respiratory Samples from a Clinical Cohort.

Viral metagenomics is increasingly applied in clinical diagnostic settings for detection of pathogenic viruses. While several benchmarking studies have been published on the use of metagenomic classifiers for abundance and diversity profiling of bacterial populations, studies on the comparative performance of the classifiers for virus pathogen detection are scarce. In this study, metagenomic data sets (n = 88) from a clinical cohort of patients with respiratory complaints were used for comparison of the performance of five taxonomic classifiers: Centrifuge, Clark, Kaiju, Kraken2, and Genome Detective. A total of 1144 positive and negative PCR results for a total of 13 respiratory viruses were used as gold standard. Sensitivity and specificity of these classifiers ranged from 83 to 100% and 90 to 99%, respectively, and was dependent on the classification level and data pre-processing. Exclusion of human reads generally resulted in increased specificity. Normalization of read counts for genome length resulted in a minor effect on overall performance, however it negatively affected the detection of targets with read counts around detection level. Correlation of sequence read counts with PCR Ct-values varied per classifier, data pre-processing (R2 range 15.1-63.4%), and per virus, with outliers up to 3 log10 reads magnitude beyond the predicted read count for viruses with high sequence diversity. In this benchmarking study, sensitivity and specificity were within the ranges of use for diagnostic practice when the cut-off for defining a positive result was considered per classifier.

Longitudinal Monitoring of DNA Viral Loads in Transplant Patients Using Quantitative Metagenomic Next-Generation Sequencing.

INTRODUCTION: Immunocompromised patients are prone to reactivations and (re-)infections of multiple DNA viruses. Viral load monitoring by single-target quantitative PCRs (qPCR) is the current cornerstone for virus quantification. In this study, a metagenomic next-generation sequencing (mNGS) approach was used for the identification and load monitoring of transplantation-related DNA viruses. METHODS: Longitudinal plasma samples from six patients that were qPCR-positive for cytomegalovirus (CMV), Epstein-Barr virus (EBV), BK polyomavirus (BKV), adenovirus (ADV), parvovirus B19 (B19V), and torque teno-virus (TTV) were sequenced using the quantitative metagenomic Galileo Viral Panel Solution (Arc Bio, LLC, Cambridge, MA, USA) reagents and bioinformatics pipeline combination. Qualitative and quantitative performance was analysed with a focus on viral load ranges relevant for clinical decision making. RESULTS: All pathogens identified by qPCR were also identified by mNGS. BKV, CMV, and HHV6B were additionally detected by mNGS, and could be confirmed by qPCR or auxiliary bioinformatic analysis. Viral loads determined by mNGS correlated with the qPCR results, with inter-method differences in viral load per virus ranging from 0.19 log10 IU/mL for EBV to 0.90 log10 copies/mL for ADV. TTV, analysed by mNGS in a semi-quantitative way, demonstrated a mean difference of 3.0 log10 copies/mL. Trends over time in viral load determined by mNGS and qPCR were comparable, and clinical thresholds for initiation of treatment were equally identified by mNGS. CONCLUSIONS: The Galileo Viral Panel for quantitative mNGS performed comparably to qPCR concerning detection and viral load determination, within clinically relevant ranges of patient management algorithms.

Translating genomic exploration of the family Polyomaviridae into confident human polyomavirus detection.

The Polyomaviridae is a family of ubiquitous dsDNA viruses that establish persistent infection early in life. Screening for human polyomaviruses (HPyVs), which comprise 14 diverse species, relies upon species-specific qPCRs whose validity may be challenged by accelerating genomic exploration of the virosphere. Using this reasoning, we tested 64 published HPyV qPCR assays in silico against the 1781 PyV genome sequences that were divided in targets and nontargets, based on anticipated species specificity of each qPCR. We identified several cases of problematic qPCR performance that were confirmed in vitro and corrected through using degenerate oligos. Furthermore, our study ranked 8 out of 52 tested BKPyV qPCRs as remaining of consistently high quality in the wake of recent PyV discoveries and showed how sensitivity of most other qPCRs could be rescued by annealing temperature adjustment. This study establishes an efficient framework for ensuring confidence in available HPyV qPCRs in the genomic era.

Viral metagenomic sequencing in a cohort of international travellers returning with febrile illness.

BACKGROUND: Diagnosis of infections in returning international travellers can be challenging because of the broad spectrum of potential infectious etiologies potentially involved. Viral metagenomic next-generation sequencing (mNGS) has the potential to detect any virus present in a patient sample and is increasingly being used for difficult to diagnose cases. The aim of this study was to analyze the performance of mNGS for viral pathogen detection in the clinical setting of international travellers returning with febrile illness. METHODS: Thirty-eight serum samples from international travellers returning with febrile illness and presenting at the outpatient clinic of the Leiden University Medical Center in the Netherlands in the time period 2015-2016 were selected retrospectively. Samples were processed for viral metagenomic sequencing using a probe panel capturing all known vertebrate viruses. Bioinformatic analysis was performed using Genome Detective software for metagenomic virus detection. Metagenomic virus findings were compared with viral pathogen detection using conventional methods. RESULTS: In 8 out of the 38 patients (21%), a pathogenic virus was detected by mNGS. All viral pathogens detected by conventional assays were also detected by mNGS: dengue virus (n=4 patients), Epstein-Barr virus (n=2), hepatitis B virus (n=1). In addition, mNGS resulted in additional pathogenic findings in 2 patients (5%): dengue virus (n=1), and hepatitis C virus (n=1). Non-pathogenic viruses detected were: GB virus C (n=1) and torque teno viruses (n=3). High genome coverage and depth using capture probes enabled typing of the dengue viruses detected. CONCLUSIONS: Viral metagenomics has the potential to assist the detection of viral pathogens and co-infections in one step in international travellers with a febrile syndrome. Furthermore, viral enrichment by probes resulted in high genome coverage and depth which enabled dengue virus typing.

Measurement and 3D Visualization of the Human Internal Heat Field by Means of Microwave Radiometry.

The possibility of non-invasive determination of the depth of the location and temperature of a cancer tumor in the human body by multi-frequency three-dimensional (3D) radiothermography is considered. The models describing the receiving of the human body's own radiothermal field processes are presented. The analysis of the possibility of calculating the desired parameters based on the results of measuring antenna temperatures simultaneously in two different frequency ranges is performed. Methods of displaying on the monitor screen the three-dimensional temperature distribution of the subcutaneous layer of the human body, obtained as a result of data processing of a multi-frequency multichannel radiothermograph, are considered. The possibility of more accurate localization of hyperthermia focus caused by the presence of malignant tumors in the depth of the human body with multi-frequency volumetric radiothermography is shown. The results of the study of various methods of data interpolation for displaying the continuous intrinsic radiothermal field of the human body are presented. Examples of displaying the volumetric temperature distribution by the moving plane method based on digital models and the results of an experimental study of the thermal field of the human body and head are given.

Benchmark of thirteen bioinformatic pipelines for metagenomic virus diagnostics using datasets from clinical samples.

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

Torque teno virus loads after kidney transplantation predict allograft rejection but not viral infection.

The main challenge of immunosuppressive therapy after solid organ transplantation is to create a new immunological balance that prevents organ rejection and does not promote opportunistic infection. Torque teno virus (TTV), a ubiquitous and non-pathogenic single-stranded DNA virus, has been proposed as a marker of functional immunity in immunocompromised patients. Here we investigate whether TTV loads predict the risk of common viral infection and allograft rejection in kidney transplantation recipients. In a retrospective cohort of 389 kidney transplantation recipients, individual TTV loads in were measured by qPCR in consecutive plasma samples during one year follow-up. The endpoints were allograft rejection, BK polyomavirus (BKPyV) viremia and cytomegalovirus (CMV) viremia. Repeated TTV measurements and rejection and infection survival data were analysed in a joint model. During follow-up, TTV DNA detection in the transplant recipients increased from 85 to 100%. The median viral load increased to 107 genome copies/ml within three months after transplantation. Rejection, BKPyV viremia and CMV viremia occurred in 23%, 27% and 17% of the patients, respectively. With every 10-fold TTV load-increase, the risk of rejection decreased considerably (HR: 0.74, CI 95%: 0.71-0.76), while the risk of BKPyV and CMV viremia remained the same (HR: 1.03, CI 95%: 1.03-1.04 and HR: 1.01, CI 95%: 1.01-1.01). In conclusion, TTV load kinetics predict allograft rejection in kidney transplantation recipients, but not the BKPyV and CMV infection. The potential use of TTV load levels as a guide for optimal immunosuppressive drug dosage to prevent allograft rejection deserves further validation.

Haem is crucial for medium-dependent metronidazole resistance in clinical isolates of Clostridioides difficile.

BACKGROUND: Until recently, metronidazole was the first-line treatment for Clostridioides difficile infection and it is still commonly used. Though resistance has been reported due to the plasmid pCD-METRO, this does not explain all cases. OBJECTIVES: To identify factors that contribute to plasmid-independent metronidazole resistance of C. difficile. METHODS: Here, we investigate resistance to metronidazole in a collection of clinical isolates of C. difficile using a combination of antimicrobial susceptibility testing on different solid agar media and WGS of selected isolates. RESULTS: We find that nearly all isolates demonstrate a haem-dependent increase in the MIC of metronidazole, which in some cases leads to isolates qualifying as resistant (MIC >2 mg/L). Moreover, we find an SNP in the haem-responsive gene hsmA, which defines a metronidazole-resistant lineage of PCR ribotype 010/MLST ST15 isolates that also includes pCD-METRO-containing strains. CONCLUSIONS: Our data demonstrate that haem is crucial for medium-dependent metronidazole resistance in C. difficile.

Recommendations for the introduction of metagenomic next-generation sequencing in clinical virology, part II: bioinformatic analysis and reporting.

Metagenomic next-generation sequencing (mNGS) is an untargeted technique for determination of microbial DNA/RNA sequences in a variety of sample types from patients with infectious syndromes. mNGS is still in its early stages of broader translation into clinical applications. To further support the development, implementation, optimization and standardization of mNGS procedures for virus diagnostics, the European Society for Clinical Virology (ESCV) Network on Next-Generation Sequencing (ENNGS) has been established. The aim of ENNGS is to bring together professionals involved in mNGS for viral diagnostics to share methodologies and experiences, and to develop application guidelines. Following the ENNGS publication Recommendations for the introduction of mNGS in clinical virology, part I: wet lab procedure in this journal, the current manuscript aims to provide practical recommendations for the bioinformatic analysis of mNGS data and reporting of results to clinicians.

Rhinovirus Detection in the Nasopharynx of Children Undergoing Cardiac Surgery Is Not Associated With Longer PICU Length of Stay: Results of the Impact of Rhinovirus Infection After Cardiac Surgery in Kids (RISK) Study.

OBJECTIVES: To determine whether children with asymptomatic carriage of rhinovirus in the nasopharynx before elective cardiac surgery have an increased risk of prolonged PICU length of stay. STUDY DESIGN: Prospective, single-center, blinded observational cohort study. SETTING: PICU in a tertiary hospital in The Netherlands. PATIENTS: Children under 12 years old undergoing elective cardiac surgery were enrolled in the study after informed consent of the parents/guardians. INTERVENTIONS: The parents/guardians filled out a questionnaire regarding respiratory symptoms. On the day of the operation, a nasopharyngeal swab was obtained. Clinical data were collected during PICU admission, and PICU/hospital length of stay were reported. If a patient was still intubated 3 days after operation, an additional nasopharyngeal swab was collected. Nasopharyngeal swabs were tested for rhinovirus and other respiratory viruses with polymerase chain reaction. MEASUREMENTS AND MAIN RESULTS: Of the 163 included children, 74 (45%) tested rhinovirus positive. Rhinovirus-positive patients did not have a prolonged PICU length of stay (median 2 d each; p = 0.257). Rhinovirus-positive patients had a significantly shorter median hospital length of stay compared with rhinovirus-negative patients (8 vs 9 d, respectively; p = 0.006). Overall, 97 of the patients (60%) tested positive for one or more respiratory virus. Virus-positive patients had significantly shorter PICU and hospital length of stay, ventilatory support, and nonmechanical ventilation. Virus-negative patients had respiratory symptoms suspected for a respiratory infection more often. In 31% of the children, the parents reported mild upper respiratory complaints a day prior to the cardiac surgery, this was associated with postextubation stridor, but no other clinical outcome measures. CONCLUSIONS: Preoperative rhinovirus polymerase chain reaction positivity is not associated with prolonged PICU length of stay. Our findings do not support the use of routine polymerase chain reaction testing for respiratory viruses in asymptomatic children admitted for elective cardiac surgery.

Coronavirus discovery by metagenomic sequencing: a tool for pandemic preparedness.

INTRODUCTION: The SARS-CoV-2 pandemic of 2020 is a prime example of the omnipresent threat of emerging viruses that can infect humans. A protocol for the identification of novel coronaviruses by viral metagenomic sequencing in diagnostic laboratories may contribute to pandemic preparedness. AIM: The aim of this study is to validate a metagenomic virus discovery protocol as a tool for coronavirus pandemic preparedness. METHODS: The performance of a viral metagenomic protocol in a clinical setting for the identification of novel coronaviruses was tested using clinical samples containing SARS-CoV-2, SARS-CoV, and MERS-CoV, in combination with databases generated to contain only viruses of before the discovery dates of these coronaviruses, to mimic virus discovery. RESULTS: Classification of NGS reads using Centrifuge and Genome Detective resulted in assignment of the reads to the closest relatives of the emerging coronaviruses. Low nucleotide and amino acid identity (81% and 84%, respectively, for SARS-CoV-2) in combination with up to 98% genome coverage were indicative for a related, novel coronavirus. Capture probes targeting vertebrate viruses, designed in 2015, enhanced both sequencing depth and coverage of the SARS-CoV-2 genome, the latter increasing from 71% to 98%. CONCLUSION: The model used for simulation of virus discovery enabled validation of the metagenomic sequencing protocol. The metagenomic protocol with virus probes designed before the pandemic, can assist the detection and identification of novel coronaviruses directly in clinical samples.

STROBE-metagenomics: a STROBE extension statement to guide the reporting of metagenomics studies.

The term metagenomics refers to the use of sequencing methods to simultaneously identify genomic material from all organisms present in a sample, with the advantage of greater taxonomic resolution than culture or other methods. Applications include pathogen detection and discovery, species characterisation, antimicrobial resistance detection, virulence profiling, and study of the microbiome and microecological factors affecting health. However, metagenomics involves complex and multistep processes and there are important technical and methodological challenges that require careful consideration to support valid inference. We co-ordinated a multidisciplinary, international expert group to establish reporting guidelines that address specimen processing, nucleic acid extraction, sequencing platforms, bioinformatics considerations, quality assurance, limits of detection, power and sample size, confirmatory testing, causality criteria, cost, and ethical issues. The guidance recognises that metagenomics research requires pragmatism and caution in interpretation, and that this field is rapidly evolving.