key: cord-0747705-jk0p61fy
authors: Brinkmann, Annika; Uddin, Steven; Ulm, Sophie-Luisa; Pape, Katharina; Förster, Sophie; Enan, Khalid; Nourlil, Jalal; Krause, Eva; Schaade, Lars; Michel, Janine; Nitsche, Andreas
title: RespiCoV: Simultaneous identification of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and 46 respiratory tract viruses and bacteria by amplicon-based Oxford-Nanopore MinION sequencing
date: 2022-03-09
journal: PLoS One
DOI: 10.1371/journal.pone.0264855
sha: 597ad1a8fafc414a25036c623684f48da52a6c0d
doc_id: 747705
cord_uid: jk0p61fy

Since December 2019 the world has been facing the outbreak of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Identification of infected patients and discrimination from other respiratory infections have so far been accomplished by using highly specific real-time PCRs. Here we present a rapid multiplex approach (RespiCoV), combining highly multiplexed PCRs and MinION sequencing suitable for the simultaneous screening for 41 viral and five bacterial agents related to respiratory tract infections, including the human coronaviruses NL63, HKU1, OC43, 229E, Middle East respiratory syndrome coronavirus, SARS-CoV, and SARS-CoV-2. RespiCoV was applied to 150 patient samples with suspected SARS-CoV-2 infection and compared with specific real-time PCR. Additionally, several respiratory tract pathogens were identified in samples tested positive or negative for SARS-CoV-2. Finally, RespiCoV was experimentally compared to the commercial RespiFinder 2SMART multiplex screening assay (PathoFinder, The Netherlands).

Infections of the respiratory tract range from the mild, self-limiting common cold to lifethreatening illnesses and epidemics caused by influenza viruses, severe acute respiratory syndrome coronavirus (SARS-CoV), or Middle East respiratory syndrome coronavirus (MERS) [1, 2] . Recently, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been causing an ongoing global pandemic with more than 281,808,270 diagnosed cases of 29 December 2021 (https://covid19.who.int/). Shortly after the identification and whole-genome sequencing of the novel emerging virus, specific real-time PCRs for SARS-CoV-2 diagnostics have been developed and deployed extensively [3] [4] [5] [6] . Although up-to-date, fast, reliable, and specific real-time PCR-based SARS-CoV-2 diagnostics has the highest priority for control and containment of the Covid-19 pandemic, the identification and possible relevance of viral or bacterial co-infections for the severity of the course of Covid-19 have been addressed in many studies [7] [8] [9] [10] . During the SARS-CoV pandemic in 2003, reports on dual infections were scarce [11, 12] . However, it has been shown by a systematic review that 19% of patients with COVID-19 have bacterial and viral co-infections which are associated with poorer outcomes [13] . As patients with symptoms described for Covid-19 are usually exclusively tested for SARS-CoV-2, many patients with negative results remain undiagnosed for co-infections, which can lead to non-specific treatment or incorrect treatment of hospitalized patients as well as uncertainty regarding the patients' health status and a needless placing in quarantine.

In contrast to specific real-time PCRs, Illumina and Nanopore shotgun sequencing enable the unbiased detection of one or several pathogens simultaneously from sputum or swab samples and have previously been performed for identification of respiratory tract pathogens, including Streptococcus pneumoniae and influenza virus [14] [15] [16] . However, shotgun sequencing generates a high amount of data accompanied by high costs, involuntary sequencing of the hosts' DNA which conflicts with personal data protection, and low sensitivity of virus identification.

Here we present an amplicon-based MinION sequencing approach (referred to as Respi-CoV) with 114 primers for simultaneous diagnostics of SARS-CoV-2 and further 40 viral and five bacterial agents related to respiratory tract infections (Table 1) . This approach can contribute to the detection of co-infections in patients infected with SARS-CoV-2 and aid in differential diagnostics of patients tested negative for SARS-CoV-2.

The targeted common upper respiratory tract viruses and bacteria for RespiCoV were chosen based on the publications by Hodinka et al. [17] and Jain et al. [18] . Additionally, herpes simplex virus type 1 and Epstein-Barr virus have been associated with upper respiratory tract infections in critically ill patients, and viruses were included as targets of the RespiCoV assay [19, 20] . Varicella zoster virus and herpes simplex virus type 2 were included in the RespiCoV Panel for validation of the method and for amplification as sequencing controls. Because generic markers like 16S rRNA can identify bacteria, we included only the most prominent bacteria.

The 114 primers (S1 Table) for the RespiCoV Panel were designed by using Primer3 v2.3.7 (https://primer3.org/). For the primer design, all references for each virus and bacteria species in Table 1 were aligned by using MAFFT v7.450, and primers were placed on conserved regions. For SARS-CoV-2, primers were designed to target three genome regions of 374, 263, and 363 bp in length, respectively (bp regions 15,814-16,188 (ORF1ab); 16,821-17,084 (ORF1ab); 29,337-29,700 (N/ORF10); based on reference MT345888). These primers can detect also SARS-CoV, with significant differences in the amplified sequence enabling a clear discrimination between SARS-CoV and SARS-CoV-2 (90.6%, 91.6%, and 88.4% referred to BetaCoV/Germany/BavPat1/2020 [GISAID] und NC_004718.3). All primers have been selected based on their melting temperature Tm (minimum Tm 58˚C, maximum Tm 62˚C, optimal Tm 60˚C), resulting amplicon length (250-600 bp), and a minimum of possible hetero-dimer formation. 

The RespiCoV method was further tested with samples provided by INSTAND e.V. for a national quality assurance exercise. INSTAND e.V. offers exercises for quality assurance for medical laboratories across Germany. This exercise focused on the genomic detection of SARS-CoV-2 and contained inactivated samples for sensitivity (SARS-CoV-2 in different concentrations) and specificity (other coronaviruses). For comparison, the samples were identified with specific real-time PCRs for SARS-CoV, SARS-CoV-2, MERS-CoV, OC43, NL63, 229E, and HKU1.

The patient samples were amplified in a single reaction with the following PCR conditions: 3 μl of viral cDNA, 1.6 μl of primer pool, 0.2 mM dNTP (Invitrogen, Karlsruhe, Germany), 4 μl of 10 x Platinum Taq buffer, 2 mM MgCl 2 , and 5 U Platinum Taq polymerase (Invitrogen) with added water to a final volume of 25 μl. Cycling conditions were 94˚C for 5 min, 45 amplification cycles at 94˚C for 20 s, 65˚C for 30 s, 72˚C for 20 s, and a final extension step for 5 min (at 72˚C). Thermal cycling was performed in an Eppendorf Mastercycler Pro (Eppendorf Vertrieb Deutschland, Wesseling-Berzdorf, Germany) with a total runtime of 64 min.

Amplified samples were processed for nanopore sequencing on the MinION (Oxford Nanopore Technologies, Oxford, United Kingdom). The libraries were prepared by using the ligation sequencing kit 1D, SQK-LSK109 (Oxford Nanopore Technologies). For combined sequencing of several samples on one flow cell, samples were barcoded with the Native Barcoding Expansion Kit (EXP-NBD104 and EXP-NBD114). Subsequently, the libraries were loaded onto Oxford Nanopore MinION SpotON Flow Cells Mk I, R9.4.1. (Oxford Nanopore Technologies). Samples were run for at least 30 min.

The Fast5 data generated during sequencing was transcribed to FastQ sequences by using Guppy v.3.4.5 (Oxford Nanopore Technologies) on the MinION IT device (MNT-001).

Computational separation of the barcoded samples was performed with Guppy v.3.4.5 for Windows. FastQ files for each sample were aligned to the reference sequences with Guppy v.4.0.11 for Linux and the resulting alignments were used for read counts. Primer sequences were soft clipped with bamclipper v.1.1.1 and all soft clippings from the BAM file were removed with custom python scripts. For species identification, consensus sequences generated from the reference alignments (Geneious prime v2020.2.3) were validated using online blast. As read counts can differ between runs, samples were only rated positive when the following parameters were met: number of total reads for each sample > 0.5% of the total reads from the run, number of reads for SARS-CoV-2 > 0.5% of all total reads of SARS-CoV-2 from the run plus the reads of SARS-CoV-2 identified in the negative control, and number of reads for SARS-CoV-2 > 50.

The studies involving human participants were reviewed and approved by the Ä rztekammer Berlin (Berlin Medical Association; #Eth 20/40). The patients/participants provided their written informed consent to participate in this study.

In one of the samples tested negative with the RespiFinder 2SMART, herpes simplex virus type 1 could be identified with the RespiCoV Panel (67,321 specific amplicons), which is not targeted by the RespiFinder 2SMART. Furthermore, in three of the patient samples tested positive with the RespiFinder 2SMART, additional pathogens could be identified with the RespiCoV Panel. Streptococcus pneumoniae, which is not targeted by the RespiFinder 2SMART, could be identified additionally in two of the samples, and Rhinovirus A could be identified in one of the samples. For three of the samples identified as positive with both methods, additional species/strain information could be gained by the sequence information obtained with the Respi-CoV Panel. For example, human adenovirus could be specified further to human adenovirus type B and the lineage of influenzavirus B could be identified as Yamagata. For two of the samples tested positive with both methods, read numbers after MinION sequencing were very low (55 reads for Human respiratory syncytial virus B and 676 reads for Human metapneumovirus). For the remaining samples, 10,937-192,431 reads were sequenced in one hour, providing sufficient viral reads for identification within the first minutes of sequencing ( Table 2 ).

Of the 150 clinical samples, 66 samples were identified as negative and 84 samples were identified as positive for SARS-CoV-2 with a specific SARS-CoV-2 real-time PCR in our routine diagnostics (Cq range of 18-38, Table 3 ). With RespiCoV, 65 of the 66 negative samples were correctly identified as negative for SARS-COV-2, whereas one sample was identified as positive for SARS-CoV-2 with low read numbers of SARS-CoV-2 amplicons after sequencing (n = 4000, mean read numbers for samples within Cq range 18-28: 35,000; and 19,000 within Cq range 29-33). However, the patient had been tested negative by specific real-time PCR previously, but after a series of positive tests.

Of As shown in Fig 1, there is a good correlation between virus genome load represented by the Cq value and the read number within one sequencing run, but not between different runs (shown for three different runs).

In the 150 samples, sequences of 32 pathogens other than SARS-CoV-2 could be identified, with 23 

The RespiCoV method was further tested with samples provided by INSTAND e.V. for quality assurance of SARS-CoV-2 diagnostics in medical laboratories across Germany. The results obtained with RespiCoV were identical when compared with the real-time PCR results (Table 4 ). After 30 min of sequencing, 137,295 reads were obtained from the samples with high SARS-CoV-2 concentration (Cq 21.4). For the samples with a low concentration of SARS-CoV-2, only 4,657 reads were sequenced, but the read number was sufficient for identification of the virus within the first minutes of sequencing. 

In this study, we introduce an amplicon-based MinION sequencing approach, referred to as RespiCoV, which is able to identify and differentiate 41 viral and five bacterial species related to respiratory tract infections, including the human coronaviruses 229E, HKU1, NL63, OC43, MERS, SARS-CoV, and SARS-CoV-2, the latter challenging the world in an ongoing pandemic since 2019. We could show that the RespiCoV Panel is able to identify several viral and bacterial species in patients with symptoms of respiratory tract infections. Furthermore, samples from patients with diagnosed infections with SARS-CoV-2 were identified with the Respi-CoV Panel, even if viral load was low (up to a Cq value of 33). Although the identification was not shown experimentally for all viral and bacterial targets of the RespiCoV Panel, the performance of the method was shown for several pathogens, including influenza A virus, influenza B virus, human coronavirus OC43 and 229E, human adenovirus B, human bocavirus, human metapneumovirus, human respiratory syncytial virus, human parainfluenza virus types 2, herpes simplex virus type 1, S. pneumoniae, and SARS-CoV-2. Compared with the extensively used and validated RespiPanel 2SMART, we could show that the RespiCoV Panel can be used as an approach for the simultaneous identification of respiratory tract pathogens. Just in one case, only low read numbers of Human respiratory syncytial virus could be identified with the RespiCoV Panel, which may be the result of low virus concentration or the primer design, that could be adapted by integrating additional primers into the RespiCoV primer pool.

Although reliable, fast, and accurate real-time PCR is the gold standard for SARS-CoV-2 detection, the method described here can further contribute to the diagnostics and differential diagnostics of patients with symptoms described for Covid-19. Identification of viral and bacterial co-infections has been performed in several studies with real-time PCR, but the abundance and potential impact of these infections remained unknown. In the 2009 H1N1 influenza outbreak, co-infections of patients with H1N1 and a second respiratory virus were associated with an increased risk of complications [21] . Furthermore, in children co-infections with respiratory syncytial virus and metapneumovirus or rhinovirus were associated with a 10-fold greater risk of Pediatric Intensive Care Unit level of care [22, 23] . In contrast, other studies have found less severe clinical outcomes with viral co-infection or showed no correlation of co-infections and severity of disease [24, 25] . Co-infections of patients diagnosed for SARS-CoV-2 identified by specific real-time PCRs performed in two independent studies also included common respiratory viruses (influenza A virus, rhinovirus, human respiratory syncytial virus, human coronavirus HKU1, human parainfluenzavirus type 1, and human metapneumovirus), but infection rates were low (5.8% and 3.2%, respectively) [7, 9, 26] . Another study reported 22.4% of all patients assigned to the emergency department to be infected with both SARS-CoV-2 and a second viral pathogen (Editor's note in [27] ). In our study, for some of the samples diagnosed as positive for SARS-CoV-2 by specific real-time PCR and the RespiCoV Panel, viral co-infections with herpes simplex virus type 1 and Epstein-Barr virus could be identified, both of which are usually not included in screening of patients with respiratory tract infections. However, herpes simplex virus type 1 infection or reactivation in the lower and upper respiratory tract has been recorded in patients in intensive care and has increasingly been associated with pulmonary diseases with poor outcome [19, 28] . Although quantification with the RespiCoV Panel is not validated, low read numbers of herpes simplex virus type 1 and Epstein-Barr virus indicate low viral concentration in the throat.

Furthermore, 24.7% of patients infected with H1N1 during the influenza pandemic showed co-infection with bacteria, mainly Staphylococcus aureus and Streptococcus pneumoniae [29] . S. pneumoniae has also been identified as a co-infection in patients infected with influenza during the pandemic 1918-1919 and during the Asian and Hong Kong influenza pandemics of 1957 and 1968 [30, 31] .

In direct comparison, the RespiCoV Panel was shown to be less sensitive than specific realtime PCRs for SARS-CoV-2, but able to identify SARS-CoV-2 from patient samples with a Cq up to 33. Hands-on and sequencing take several hours and costs can be higher than commercial multiplex-PCR; however, additional information about the identified pathogen, including species and strain, can be obtained by the method. Due to the generation of specific amplicons, no sequence information of the host is generated which could be conflicting with personal data protection for shotgun sequencing.

Since the ongoing outbreak of SARS-CoV-2 starting in 2019, specific real-time PCR diagnostics has been contributing to the elucidation and containment of the pandemic. However, differential diagnostics and identification of Covid-19 co-infections might contribute to health care management and provide further understanding of Covid-19 courses of diseases. With RespiCoV, we have introduced an approach of highly multiplexed PCRs and MinION sequencing which can be used for rapid and comprehensive simultaneous screening for many pathogens.

Supporting information S1 

The 2009 A (H1N1) influenza virus pandemic: A review

SARS and MERS: Recent insights into emerging coronaviruses

Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR

Development of Genetic Diagnostic Methods for Novel Coronavirus 2019 (nCoV-2019) in Japan

Molecular Diagnosis of a Novel Coronavirus (2019-nCoV) Causing an Outbreak of Pneumonia

Resource-efficient internally controlled in-house real-time PCR detection of SARS-CoV-2

Co-infections of SARS-CoV-2 with multiple common respiratory pathogens in infected patients

SARS-CoV-2 respiratory co-infections: incidence of viral and bacterial co-pathogens

Clinical diagnosis of 8274 samples with 2019-novel coronavirus in Wuhan. medRxiv

Prevalence of Co-Infections with Respiratory Viruses in Individuals Investigated for SARS-CoV-2 in Ontario

Human metapneumovirus detection in patients with severe acute respiratory syndrome

SARS and common viral infections

Prevalence and outcomes of co-infection and superinfection with SARS-CoV-2 and other pathogens: A systematic review and meta-analysis

A method to identify respiratory virus infections in clinical samples using next-generation sequencing

Nanopore metagenomics enables rapid clinical diagnosis of bacterial lower respiratory infection

Metagenomic nanopore sequencing of influenza virus direct from clinical respiratory samples

Upper respiratory tract infections

Herpes simplex virus in the respiratory tract of critical care patients: A prospective study

Human herpesviruses respiratory infections in patients with acute respiratory distress (ARDS)

Clinical Characteristics and Outcomes in Hospitalized Patients with Respiratory Viral Co-Infection during the 2009 H1N1 Influenza Pandemic

Mixed Respiratory Virus Infections

Dual Infection of Infants by Human Metapneumovirus and Human Respiratory Syncytial Virus Is Strongly Associated with Severe Bronchiolitis

Multiple versus single virus respiratory infections: Viral load and clinical disease severity in hospitalized children

Circulation of other respiratory viruses and viral co-infection during the 2009 pandemic influenza

Co-infection with SARS-CoV-2 and Influenza A Virus in Patient with Pneumonia

Co-infection with SARS-CoV-2 and Human Metapneumovirus

Clinical impact of HSV-1 detection in the lower respiratory tract from hospitalized adult patients

The impact of bacterial and viral coinfection in severe influenza

Deaths from bacterial pneumonia during 1918-19 influenza pandemic

Predominant Role of Bacterial Pneumonia as a Cause of Death in Pandemic Influenza: Implications for Pandemic Influenza Preparedness

We kindly thank Ursula Erikli for copy-editing and Ute Kramer for sequencing.