Respiratory virus detection and sequencing from negative SARS-CoV-2 rapid antigen tests

Introduction

The SARS-CoV-2 pandemic highlighted the importance of genomic epidemiology in understanding virus transmission and evolution, informing essential countermeasures from non-pharmaceutical interventions to vaccines. Massive global efforts in SARS-CoV-2 genomic surveillance were made possible by widespread diagnostic testing and the growth of new infrastructure and methods for sequencing and analysis (1). Most genomic surveillance pipelines in the U.S. obtained residual SARS-CoV-2 positive samples from clinical, public health, and commercial laboratories. This strategy was effective during the pandemic but difficult to maintain with the rise of at-home rapid antigen tests (2, 3). As traditional sample sources declined, our group and others demonstrated that residual samples from rapid antigen tests could be used to generate and analyze full SARS-CoV-2 sequences for genomic surveillance (4-6).

Here, we build upon this work by identifying, sequencing, and analyzing other respiratory viruses using residual swab samples from negative BinaxNow™ COVID-19 antigen tests. This multi-virus approach is important as SARS-CoV-2 has transitioned to an endemic virus whose symptoms resemble those of other respiratory viruses (7). Thus, there is both a need for broad testing and an opportunity to expand genomic surveillance for respiratory viruses using self-collected samples.

Methods

Detailed laboratory and analysis methods are provided in the Appendix. Briefly, participants were enrolled in a parent study evaluating novel viral diagnostic tests through the RADx program at the Atlanta Center for Microsystems Engineered Point-of-Care Technologies. The study protocol was approved by the Emory Institutional Review Board and the Grady Research Oversight Committee. We performed RNA metagenomic sequencing as described (8), obtaining a median of 5.8 million reads per sample (Supplementary Data). We used a three-step bioinformatic approach to detect viruses (Supplementary Figure 1) using KrakenUniq, blastn, and reference mapping, with a final criterion requiring coverage of at least 3 distinct genome regions, based on clinical diagnostic criteria for metagenomic sequencing (9).

Results

We collected negative BinaxNOW™ test samples from 53 individuals between April-August 2023 (Supplementary Table 1), a period during which 68% of the BinaxNOW™ tests in the parent study were negative. All individuals were symptomatic at the time of testing (Table 1), and the median interval between symptom onset and testing was 2 days (range 0-9). RT-PCR was positive for influenza B in three samples and negative for influenza A and SARS-CoV-2 in all samples (Supplementary Data).

Metagenomic sequencing identified a low level of SARS-CoV-2 in one sample and a different pathogenic human respiratory virus in 17 of the other 52 samples (33%) (Supplementary Data). The following viruses were detected: parainfluenza viruses (N=7), rhinoviruses (N=5), influenza B (N=3), seasonal coronaviruses (N=2), and adenovirus (N=1) (Figure 1). In one sample, both influenza B and parainfluenza 2 were detected. In another sample positive for influenza B by RT-PCR, metagenomic sequencing did not identify influenza but identified human mastadenovirus E. Thus, excluding SARS-CoV-2, a total of 18 viruses were detected across 17 samples. There was no difference in the total number of reads obtained for samples with and without viruses detected (Mann Whitney U test, p=0.29).

• Download figure
• Open in new tab

Figure 1: Frequency of human pathogenic respiratory viruses found in 53 residual samples from BinaxNOW™ tests that were negative for SARS-CoV-2.

Pie charts indicate the number of samples positive for each virus among all samples (left panel) and among the 18 positive samples (right panel). Numbers indicate the number of samples with a virus identified, followed in parentheses by the number of samples with a >90% complete genome sequence assembled.

We observed potential differences in symptom frequencies between individuals with and without viruses detected, but none were statistically significant (Table 1).

Of the 18 viruses detected, we generated full viral genome sequences from 11 (61%), with >90% coverage and 71-24,000 fold depth (Supplementary Data). These included parainfluenza 3 (4/4 samples), parainfluenza 2 (1/2), rhinovirus (5/5), and influenza B (1/3).

We performed phylogenetic analysis of parainfluenza 3 as a proof-of-concept for genomic epidemiology studies and found substantial diversity. Using the lineage classification system described in (10), two of our sequences clustered with Lineage A1 sequences from 2019-2023 (Figure 2A), another clustered with Lineage C sequences from Japan in 2023, and the fourth with Lineage C sequences from the U.S. collected between 2015-2017 (Figure 2B), all with high bootstrap support (Supplementary Figure 2). Of note, there are only about 450 complete parainfluenza 3 virus sequences available; the data from our small study represent nearly 1% of this number, underscoring the opportunity to easily expand genomic surveillance using this approach.

• Download figure
• Open in new tab

Figure 2: Phylogenetic analysis of parainfluenza 3 virus sequences.

The names of sequences obtained in this study are bold and in red, and reference sequences are in black. The outer ring indicates virus lineage: A) contains representative lineage A1 sequences, and B) contains representative sequences from lineages C, E, F, and G. Each tree is a maximized parsimony subtree using down sampled data from the full analysis in Supplementary Figure 2, for ease of visualization.

In addition to human pathogenic respiratory viruses, we detected over 100 viruses of no clinical significance, including bacteriophages and plant viruses, many of which were also detected in our negative controls (Figure 3). Similarly, mastadenovirus C was found in many samples and negative controls. These are all consistent with environmental or reagent contaminants. Herpesviruses were found in many samples by KrakenUniq and blastn, but were not confirmed by mapping to a reference sequence with coverage of at least 3 regions. Overall, 1,367 viral taxa were identified by KrakenUniq, only 254 (18.6%) were confirmed by BLAST, and only 137 (53.9% of these, 10% of total) met our criteria for detection, highlighting the importance of confirmatory steps in metagenomic analysis.

• Download figure
• Open in new tab

Figure 3: Plot of the viral taxa (rows) that were detected in each sample (columns).

Blue boxes indicate viruses that were detected by both KrakenUniq and blastn but not confirmed by reference mapping, while red boxes indicate viruses that were detected by both KrakenUniq and blastn and were confirmed by reference mapping (covering at least 3 distinct regions of the reference virus genome).

Discussion

Our study demonstrates that RNA metagenomic sequencing of residual swab samples from negative BinaxNOW™ tests can be used to detect a broad range of respiratory viruses, including rhinoviruses, parainfluenza viruses, influenza B, seasonal coronaviruses, and adenovirus. All of these have overlapping symptoms with one another and with SARS-CoV-2, underscoring the need for multi-virus testing approaches. Although our study was not designed for clinical diagnosis, metagenomic sequencing is increasingly used clinically, and our results illustrate the need for rigorous analysis techniques and careful interpretation.

It is notable that only 33% of samples had a human pathogenic respiratory virus. This is similar to our prior study detecting alternative respiratory viruses in only 40% of SARS-CoV-2 negative individuals using residual clinical samples early in the pandemic (8). Possible explanations include individuals with a non-infectious syndrome, a bacterial or other non-viral infection, or a virus present at a low level. It is also possible that some individuals were infected with a DNA virus not optimally captured by RNA sequencing. However, we detected adenovirus, the most prevalent respiratory DNA virus. Among common RNA viruses, we did not detect influenza A or RSV, which we attribute to the winter-predominant seasonality of these viruses compared to our sample collection in spring and summer.

Importantly, of the 18 viruses detected, we were able to generate full viral genome sequences from 11 (61%) using moderate sequencing depths. Thus, the single laboratory technique of metagenomic sequencing can not only identify diverse respiratory viruses but also contribute to their genomic surveillance. The surprisingly high depth of genome coverage achieved for many sequences indicates that throughput and cost can be improved by reducing total sequencing reads from each sample in future studies.

By combining metagenomic sequencing with the use of residual antigen test samples, we demonstrate a mechanism for convenient and broad respiratory virus surveillance. Our study used BinaxNOW™ tests, which conveniently preserve the used swab within the kit cassette; future work is needed to evaluate this approach using rapid antigen test strips themselves, as previously demonstrated for SARS-CoV-2 sequencing (5). Additionally, future studies would benefit from a regulatory framework in which results can be returned to study participants, who are likely curious about the presence of other respiratory viruses when rapid antigen testing is negative.

In conclusion, our study illustrates that residual samples from self-collected antigen tests can be a powerful sample source for investigating the genomic epidemiology of a broad range of respiratory viruses, building upon the strong foundations for viral surveillance established during the SARS-CoV-2 pandemic.

Data Availability

All raw sequencing data (cleaned of human reads) is available in NCBI SRA under BioProject PRJNA1144955, and assembled virus genome sequences are available in NCBI GenBank with accession numbers listed in the Supplementary Data file.

Disclosures

All authors report no conflicts of interest to disclose.

Funding

This work was supported by NIH U54 EB027690 02S1, U54 EB027690 03S1, U54EB027690 03S2 UL1 TR002378 and the Centers for Disease Control and Prevention-funded Georgia Pathogen Genomics Center of Excellence contract 40500-050-23234506. This study was supported in part by the Emory Integrated Genomics Core (EIGC) (RRID:SCR_023529), which is subsidized by the Emory University School of Medicine and is one of the Emory Integrated Core Facilities. Additional support was provided by the Georgia Clinical & Translational Science Alliance of the National Institutes of Health under Award Number UL1TR002378. The content is solely the responsibility of the authors and does not necessarily reflect the official views of the National Institutes of Health.

Author bio

Ms. Jules received a Bachelor of Science in Anthropology and Human Biology from Emory University and is currently a research specialist in the Department of Pathology and Laboratory Medicine in the Emory University School of Medicine. She will be applying to medical school with the aspiration of becoming a family doctor and expanding healthcare to underserved communities.

Figures

• Download figure
• Open in new tab

Supplementary Figure 1: Analysis pipeline for metagenomic classification and confirmation.

Red boxes describe processing and quality control steps, yellow boxes describe initial metagenomic classification using KrakenUniq, green boxes describe blastn confirmation of reads classified as viral by KrakenUniq, and blue boxes describe reference-based mapping for final confirmation. LCA = lowest common ancestor.

Supplementary Figure 2: Maximum likelihood phylogenetic analysis of parainfluenza 3 virus sequences. The names of sequences obtained in this study are bold and in red, and reference sequences in black represent all unique full-length genome sequences of parainfluenza 3 available in GenBank (7/30/24). Circles indicate nodes with >95% ultrafast bootstrap support. The outer ring indicates virus lineage.