Metagenomic Analysis of Respiratory RNA Virome of Children with and without Severe Acute Respiratory Infection from the Free State, South Africa during COVID-19 Pandemic Reveals Higher Diversity and Abundance in Summer Compared with Winter Period

Ogunbayo, Ayodeji E.; Mogotsi, Milton T.; Sondlane, Hlengiwe; Nkwadipo, Kelebogile R.; Sabiu, Saheed; Nyaga, Martin M.

doi:10.3390/v14112516

Open AccessArticle

Metagenomic Analysis of Respiratory RNA Virome of Children with and without Severe Acute Respiratory Infection from the Free State, South Africa during COVID-19 Pandemic Reveals Higher Diversity and Abundance in Summer Compared with Winter Period

by

Ayodeji E. Ogunbayo

¹,

Milton T. Mogotsi

¹

,

Hlengiwe Sondlane

¹,

Kelebogile R. Nkwadipo

¹,

Saheed Sabiu

²

and

Martin M. Nyaga

^1,*

¹

Next Generation Sequencing Unit and Division of Virology, Faculty of Health Sciences, University of the Free State, Bloemfontein 9300, South Africa

²

Department of Biotechnology and Food Science, Durban University of Technology, P.O. Box 1334, Durban 4000, South Africa

^*

Author to whom correspondence should be addressed.

Viruses 2022, 14(11), 2516; https://doi.org/10.3390/v14112516

Submission received: 24 October 2022 / Revised: 4 November 2022 / Accepted: 11 November 2022 / Published: 14 November 2022

(This article belongs to the Special Issue Pediatric Respiratory Viral Infection)

Download

Browse Figures

Versions Notes

Abstract

:

Viral respiratory infections contribute to significant morbidity and mortality in children. Currently, there are limited reports on the composition and abundance of the normal commensal respiratory virome in comparison to those in severe acute respiratory infections (SARIs) state. This study characterised the respiratory RNA virome in children ≤ 5 years with (n = 149) and without (n = 139) SARI during the summer and winter of 2020/2021 seasons in South Africa. Nasopharyngeal swabs were, collected, pooled, enriched for viral RNA detection, sequenced using Illumina MiSeq, and analysed using the Genome Detective bioinformatic tool. Overall, Picornaviridae, Paramoxyviridae, Pneumoviridae, Picobirnaviridae, Totiviridae, and Retroviridae families were the most abundant viral population in both groups across both seasons. Human rhinovirus and endogenous retrovirus K113 were detected in most pools, with exclusive detection of Pneumoviridae in SARI pools. Generally, higher viral diversity/abundance was seen in children with SARI and in the summer pools. Several plant/animal viruses, eukaryotic viruses with unclear pathogenicity including a distinct rhinovirus A type, were detected. This study provides remarkable data on the respiratory RNA virome in children with and without SARI with a degree of heterogeneity of known viruses colonizing their respiratory tract. The implication of the detected viruses in the dynamics/progression of SARI requires further investigations.

Keywords:

severe acute respiratory infection; virome; metagenomics; children; winter; summer

1. Introduction

Severe acute respiratory infection (SARI) is a leading cause of paediatric hospitalisation and mortality [1,2]. Respiratory RNA viruses are mostly implicated, with potential to cause pandemics [3]. Despite intensive laboratory investigations, a substantial proportion of acute respiratory infections are of unknown aetiology [4,5]. Beside the known viral agents often implicated in acute, symptomatic respiratory infections, recent analysis of the human respiratory virome has found hitherto undescribed viruses, viruses with unclear pathogenicity, viruses that induces symptoms but are relatively uncommon respiratory tract pathogens bacteriophages, and retroviral elements [6,7,8,9,10,11,12,13]. The virome composition is also reported to vary in response to environment, temperature/humidity, age, and immune status; and may be distinct in health and disease [12,14].

Recent studies suggests that the respiratory microbiome is linked to airway health and may influence the development of upper and lower respiratory tract illnesses (RTIs) [15,16,17]. Characterising detectable viral populations in the human respiratory tract is essential for comprehending the role of the respiratory virome to diseases affecting the respiratory tract, as this may reveal information on a large number of undetectable virus-induced RTIs [17,18]. Additionally, it may lead to the discovery of new viruses, the identification of various viral variations, and discovery of other viruses with hitherto undiscovered tropisms [17,18,19]. Moreover, analysis of the human respiratory virome may reveal differences in viral species between patients with high/low disease severity and healthy individuals [12,14,20]. Additionally, comprehensive data of viruses present in the respiratory tract of children during different weather seasons may provide baseline information on which pathogenic viruses are predominant in various weather conditions [19]. The knowledge garnered can guide seasonal surveillance, optimal schedule for prophylactics through vaccination, and the rational inventory and use of antivirals [21,22,23]. All of which can assist in prevention and control of RTIs.

The upscaled use of metagenomics next-generation sequencing (mNGS) has increasingly provided the invaluable comprehensive genomic profile of diverse micro-organisms present in clinical samples [14,24,25]. Beside detailed virome characterization, mNGS immediately offers additional information on virulence markers, epidemiology, and molecular genotyping [25]. Although the results from investigations of children’s respiratory virome in health and/or disease states have been reported [9,14,20,26,27,28], there is little information on the RNA respiratory virome and its composition during different weather seasonality in children with and without SARI. Therefore, using clinically informative data, this study provides pioneering comprehensive data on the respiratory virome of children in health and disease during varied weather conditions in the Free State Province of South Africa, with a view to afford baseline data that can guide further studies seeking to expand/understand viral respiratory diseases in children.

2. Materials and Methods

2.1. Study Settings and Demography

The patients were recruited from Botshabelo District Hospital, Pelonomi Regional Hospital, and National District Hospital in the Free State Province, South Africa. The patient population are children ≤ 5 years of age admitted with SARI (as per the World Health Organization (WHO) definition of SARI) and without SARI. The WHO case definitions for SARI were implemented as follows: children presenting with acute respiratory infection with a history of fever or measured fever of ≥38 °C and cough; with onset within the last ten days and requiring hospitalisation [29]. For this study, non-SARI children were defined as any child ≤ 5 years admitted/attended at the hospital for other forms of illnesses without any reported/documented history of illness with signs/symptoms relating to the respiratory tract in the last 14 days and at the time of admission/sample collection.

2.2. Sampling and Sample Collection

The sampling for this study was performed during two seasons, summer 2020/2021 and the winter of 2021. The summer sampling started from January to March 2020 and then commenced again between December 2020 to March 2021. This was due to national lockdown measures implemented during the COVID-19 pandemic. The winter sampling was performed between May 2021 to August 2021. In total, (n = 105 with SARI versus n = 80 without SARI) and (n = 44 with SARI versus n = 59 without SARI) were enrolled during the summer and winter periods, respectively. Upon admission, qualified medical personnel collected nasopharyngeal swabs (BD Diagnostics, Franklin Lakes, NJ, USA). After collection, the swab was immediately inserted into a viral transport media (VTM) (BD Diagnostics, Franklin Lakes, NJ, USA). Samples were labelled and transported to the University of the Free State-Next Generation Sequencing Unit (UFS-NGS Unit) Bloemfontein, Free State, South Africa, via cold chain transportation. Subsequently, the samples were registered, and stored at −80 °C until processing.

2.3. Sample Pooling

The respiratory samples collected per child (nasopharyngeal/oropharyngeal swabs) were processed as per in-house protocol. Briefly, the samples in VTM BD Diagnostics (Franklin Lakes, NJ, USA) were vortexed for 12 s and processed in pools using 1000 µL of each sample. For the summer samples, n = 105 children with SARI were processed into five pools (5 pools of n = 3, n = 8, n = 23, n = 55, and n = 16) and n = 80 children without SARI were processed into four pools (4 pools of n = 6, n = 16, n = 38, and n = 20 samples per pool). In winter, n = 44 children with SARI were processed into three pools (3 pools of n = 15, n = 14, and n = 15 samples per pool), and n = 59 children without SARI were also processed into three pools (3 pools of n = 20, n = 20, and n = 19 samples per pool). Beside the weather-guided sampling, pooling was further made based on the different waves of COVID-19 for the summer samples and based on monthly collection for the winter samples (considering the third wave was experienced through the winter periods). The adopted pooling strategy is presented in Supplementary Table S1.

2.4. Sample Processing, Enrichment and Extraction

During the summer season, the SARI group was processed into five pools (A1 to E1), and the non-SARI group was processed into four pools (A2 to D2). During the winter period, the SARI group was processed into three pools (E1 to G1), and the non-SARI group was processed into three pools (E2 to G2).

For each pool, the VTM was centrifuged at 10,000× g for 10 min to remove cellular debris. The supernatant was filtered through a 0.22 µm filter to remove the remaining possible host cellular debris and bacteria. The resulting filtrate was treated with a nuclease mixture of 0.1U µL⁻¹ Turbo DNAse (Life Technologies, Carlsbad, CA, USA), 0.1U µL⁻¹ RNAse One (Promega, Fitchburg, WI, USA), and 1X DNAse buffer (Life Technologies, Carlsbad, CA, USA) and incubation at 37 °C for 90 min, to remove non-incorporating nucleic acid. Nucleic acid extraction was performed using the PureLink viral RNA/DNA mini kit (Thermofisher Scientific, Waltham, MA, USA), following the manufacturer’s instructions but without using the carrier RNA. Extracted viral RNA was eluted in 50 µL of RNase-free water.

2.5. Positive and Negative Controls for Evaluation of Viral Recovery and Workflow Contamination

Clinical respiratory samples positive for RNA viruses (SARS-CoV-2 and rhinovirus A) (the viral transport media (VTM) tested positive for the different viruses during routine diagnostic testing and were provided by the National Health Laboratory Service, Universitas Bloemfontein, Free State, South Africa) and no-template control (Nuclease-free water) were subjected to nucleic acid extraction and processed through the same mNGS workflow as the pooled samples. The positive controls were included to evaluate the viral recovery of the workflow, and the negative control was included to assess the presence of cross-contamination or kitome contamination from reagents.

2.6. Qubit Quantification of the Extracted RNA

Extracted RNA was quantified using the Qubit RNA High Sensitivity (HS) Assay Kit (Thermofisher Scientific, Waltham, MA, USA) with the Qubit^® 3.0 Fluorometer (Thermofisher Scientific, Waltham, MA, USA) as per the manufacturer’s instructions.

2.7. DNase Treatment on Extracted RNA and Purification

Extracted RNA was treated with TURBO DNA-free™ reagents (Thermofisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. The DNase-treated samples were purified using the RNeasy Mini Kit (Qiagen, Hilden, Germany) as per the manufacturer’s instructions.

2.8. Depletion of Ribosomal RNA (rRNA)

For each pool, rRNA depletion was performed to decrease the human, mouse, and rat rRNA by using NEBNext rRNA depletion kit (New England Biolabs, Ipswich, MA, USA) according to the manufacturer’s instructions.

2.9. Reverse Transcription and Whole Transcriptome Amplification

For library preparation purposes, the starting material was first converted to DNA. Therefore, the enriched, isolated viral RNA was reverse transcribed to generate cDNA. For the cDNA synthesis and random amplification, a QIASeq FX Single Cell RNA Library Preparation Kit (Qiagen, Hilden, Germany) was used as per the manufacturer’s instructions. Amplified cDNA was subsequently quantified using the Qubit™ 1X dsDNA HS Assay Kits with the Qubit 3.0 Fluorometer (Thermofisher Scientific, Waltham, MA, USA) as per the manufacturer’s instructions and normalised.

2.10. Library Preparation and Next-Generation Sequencing

Libraries were prepared using the QIASeq FX Single Cell RNA Library Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The DNA libraries were analysed for fragment size distribution using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Metagenomic sequencing was performed on the Illumina MiSeq system with the reagent kit v3 (Illumina, San Diego, CA, USA) for 600 cycles to generate 2 × 250 bp paired-end reads.

2.11. Bioinformatic Analysis

The generated paired-end reads were analysed using Genome Detective (https://www.genomedetective.com/, (accessed on 22 February 2022)), an automated web system for virus identification from high-throughput sequencing data [30,31]. The characteristic of this tool is presented in Supplementary Table S2. Confirmatory analyses were further performed for detections with low viral reads using BLAST.

2.12. Phylogenetic Analysis

To analyse the genetic variation in selected respiratory viruses, nucleotide sequences with close to full genome coverage detected in this study (at least 99%) were compared with hits entries in the GenBank database. Two Human rhinovirus (HRV) representative sequences, HRV A C-series and HRV A B-series, with 99.2% and 99.8% genome coverage, respectively, were selected. For the HRV analysis, the 100 closest hits full-length genomes of HRV A were retrieved for analysis, and a full-length genome of HRV B was also retrieved for use as an outgroup. Phylogenetic analysis was conducted using the neighbor-joining method using MEGA 6.0 with a bootstrap value of 1000 [32].

2.13. Statistical Analysis

Continuous variables are reported as median, whereas categorical variables are reported as percentages and calculated using Microsoft excel version 365 version 2206.

3. Results

3.1. General Demographics and Clinical Presentations

The average age for the SARI and non-SARI groups were 14.5 and 19 months, respectively. A total of 27 (18.1%) and 12 (8.5%) in the SARI and non-SARI group had HIV infection, respectively. Notably, none of the SARI group required the need for intensive care unit. Complete demographic and clinical data are presented as Supplementary Table S3.

3.2. Sequencing Data

After quality filtering, the total number of reads per pool during summer season ranged from 140,444 to 1,703,176. The proportion of viral reads in each pool ranged from 2% to 12% (Table 1). During winter season, the total number of reads per sample pool after quality filtering ranged from 13,308 to 293,394 and the proportion of virus-specific reads in each pool ranged from 2% to 9% (Table 1). It was noteworthy that the highest percentage of viral reads were observed in the non-SARI groups across the two investigated seasons (Table 1).

3.3. Respiratory Tract Virome Analyses from Children with and without SARI during Summer and Winter

Bioinformatic analyses revealed diverse viruses in the respiratory tracts of children with SARI across the two seasons investigated (Table 2 and Table 3, respectively). During summer season, viral families with the most abundant viral reads can be grouped into: Picornaviridae, Retroviridae, Totiviridae, Picobirnaviridae, and Pneumoviridae (Figure 1). During winter, the viral family with the most abundant reads can be grouped into Picornaviridae, Pneumoviridae, and Picobirnaviridae (Figure 2). Notably, the Betaflexiviridae, Partitiviridae, Chrysoviridae, Herpesviridae, Phycodnaviridae, and Circoviridae were exclusively detected in the summer period in this group.

The viral composition in the non-SARI group also contains diverse viral families (Table 4 and Table 5). The viral families with the most abundant reads in this group during summer can be grouped into Picornaviridae, Totiviridae, Paramoxyviridae, and Retroviridae (Figure 3). During winter, the viral families with the most reads in this group were mainly Retroviridae, Bromoviridae, and Reoviridae (Figure 4). Notably, the Reoviridae and Bromoviridae viral families were only detected in this group during winter period. Comparatively, Totiviridae, Coronaviridae, Partitiviridae, Picobirnaviridae, Endornaviridae, and Herpesviridae were exclusively detected during the summer period.

3.4. Viral Composition during the Seasons Investigated and between Both SARI and Non-SARI Group

Based on viral reads, mammalian eukaryotic viruses were more abundant than plant and other viruses in both groups during the two seasons (Figure 5 and Figure 6). More so, HRV and human endogenous retrovirus K113 were the only viruses detected in most pools from the SARI and non-SARI groups during both the summer and winter seasons. Comparatively, regardless of having SARI or not, children in the summer group (both SARI and Non-SARI) had higher viral reads and abundance than the participants recruited during the winter period. Notably, the Pneumoviridae viral family (RSV) was only detected in the SARI group during both periods investigated. Moreover, viral reads for Human Rotavirus and Human parainfluenza viruses 2 and 3 were exclusively detected in the non-SARI group. Besides detection of reads matching human and eukaryotic viruses, several phage-related reads (matching Escherichia phage, Enterobacteria phage, Streptococcus phage, Klebsiella phage, Staphylococcus phage, Stx2-converting phage, Lactococcus phage, Vibrio phage, Salmonella phage, Proteus phage, and Aggregatibacter phage) were detected from both groups (Supplementary Table S4).

3.5. Exclusive Viral Reads Detection in Pools

Beside the detection of RSV exclusively in patients with SARI. Its exclusive detection was also noted between second and third waves upward (Table 2). Rotavirus was also exclusively detected only in the July/August Non-SARI samples (Table 5). Of interest, the only SARS-CoV-2 reads detected were from the pools between COVID-19 lockdown declarations through the first wave.

3.6. Viral Recovery and Evaluation for Contaminants in the mNGS Workflow

As presented in supplementary data (Table S5 and Figure S1a,b), the recovery and distribution of reads on the viruses included in the run as positive control suggest that RNA viral genomes can be successfully obtained through the workflow adopted in this study. Notably, only bacteriophages were detected at less than 0.6% genome coverage in the no template control included in the sequencing run.

3.7. Phylogenetic Analysis

The phylogenetic analysis for the contigs built identified a distinct type of human rhinovirus A named (RvA-Cseries) which formed a distinct cluster with the 100 closest hits on GenBank (mostly sequences isolated in 2021), while RvA-Bseries formed a moderate cluster with the closest hits on GenBank. The RvA-Bseries series exhibited nucleotide/amino acid similarities with the RvA-Cseries ranging of 70.6%; 64.7%. In comparison with the 100 hits from GenBank, the RvA-Bseries displayed nucleotide and amino acid identities of 95.7–96.0% and 98.0–98.3% respectively, while the RvA-Cseries exhibited nucleotide and amino acid identities of 70.6–72.7% and 65.9–66.1%, respectively. The 100 HRV A hits from GenBank are collapsed in Figure 7 but complete data are provided in Table S3. The two complete genomes of HRV A (RvA-Cseries SARI and RvA-Bseries SARI), and one partial genome of HRV C (RvC-Cseries SARI) sequences were logged in the GenBank (Accession numbers: OP114090; OP114091; OP114092).

4. Discussion

This study demonstrated how mNGS is an important tool to decipher the respiratory RNA virome with significant depth. The diversity of detected reads originating from ~59 species in the pooled samples in this study mirrors the relative ease of exposure of the human respiratory system to several organisms. The majority of known respiratory RNA viruses detected in this study belonged to four families, namely Picornaviridae, Pneumoviridae, Coronaviridae, and Paramoxyviridae, and interestingly, these families have all been previously detected in the human respiratory tract [14,18,27,28]. The detected reads for several plants and animal viruses and other human non-respiratory viruses in this study for which pathogenicity in the respiratory tract is yet to be established can be either regarded as temporary commensals, or could have been introduced from the environment, contact with animals, food, and water [26,28]. While this study further reports on viral reads for several RNA viral families/species which have not been reported elsewhere, some viruses reported from other studies were not detected. For instance, some studies [14,33] detected a high level of viruses in the family Anelloviridae, which were not detected in the current study. The differences noted in viral detection in different studies may be attributed to differences in geographical location, viral exposure, or methodological choices [18].

In terms of virome composition in children with and without SARI, regardless of the sampling period (winter or summer), the proportion and abundance of viruses (in terms of family diversity and reads number) in the SARI group were higher than that of the control group. This observation is consistent with previous studies [14,20] on the virome composition in febrile versus afebrile children. Specifically, a study by Wylie et al. [20] reported higher viral sequences in the nasopharyngeal swabs samples from febrile children compared with afebrile ones. Wang et al. [14] also reported higher viral abundance in the respiratory tracts of children with SARI compared with children without SARI. Similar to other metagenomic studies [14,28], and PCR-based studies [34,35], the samples from the subjects in this study without SARI also contains reads from known epidemic respiratory viruses, such as HRV A, B, and C, Coronavirus NL63, PIV 2, and 3. The presence of these viruses in children without SARI may be transient, or the particular viral strain does not induce clinical symptoms, especially with HRV, which has been reported to be commonly present in the respiratory tracts of young children [20,36,37,38,39,40]. Moreover, the detection of the known pathogenic viruses in the non-SARI group further makes it challenging to use the virome data in establishing/judging the possible causative agents (in the SARI group) in this study as previously performed [14]. Nonetheless, observations from this study highlight the need for comprehensive studies of the dynamics of the airway virome over time to evaluate the possibility of a transient infection becoming symptomatic and to understand the future contribution of asymptomatic viral infection to community transmission.

Regarding the contribution of HIV to virome composition, a previous study reported alteration of microbiome composition in individuals with advanced HIV infection [41], while Beck et al. [42] reported that the respiratory microbiome measured in whole bronchoalveolar lavage (BAL) was indistinguishable between an HIV-infected and uninfected population. Similarly, Monaco and colleagues reported that in the absence of immunodeficiency, HIV has a minimal effect on the enteric DNA virome and bacterial microbiome. Rather, AIDS and the resultant immunodeficiency were associated with notable alterations [43]. In the current study, even though samples were pooled regardless of HIV status, the HIV-infected children were on antiretroviral therapy (ART), and as such their inclusion is not expected to contribute to alteration in the pooled virome composition towards increased viral detection.

Furthermore, the detection of sequence reads belonging to the human endogenous retrovirus K113 (HERV-K) in this study can be attributed to the recent entry of this group into the human genome. The group includes numerous retroviruses with full-length intact proviruses [44]. Typically, these proviruses only occasionally express in patients with cancer or autoimmune disease [45,46]. However, a recent study by Ferravente and colleagues used NPS to characterise the respiratory virome in SARS-CoV-2 patients; they detected HERV-K in 10 patients who were characterised with severe outcomes [47]. Similarly, Temerozo et al. [13] characterised the virome of tracheal aspirates of severe COVID-19 patients and linked the detection of HERV-K in their lower airways with early mortality. In this study, the detection of reads for the HERV-K113 in the sample pools for both SARI and non-SARI cases during winter and summer cannot be exclusively associated with severe outcomes, especially in patients in the SARI group, and as such the detection of this provirus in both groups may warrant further studies.

The detection of HRV in most sample pools in both SARI and non-SARI cases throughout the sampling period (winter and summer) further corroborates reports from another study in South Africa [48] and other countries [49,50,51] that HRV is an all-year-round virus that may circulate in both winter and summer. The detection was reportedly maintained or heightened post-COVID-19-targeted NPIs [52,53]. In addition, we identified a distinct type of HRV A, which from the phylogenetic analysis formed a distinct cluster with the closest hits on GenBank. These closest hits were majorly samples detected and deposited in the year 2021 from the USA. Despite being the closest 100 hits on GenBank, the detected HRV A strain only had a nucleotide/amino acid identity similarity that ranged from 70.6–72.7% to 65.9–66.1% with these hits. This detection may warrant future genomic surveillance of HRV in the study setting. Although, the presence of RSV in healthy children has been reported previously [14,18,28], in this study, reads for RSV were exclusively detected only in the SARI pools. More so, the total absence of the seasonal influenza virus was noted and can be due to the sampling period (during COVID-19 waves), where the complete absence/limited presence of influenza viral circulation has been reported [54,55].

In terms of viral diversity/abundance during the summer versus winter period, both the SARI and non-SARI groups had a relatively higher viral diversity and abundance (in terms of reads number) of known respiratory viruses and other viruses during the summer than the winter period. The differences in viral abundance can be due to the number of children recruited during the winter period, as fewer children presented to the hospital with SARI and non-SARI cases during this time, resulting in fewer participants per pool, thus less viral abundance. It can also be due to viral exposure, considering the sampling period, where several non-pharmaceutical interventions (NPIs) such as wearing of mask and constant hand sanitization were in place. Moreover, there was a build-up of social restrictions alongside other NPIs into the winter period which may have affected viral exposure and circulation, especially of respiratory viruses. Hence, there was a low pocket of circulating and limited detection of respiratory viruses in the winter period.

Beside the weather-based sampling, the samples were also pooled based on pre and during waves of COVID-19 experienced in the country (for the summer samples) and based on months (for the winter samples). The rationale was to evaluate compositional changes, if any, in the respiratory virome prior to COVID-19 declaration, over the different COVID-19 waves and different winter months. Notably, no RSV reads were detected in the SARI sample pools prior to COVID-19 lockdown declaration through the first wave. Reads for RSV were only detected from second wave upwards. This agrees with a previous study [56] where the surge in reported cases of RSV were noted from summer 2020 (during second wave). Furthermore, the only reads detected for SARS-CoV-2 were from the pools between COVID-19 declaration through first wave. Reads for SARS-CoV-2 would be expected in at least more than one pool considering the study period; however, relatively fewer reads for SARS-CoV-2 were detected. From our recent reports on the pathogen profile of children with SARI using a multiplex real-time PCR in the same study settings and period, comparatively higher Ct values (average of ≥33) were generally noted for SARS-CoV-2 (inversely correlated with low viral load) in the child population [57]. If this is also the case across the child population in this study in those possibly infected with COVID-19, then suggestively SARS-CoV-2 viral reads recovery could have been impacted or even masked by the more abundant host reads and could have resulted in missed detection by mNGS [58,59,60], hence fewer reads for SARS-CoV-2. Moreover, Ct values of ≤30 were previously reported for optimal detection of SARS-CoV-2 virus in whole genome sequencing [61].

A major limitation of the study was the inability to attribute the occurrence of each virus to individual participants due to the pooling of samples. In this regard, a larger individual-based longitudinal sampling can be performed. This would allow a comprehensive elucidation of the length of stay of these viruses in the respiratory tract in health and disease, and the existence or not of a permanent viral community in the respiratory niche in each patient. Moreover, relatively fewer participants were recruited during the winter period, a situation beyond control as fewer admissions were seen for both SARI and Non-SARI cases. Furthermore, this study only focused on enrichment for RNA viruses and thus cannot account for missed DNA viruses; thereby necessitating the need for further studies on the DNA respiratory virome in children with and without SARI. Lastly, using a 0.22 µm filter may have resulted in the loss of larger viruses that could have been detected. Additionally, the non-SARI group included diverse groups of patients, with varying conditions including meningitis, neurological disorders, and diarrhoea and may thus present with detectable viruses (such as enterovirus and rotavirus) in their nasopharynx. In this regard, some virome characteristics described in the non-SARI group may be confounded by the diversity of some disorders in these patients.

5. Conclusions

This study provided remarkable pioneering data of the RNA respiratory virome of children with and without SARI in South Africa during winter and summer when the COVID-19 pandemic was at the peak before vaccination interventions. The diverse distribution of viral reads in each sample pool, highlights similarities and differences both within and between children with and without SARI. Despite the study’s exploratory nature, it raised concerns about whether some of the viruses detected with unknown pathogenicity may exacerbate clinical course or contribute to symptoms manifestation. It also raises the question of whether the several known epidemic viruses detected in the non-SARI group are transient or biomarkers for future respiratory infection. Moreover, the detection of a distinct HRV A type in this study highlights the need for continuous genomic surveillance of the respiratory virome for possible detection and characterization of other distinct/novel strains which may have the potential for outbreak.

Of importance, considering the study sampling period, the virome composition could have been impaired due to various NPIs against COVID-19 such as restriction of social gatherings and wearing of mask; all of which significantly impacted viral circulation of respiratory viruses. Consequently, there may be a need for other studies to evaluate the spectrum/composition of respiratory virome in health and disease post-COVID-19; of which this study would serve as an invaluable baseline. Lastly, the detection of expected/unexpected pathogenic viruses in both SARI and non-SARI groups, the detection of viruses with potential to worsen the course of respiratory infection, and the detection of other pathogenic viruses with unknown roles in respiratory infection further highlights the complexity of this niche in health and disease and may contribute further to the understanding and management of SARI in children.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14112516/s1, Table S1: Pooling criteria for the SARI and non-SARI samples in winter and summer; Table S2: Characteristics of the Genome Detective analysis tool; Table S3: Demographic and clinical data of children recruited for the study; Table S4: Detected bacteriophages across all sample pools in both seasons; Table S5: Viral recovery and coverage of the included positive controls; Figure S1: Viral genome coverage and sequencing depth as retrieved from Genome Detective web result. The two representative viruses are (a) SARS-Co-V-2 (b) rhinovirus A; Figure S2: Complete accession numbers of the 100 rhinovirus A sequences retrieved from GenBank.

Author Contributions

Conceptualization, M.M.N., S.S. and A.E.O.; Supervision, and Project administration, M.M.N. and S.S.; sample delivery, A.E.O., M.T.M. and K.R.N.; Data curation, A.E.O., K.R.N. and H.S.; Investigation, A.E.O., M.T.M. and H.S.; Validation, A.E.O.; Methodology, A.E.O.; Formal analysis, A.E.O.; Visualization and Original draft, A.E.O.; Writing—review and editing, all co-authors; Funding acquisition, M.M.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Bill and Melinda Gates Foundation, grant number BMGF-OPP1180423_2017, the National Research Foundation (NRF) South Africa, grant number 120814, and the Poliomyelitis Research Foundation (PRF) South Africa, grant number 19/16.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Health Science Research Ethics Committee (HSREC) of the University of the Free State, South Africa (Approval no UFS-HSD2019/1129/2910 and approval date of 24 October 2019).

Informed Consent Statement

Informed consent was obtained from the parents of all participants involved in the study.

Data Availability Statement

The data presented in this study are available in the supplementary information.

Acknowledgments

We would like to thank the following senior nurses, Palesa Kortman, Late Yolanda Porota, and Mahlapenyane Botsane for their assistance with the patient recruitment and sample collection. Our appreciation also goes to Peter Nthiga Mwangi for his assistance with the phylogenetic analysis. We also appreciate the National Health Laboratory Service, Universitas, Virology Section, for providing us with the positive controls used in validation of the workflow.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

Gessner, B.D. Severe acute respiratory illness in Sub-Saharan Africa. J. Infect. Dis. 2015, 212, 843–844. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ayeni, O.A.; Walaza, S.; Tempia, S.; Groome, M.; Kahn, K.; Madhi, S.A.; Cohen, A.L.; Moyes, J.; Venter, M.; Pretorius, M.; et al. Mortality in children aged <5 years with severe acute respiratory illness in a high HIV- prevalence urban and rural areas of South Africa, 2009–2013. PLoS ONE 2021, 16, e0255941. [Google Scholar] [CrossRef]
Tang, J.W.T.; Lam, T.T.; Zaraket, H.; Lipkin, W.I.; Drews, S.J.; Hatchette, T.F.; Heraud, J.M.; Koopmans, M.P.; Abraham, A.M.; Baraket, A.; et al. Global epidemiology of non-influenza RNA respiratory viruses: Data gaps and a growing need for surveillance. Lancet Infect. Dis. 2017, 17, e320–e326. [Google Scholar] [CrossRef]
Aman, A.T.; Wibawa, T.; Kosasih, H.; Asdie, R.H.; Safitri, I.; Intansari, U.S.; Mawarti, Y.; Sudarmono, P.; Arif, M.; Puspitasari, D.; et al. Etiologies of severe acute respiratory infection (SARI) and misdiagnosis of influenza in Indonesia, 2013–2016. Influenza Other Respir. Viruses 2021, 15, 34–44. [Google Scholar] [CrossRef] [PubMed]
Anders, K.L.; Nguyen, H.L.; Nguyen, N.M.; Van Thuy, N.T.; Hong Van, N.T.; Hieu, N.T.; Hong Tham, N.T.; Thanh Ha, P.T.; Lien, L.B.; Chau, N.V.V.; et al. Epidemiology and virology of acute respiratory infections during the first year of life: A birth cohort study in Vietnam. Pediatr. Infect. Dis. J. 2015, 34, 361–370. [Google Scholar] [CrossRef] [Green Version]
Van Den Hoogen, B.G.; De Jong, J.C.; Groen, J.; Kuiken, T.; De Groot, R.; Fouchier, R.A.M.; Osterhaus, A.D.M.E. A newly discovered human pneumovirus isolated from young children with respiratory tract disease. Nat. Med. 2001, 7, 719–724. [Google Scholar] [CrossRef]
Ksiazek, T.G.; Erdman, D.; Goldsmith, C.S.; Zaki, S.R.; Peret, T.; Emery, S.; Tong, S.; Urbani, C.; Comer, J.A.; Lim, W.; et al. A novel coronavirus associated with severe acute respiratory syndrome. N. Engl. J. Med. 2003, 348, 1953–1966. [Google Scholar] [CrossRef]
Allander, T.; Tammi, M.T.; Eriksson, M.; Bjerkner, A.; Tiveljung-Lindell, A.; Andersson, B. Cloning of a human parvovirus by molecular screening of respiratory tract samples. Proc. Natl. Acad. Sci. USA 2005, 102, 12891–12896. [Google Scholar] [CrossRef] [Green Version]
Zhou, Y.; Fernandez, S.; Yoon, I.K.; Simasathien, S.; Watanaveeradej, V.; Yang, Y.; Marte-Salcedo, O.A.; Shuck-Lee, D.J.; Thomas, S.J.; Hang, J.; et al. Metagenomics study of viral pathogens in undiagnosed respiratory specimens and identification of human enteroviruses at a Thailand hospital. Am. J. Trop. Med. Hyg. 2016, 95, 663–669. [Google Scholar] [CrossRef] [Green Version]
Jartti, T.; Lehtinen, P.; Vuorinen, T.; Österback, R.; Van Den Hoogen, B.; Osterhaus, A.D.M.E.; Ruuskanen, O. Respiratory picornaviruses and respiratory syncytial virus as causative agents of acute expiratory wheezing in children. Emerg. Infect. Dis. 2004, 10, 1095–1101. [Google Scholar] [CrossRef]
Young, N.S.; Brown, K.E. Parvovirus B19. N. Engl. J. Med. 2004, 350, 586–597. [Google Scholar] [CrossRef] [PubMed]
Zárate, S.; Taboada, B.; Yocupicio-Monroy, M.; Arias, C.F. Human virome. Arch. Med. Res. 2017, 48, 701–716. [Google Scholar] [CrossRef] [PubMed]
Temerozo, J.R.; Fintelman-Rodrigues, N.; dos Santos, M.C.; Hottz, E.D.; Sacramento, C.Q.; de Paula Dias da Silva, A.; Mandacaru, S.C.; dos Santos Moraes, E.C.; Trugilho, M.R.O.; Gesto, J.S.M.; et al. Human endogenous retrovirus K in the respiratory tract is associated with COVID-19 physiopathology. Microbiome 2022, 10, 65. [Google Scholar] [CrossRef] [PubMed]
Wang, Y.; Zhu, N.; Li, Y.; Lu, R.; Wang, H.; Liu, G.; Zou, X.; Xie, Z.; Tan, W. Metagenomic analysis of viral genetic diversity in respiratory samples from children with severe acute respiratory infection in China. Clin. Microbiol. Infect. 2016, 22, 458.e1–458.e9. [Google Scholar] [CrossRef] [Green Version]
Louten, J. Virus transmission and epidemiology. In Essential Human Virology; Academic Press: Cambridge, MA, USA, 2016; pp. 71–92. [Google Scholar]
Man, W.H.; Piters, W.A.A.D.S.; Bogaert, D. The microbiota of the respiratory tract: Gatekeeper to respiratory health. Nat. Rev. Microbiol. 2017, 15, 259–270. [Google Scholar] [CrossRef]
Tay, C.J.X.; Ta, L.D.H.; Ow Yeong, Y.X.; Yap, G.C.; Chu, J.J.H.; Lee, B.W.; Tham, E.H. Role of upper respiratory microbiota and virome in childhood rhinitis and wheeze: Collegium Internationale Allergologicum Update 2021. Int. Arch. Allergy Immunol. 2021, 182, 265–276. [Google Scholar] [CrossRef]
Wylie, K.M. The virome of the human respiratory tract. Clin. Chest Med. 2017, 38, 11–19. [Google Scholar] [CrossRef]
Romero-espinoza, J.A.; Moreno-valencia, Y.; Coronel-tellez, R.H.; Castillejos-lopez, M.; Hernandez, A.; Dominguez, A.; Miliar-garcia, A.; Barbachano-guerrero, A.; Perez-padilla, R.; Alejandre-garcia, A.; et al. Virome and bacteriome characterization of children with pneumonia and asthma in Mexico City during winter seasons 2014 and 2015. PLoS ONE 2018, 13, e0192878. [Google Scholar] [CrossRef] [Green Version]
Wylie, K.M.; Mihindukulasuriya, K.A.; Sodergren, E.; Weinstock, G.M.; Storch, G.A. Sequence analysis of the human virome in Febrile and Afebrile children. PLoS ONE 2012, 7, e27735. [Google Scholar] [CrossRef] [Green Version]
Famoroti, T.; Sibanda, W.; Ndung, T. Prevalence and seasonality of common viral respiratory pathogens, including Cytomegalovirus in children, between 0–5 years of age in KwaZulu-Natal, an HIV endemic province in South Africa. BMC Pediatr. 2018, 18, 240. [Google Scholar] [CrossRef]
Khor, C.S.; Sam, I.C.; Hooi, P.S.; Quek, K.F.; Chan, Y.F. Epidemiology and seasonality of respiratory viral infections in hospitalized children in Kuala Lumpur, Malaysia: A retrospective study of 27 years. BMC Pediatr. 2012, 12, 32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Lei, C.; Yang, L.; Lou, C.T.; Yang, F.; SiTou, K.I.; Hu, H.; Io, K.; Cheok, K.T.; Pan, B.; Ung, C.O.L. Viral etiology and epidemiology of pediatric patients hospitalized for acute respiratory tract infections in Macao: A retrospective study from 2014 to 2017. BMC Infect. Dis. 2021, 21, 306. [Google Scholar] [CrossRef] [PubMed]
Parker, J.; Chen, J. Application of next generation sequencing for the detection of human viral pathogens in clinical specimens. J. Clin. Virol. 2017, 86, 20–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
van Boheemen, S.; van Rijn, A.L.; Pappas, N.; Carbo, E.C.; Vorderman, R.H.P.; Sidorov, I. Retrospective validation of a metagenomic sequencing protocol for combined detection of RNA and DNA viruses using respiratory samples from pediatric patients. J. Mol. Diagn. 2020, 22, 196–207. [Google Scholar] [CrossRef] [Green Version]
Rajagopala, S.V.; Bakhoum, N.G.; Pakala, S.B.; Shilts, M.H.; Rosas-Salazar, C.; Mai, A.; Boone, H.H.; McHenry, R.; Yooseph, S.; Halasa, N.; et al. Metatranscriptomics to characterize respiratory virome, microbiome, and host response directly from clinical samples. Cell Rep. Methods 2021, 1, 100091. [Google Scholar] [CrossRef]
Lysholm, F.; Wetterbom, A.; Lindau, C.; Darban, H.; Bjerkner, A.; Fahlander, K.; Lindberg, A.M.; Persson, B.; Allander, T.; Andersson, B. Characterization of the viral microbiome in patients with severe lower respiratory tract infections, using metagenomic sequencing. PLoS ONE 2012, 7, e30875. [Google Scholar] [CrossRef]
Li, Y.; Fu, X.; Ma, J.; Zhang, J.; Hu, Y.; Dong, W.; Wan, Z.; Li, Q.; Kuang, Y.Q.; Lan, K.; et al. Altered respiratory virome and serum cytokine profile associated with recurrent respiratory tract infections in children. Nat. Commun. 2019, 10, 228. [Google Scholar] [CrossRef] [Green Version]
Fitzner, J.; Qasmieh, S.; Mounts, A.W.; Alexander, B.; Besselaar, T.; Briand, S.; Brown, C.; Clark, S.; Dueger, E.; Gross, D.; et al. Revision of clinical case definitions: Influenza-like illness and severe acute respiratory infection. Bull. World Health Organ. 2018, 96, 122–128. [Google Scholar] [CrossRef]
Vilsker, M.; Moosa, Y.; Nooij, S.; Fonseca, V.; Ghysens, Y.; Dumon, K.; Pauwels, R.; Alcantara, L.C.; Vanden Eynden, E.; Vandamme, A.M.; et al. Genome Detective: An automated system for virus identification from high-throughput sequencing data. Bioinformatics 2019, 35, 871–873. [Google Scholar] [CrossRef] [Green Version]
Cleemput, S.; Dumon, W.; Fonseca, V.; Karim, W.A.; Giovanetti, M.; Alcantara, L.C.; Deforche, K.; Oliveira, T. Genome detective coronavirus typing tool for rapid identification and characterization of novel coronavirus genomes. Bioinformatics 2020, 36, 3552–3555. [Google Scholar] [CrossRef]
Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Lewandowska, D.W.; Schreiber, P.W.; Ruehe, B.; Zagordi, O.; Bayard, C.; Greiner, M.; Benden, C.; Geissberger, D.; Capaul, R.; Zbinden, A.; et al. Metagenomic sequencing complements routine diagnostics in identifying viral pathogens in lung transplant recipients with unknown etiology of respiratory infection. PLoS ONE 2017, 12, e0177340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Moe, N.; Pedersen, B.; Nordbø, S.A.; Skanke, L.H.; Krokstad, S.; Smyrnaios, A.; Døllner, H. Respiratory virus detection and clinical diagnosis in children attending day care. PLoS ONE 2016, 11, e0159196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Douros, K.; Kotzia, D.; Kottaridi, C.; Giotas, A.; Boutopoulou, B.; Bozas, E.; Matziou, V.; Priftis, K.; Papaevangelou, V. Evidence for respiratory viruses interactions in asymptomatic preschool-aged children. Allergol. Immunopathol. 2019, 47, 260–264. [Google Scholar] [CrossRef]
Advani, S.; Sengupta, A.; Forman, M.; Valsamakis, A.; Milstone, A. Detecting respiratory viruses in asymptomatic children. Pediatr. Infect. Dis. J. 2012, 23, 1221–1226. [Google Scholar] [CrossRef] [Green Version]
Walker, G.J.; Stelzer-Braid, S.; Shorter, C.; Honeywill, C.; Wynn, M.; Willenborg, C.; Barnes, P.; Kang, J.; Pierse, N.; Crane, J.; et al. Viruses associated with acute respiratory infection in a community-based cohort of healthy New Zealand children. J. Med. Virol. 2022, 94, 454–460. [Google Scholar] [CrossRef] [Green Version]
Castro, I.A.; Costa, L.D.C.; Oliveira, A.C.R.; Souza, M.; das Dôres de Paula Cardoso, D.; Camargos, P.A.M.; Costa, P.S.S.; Fiaccadori, F.S. Circulation profile of respiratory viruses in symptomatic and asymptomatic children from Midwest Brazil. Brazilian J. Microbiol. 2020, 51, 1729–1735. [Google Scholar] [CrossRef]
Sarna, M.; Lambert, S.B.; Sloots, T.P.; Whiley, D.M.; Alsaleh, A.; Mhango, L.; Bialasiewicz, S.; Wang, D.; Nissen, M.D.; Grimwood, K.; et al. Viruses causing lower respiratory symptoms in young children: Findings from the ORChID birth cohort. Thorax 2018, 73, 969–979. [Google Scholar] [CrossRef] [Green Version]
Morikawa, S.; Hiroi, S.; Kase, T. Detection of respiratory viruses in gargle specimens of healthy children. J. Clin. Virol. 2015, 64, 59–63. [Google Scholar] [CrossRef]
Twigg, H.L.; Knox, K.S.; Zhou, J.; Crothers, K.A.; Nelson, D.E.; Toh, E.; Day, R.B.; Lin, H.; Gao, X.; Dong, Q.; et al. Effect of advanced HIV infection on the respiratory microbiome. Am. J. Respir. Crit. Care Med. 2016, 194, 226–235. [Google Scholar] [CrossRef]
Beck, J.M.; Schloss, P.D.; Venkataraman, A.; Twigg, H.; Jablonski, K.A.; Bushman, F.D.; Campbell, T.B.; Charlson, E.S.; Collman, R.G.; Crothers, K.; et al. Multicenter comparison of lung and oral microbiomes of HIV-infected and HIV-uninfected individuals. Am. J. Respir. Crit. Care Med. 2015, 192, 1335–1344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Monaco, C.L.; Gootenberg, D.B.; Zhao, G.; Handley, S.A.; Ghebremichael, M.S.; Lim, E.S.; Lankowski, A.; Baldridge, M.T.; Wilen, C.B.; Flagg, M.; et al. Altered virome and bacterial microbiome in human immunodeficiency virus-associated acquired immunodeficiency syndrome human immunodeficiency virus (HIV) infection is associated with increased intestinal translocation of microbial products and enteropathy as well as alterations in gut bacterial. Cell Host Microbe 2016, 19, 311–322. [Google Scholar] [CrossRef] [PubMed]
Moyes, D.L.; Martin, A.; Sawcer, S.; Temperton, N.; Worthington, J.; Griffiths, D.J.; Venables, P.J. The distribution of the endogenous retroviruses HERV-K113 and HERV-K115 in health and disease. Genomics 2005, 86, 337–341. [Google Scholar] [CrossRef] [PubMed]
Johanning, G.L.; Malouf, G.G.; Zheng, X.; Esteva, F.J.; Weinstein, J.N.; Wang-Johanning, F.; Su, X. Expression of human endogenous retrovirus-K is strongly associated with the basal-like breast cancer phenotype. Sci. Rep. 2017, 7, 41960. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Boller, K.; Schönfeld, K.; Lischer, S.; Fischer, N.; Hoffmann, A.; Kurth, R.; Tönjes, R.R. Human endogenous retrovirus HERV-K113 is capable of producing intact viral particles. J. Gen. Virol. 2008, 89, 567–572. [Google Scholar] [CrossRef]
Ferravante, C.; Sanna, G.; Melone, V.; Fromentier, A.; Rocco, T.; D’Agostino, Y.; Lamberti, J.; Alexandrova, E.; Pecoraro, G.; Pagliano, P.; et al. Nasopharyngeal virome analysis of COVID-19 patients during three different waves in Campania region of Italy. J. Med. Virol. 2022, 94, 2275–2283. [Google Scholar] [CrossRef]
Pretorius, M.A.; Tempia, S.; Treurnicht, F.K.; Walaza, S.; Cohen, A.L.; Moyes, J.; Hellferscee, O.; Variava, E.; Dawood, H.; Chhagan, M.; et al. Genetic diversity and molecular epidemiology of human rhinoviruses in South Africa. Influenza Other Respi. Viruses 2014, 8, 567–573. [Google Scholar] [CrossRef]
Jiang, H.; Yang, T.; Yang, C.; Lu, Y.; Yi, Z.; Zhang, Q.; Wang, W. Molecular epidemiology and clinical characterization of human rhinoviruses circulating in Shanghai, 2012–2020. Arch. Virol. 2022, 167, 1111–1123. [Google Scholar] [CrossRef]
Moriyama, M.; Hugentobler, W.J.; Iwasaki, A. Seasonality of respiratory viral infections. Annu. Rev. Virol. 2020, 7, 83–101. [Google Scholar] [CrossRef]
Toivonen, L.; Schuez-Havupalo, L.; Karppinen, S.; Teros-Jaakkola, T.; Rulli, M.; Mertsola, J.; Waris, M.; Peltola, V. Rhinovirus infections in the first 2 years of life. Pediatrics 2016, 138, e20161309. [Google Scholar] [CrossRef]
Oh, D.Y.; Buda, S.; Biere, B.; Reiche, J.; Schlosser, F.; Duwe, S.; Wedde, M.; von Kleist, M.; Mielke, M.; Wolff, T.; et al. Trends in respiratory virus circulation following COVID-19-targeted nonpharmaceutical interventions in Germany, January–September 2020: Analysis of national surveillance data. Lancet Reg. Health-Eur. 2021, 6, 100112. [Google Scholar] [CrossRef] [PubMed]
Varela, F.H.; Sartor, I.T.S.; Polese-Bonatto, M.; Azevedo, T.R.; Kern, L.B.; Fazolo, T.; de David, C.N.; Zavaglia, G.O.; Fernandes, I.R.; Krauser, J.R.M.; et al. Rhinovirus as the main co-circulating virus during the COVID-19 pandemic in children. J. Pediatr. 2022, 98, 579–586. [Google Scholar] [CrossRef]
Tempia, S.; Walaza, S.; Bhiman, J.N.; McMorrow, M.L.; Moyes, J.; Mkhencele, T.; Meiring, S.; Quan, V.; Bishop, K.; McAnerney, J.M.; et al. Decline of influenza and respiratory syncytial virus detection in facility-based surveillance during the COVID-19 pandemic, South Africa, January to October 2020. Eurosurveillance 2021, 26, 2001600. [Google Scholar] [CrossRef] [PubMed]
Groves, H.E.; Papenburg, J.; Mehta, K.; Bettinger, J.A.; Sadarangani, M.; Halperin, S.A.; Morris, S.K.; Bancej, C.; Burton, C.; Embree, J.; et al. The effect of the COVID-19 pandemic on influenza-related hospitalization, intensive care admission and mortality in children in Canada: A population-based study. Lancet Reg. Health-Am. 2022, 7, 100132. [Google Scholar] [CrossRef] [PubMed]
Bents, S.; Viboud, C.; Grenfell, B.; Hogan, A.; Tempia, S.; von Gottberg, A.; Moyes, J.; Walaza, S.; Cohen, C.; Baker, R. The impact of COVID-19 non-pharmaceutical interventions on future respiratory syncytial virus transmission in South Africa. medRxiv 2022, 12, 22271872. [Google Scholar] [CrossRef]
Ogunbayo, A.E.; Mogotsi, M.T.; Sondlane, H.; Nkwadipo, K.R.; Sabiu, S.; Nyaga, M.M. Pathogen profile of children hospitalised with severe acute respiratory infections during COVID-19 pandemic in the Free State Province, South Africa. Int. J. Environ. Res. Public Health 2022, 19, 10418. [Google Scholar] [CrossRef] [PubMed]
Babiker, A.; Bradley, H.L.; Stittleburg, V.D.; Ingersoll, J.M.; Key, A.; Kraft, C.S.; Waggoner, J.J.; Piantadosi, A. Metagenomic sequencing to detect respiratory viruses in persons under investigation for COVID-19. J. Clin. Microbiol. 2021, 59, e02142-20. [Google Scholar] [CrossRef]
Mostafa, H.H.; Fissel, J.A.; Fanelli, B.; Bergman, Y.; Gniazdowski, V.; Dadlani, M.; Carroll, K.C.; Colwell, R.R.; Simner, P.J. Metagenomic next-generation sequencing of nasopharyngeal specimens collected from confirmed and suspect covid-19 patients. mBio 2020, 11, e01969-20. [Google Scholar] [CrossRef]
de Vries, J.J.C.; Brown, J.R.; Fischer, N.; Sidorov, I.A.; Morfopoulou, S.; Huang, J.; Munnink, B.B.O.; Sayiner, A.; Bulgurcu, A.; Rodriguez, C.; et al. Benchmark of thirteen bioinformatic pipelines for metagenomic virus diagnostics using datasets from clinical samples. J. Clin. Virol. 2021, 141, 104908. [Google Scholar] [CrossRef]
Cloete, J.; Kruger, A.; Masha, M.; du Plesis, N.M.; Mawela, D.; Tshukudu, M.; Manyane, T.; Komane, L.; Venter, M.; Jassat, W.; et al. Paediatric hospitalisations due to COVID-19 during the first SARS-CoV-2 omicron (B.1.1.529) variant wave in South Africa: A multicentre observational study. Lancet Child Adolesc. Health 2022, 6, 294–302. [Google Scholar] [CrossRef]

Figure 1. An overview of representative viral family distribution from children with SARI in summer based on read abundance. The viral families shown include, Picobirnaviridae, Picornaviridae, Pneumoviridae, Totiviridae, and Retroviridae. Others comprise of Partitiviridae, Betaflexiviridae, Coronaviridae, Chrysoviridae, Virgaviridae, Herpesviridae, Phycodnaviridae, Endornaviridae, and Circoviridae.

Figure 2. An overview of representative viral family distribution from children with SARI in winter based on read abundance. Viral families with the most abundance include Picornaviridae, Pneumoviridae, and Picobirnaviridae, Retroviridae, and Totiviridae. Others include Endornaviridae, Coronaviridae, Virgaviridae, and Tombusviridae.

Figure 3. An overview of representative viral family distribution from the non-SARI group in summer based on read abundance. The major viral families detected include Picornaviridae, Totiviridae, Paramoxyviridae, Retroviridae, Picobirnaviridae, and Endornaviridae. Other viral families at low abundance include Coronaviridae, Partitiviridae, Chrysoviridae, and Herpesviridae.

Figure 4. An overview of representative viral family distribution from the non-SARI group in winter based on read abundance. The viral families include Retroviridae, Bromoviridae, Reoviridae, and Picornaviridae, and Retroviridae.

Figure 5. Different types of viruses from the SARI group in winter and summer. Mammalian virus were the most predominant followed by fungal, plant, and a small fraction of other viruses (avian viruses, yeast viruses, and algae viruses).

Figure 6. Different types of viruses from the non-SARI group in winter and summer. Mammalian viruses were predominant, followed by yeast virus, plant virus, and a small fraction of avian, fungal, and algae viruses.

Figure 7. Phylogenetic analysis of rhinovirus A strains (B series and C series SARI) with 100 hit HRV-A strains from the GenBank. The complete 100 rhinovirus A hits strains were collapsed and are available in supplementary Figure S2.

Table 1. Reads for each pooled samples after filtering with the percentage of viral reads.

Pooled Samples	Number of Samples	Reads after Quality Filtering	% of Viral Reads
SARI Summer
SARI Pool A1	3	1,079,760	2
SARI Pool B1	8	1,703,176	3
SARI Pool C1	23	8,500,586	4
SARI Pool D1	55	1,341,384	5
SARI Pool E1	16	862,754	11
Non-SARI Summer
Non-SARI Pool A2	6	140,444	5
Non-SARI Pool B2	16	312,496	12
Non-SARI Pool C2	38	366,146	4
Non-SARI Pool D2	20	1,250,024	6
SARI Winter
SARI Pool E1	15	293,394	2
SARI Pool F1	14	166,726	2
SARI Pool G1	15	272,596	3
Non-SARI Winter
Non-SARI Pool E2	20	15,434	5
Non-SARI Pool F2	20	13,308	6
Non-SARI Pool G2	19	37,010	9

SARI: severe acute respiratory infection.

Table 2. Viral family, species, reads numbers, and genome structure of viruses detected in the respiratory sample pools of children with SARI during summer period.