Skip to main content

Unique microbial landscape in the human oropharynx during different types of acute respiratory tract infections



Secondary bacterial infections and pneumonia are major mortality causes of respiratory viruses, and the disruption of the upper respiratory tract (URT) microbiota is a crucial component of this process. However, whether this URT dysbiosis associates with the viral species (in other words, is viral type-specific) is unclear.


Here, we recruited 735 outpatients with upper respiratory symptoms, identified the infectious virus types in 349 participants using multiplex RT-PCR, and profiled their upper respiratory microbiome using the 16S ribosomal RNA gene and metagenomic gene sequencing. Microbial and viral data were subsequently used as inputs for multivariate analysis aimed at revealing viral type-specific disruption of the upper respiratory microbiota. We found that the oropharyngeal microbiota shaped by influenza A virus (FluA), influenza B virus (FluB), respiratory syncytial virus (RSV), and human rhinovirus (HRV) infections exhibited three distinct patterns of dysbiosis, and Veillonella was identified as a prominent biomarker for any type of respiratory viral infections. Influenza virus infections are significantly correlated with increased oropharynx microbiota diversity and enrichment of functional metabolic pathways such as L-arginine biosynthesis and tetracycline resistance gene tetW. We used the GRiD algorithm and found the predicted growth rate of common respiratory pathogens was increased upon influenza virus infection, while commensal bacteria, such as Streptococcus infantis and Streptococcus mitis, may act as a colonization resistance to the overgrowth of these pathogens.


We found that respiratory viral infections are linked with viral type-specific disruption of the upper respiratory microbiota, particularly, influenza infections uniquely associated with increased microbial diversity and growth rates of specific pathogens in URT. These findings are essential for clarifying the differences and dynamics of respiratory microbiota in healthy participants and acute respiratory viral infections, which contribute to elucidating the pathogenesis of viral-host-bacterial interactions to provide insights into future studies on effective prevention and treatment of respiratory tract infections.

Video Abstract


Respiratory viruses, such as influenza (Flu) and human respiratory syncytial virus (RSV), are major threats to global public health [1]. Respiratory viral infections predispose patients to secondary bacterial infections and are the leading cause of death from infectious diseases [2,3,4]. It is crucial to study the major causes of bacterial infections and the factors that influence symptoms or outcomes; however, the interplay of multiple factors, such as viruses, bacteria, and the host immune system, complicates the elucidation of bacterial infections [5, 6]. Current advances in microbiome research have led to a significant change in our view of bacteria: more than simply pathogens, they are a community that can live in symbiosis with the human body [7]. Microbes form a complex system (i.e., the microbiome) closely involved in human health and various conditions [8]. Viral infections disrupt the microbial-host balance maintained by this complex system [9]. Yet, this process is so complicated due to the factors involved that we are still unclear about its exact course and whether it is viral type-specific.

The respiratory microbiome is considered the gatekeeper of the respiratory tract, although its clinical implications in respiratory diseases remain largely unknown [10]. The respiratory tract can be divided into two parts, upper and lower respiratory tracts (URT and LRT, respectively), and the former has a significantly higher microbial density than the latter. In URT, the oropharynx is an important site of infection and replication for most respiratory viruses [11]. Meanwhile, the oropharyngeal microbiota is considered a major source of lung microbiota, even in healthy adults, and is involved in secondary pneumonia [10, 12]. Among viral infections, influenza virus and RSV have been well studied for their epidemics and high mortality caused by bacterial pneumonia [13,14,15]. For example, the results from a meta-analysis show that the proportion of secondary bacterial infections following influenza infection is up to 65% [16], with the majority ranging from 11.5 to 34% [14, 17,18,19]. The secondary bacterial infections in RSV infection is up to 32.4% [20,21,22], and human rhinovirus (HRV) infection ranges from 9.1 to 18.9% [23, 24]. At the same time, in case reports of influenza infection, patients with higher influenza-mediated secondary infections also have a higher disease burden [17], and nearly 30% of severe influenza patients are exacerbated by co-infection with Streptococcus pneumoniae, Staphylococcus aureus, and Haemophilus influenzae in the lung [25, 26]. With the development of microbiome studies, the similarities and differences of dysbacteriosis characteristics upon different viral infections are anticipated to explore and further characterize the intrinsic connection between this dysbiosis and secondary bacterial infections of the lower respiratory tract.

Several studies based on 16S rRNA gene and metagenomic sequencing have reported the disruption of the respiratory microbiota with influenza infection [27,28,29,30]; however, it remains unclear whether this disruption is viral type-specific due to the complexity of clinical viral infections, such as mixed infections with various viruses. This hinders further studies on how viruses act in an inflammatory environment and then influence the respiratory microbiota, contributing to a different outcome of bacterial infections in the respiratory tract. Therefore, human studies based on large sample sizes of multiple viral infections are needed to reveal the role and the mechanism of viral infections in shaping the respiratory microbiome. In the present study, we identified 349 patients with single or multiple infections of FluA, FluB, RSV, and human rhinovirus (HRV) from 735 patients using a multiplex RT-PCR approach from March to August 2019 and recruited 98 healthy participants as the control group. We collected oropharyngeal swabs (representing URT specimens) and performed 16S rRNA gene and metagenomic sequencing to profile the microbiome and elucidate the respiratory dysbiosis.


Study participants and sample preparation

This study included 735 influenza-like cases from the outpatient and 98 healthy participants from March to August 2019 at Nanhai Hospital of Southern Medical University (Foshan, Guangdong, China). Patients with three or more clinical symptoms, including fever, cough, sore throat, runny or stuffy nose, body aches or headaches, and fatigue, would be involved in this study. All participants provided written informed consent. The demographic characteristics including age, gender, height, weight, influenza vaccination records, and chronic respiratory diseases were recorded (Table 1).

Table 1 Demographics and infection characteristics of the participants included in this study

Two oropharyngeal samples were obtained from each patient for further testing: one with a FLOQSwabs (Copan, FLOQSwabs™) for RNA extraction was diluted in 3 mL of Universal Transport Medium™ (Copan, UTM™), and the other with a sterile swab (Copan, 108CSR) for DNA extraction was flashed frozen with dry ice and transferred within 1 h to a-80 °C refrigerator. The FLOQSwabs™ swab was collected first, followed by sterile swabs. Each swab was used to wipe the root of the tongue and the pharyngeal tonsils and the posterior pharyngeal wall six times on both sides of the subject.

RNA extraction, multiplex RT-PCR, and laboratory diagnosis

Total RNA was extracted with RNAiso Plus (Takara, 9109) and then detected by multiplex RT-PCR with a PrimeScriptTMRT Master Mix (Takara, RR036A) reverse transcription kit using primers and probes (Additional file 1: Table S1) targeting FluA, FluB, HRV, and RSV as previously [31]. Synthetic plasmids targeting FluA, FluB, HRV, and RSV were used to produce standard curves (Additional file 2: Fig. S1). The multiplex RT-PCR cycles were initiated at 95 °C for 5 min, followed by 40 cycles at 95 °C for 5 s, 60 °C for 30 s. Three replicate wells were set up for each sample amplification, and specimens were considered positive when the Ct value of one well was less than 37 or two wells were between 37 and 39. Among 735 influenza-like cases, 349 samples were identified as single positive or multiple positive (two or more viruses) for FluA, FluB, HRV, and RSV; the remaining 386 samples and all of the 98 healthy controls were identified as negative. Subsequently, 349 infectious patients and 98 healthy participants were applied for 16S rRNA gene sequencing. We randomly selected 10 oropharyngeal samples for metagenomic sequencing from FluA, FluB, and FluA/FluB co-infected patients and healthy participants, respectively (n = 40).

Isolation of DNA and next-generation sequencing

DNA (including negative control samples) was isolated using the Magigene DNA Isolation Kit (Magigene Ltd., Guangzhou, China) and analyzed for DNA integrity on an Agilent Bioanalyzer 2100 (Agilent Technologies, CA, USA). DNA extracted from 447 samples (349 infectious patients and 98 healthy participants), negative control samples (1 sampling negative control, 5 DNA extraction controls, and 3 PCR amplification controls), and together with the mock community as positive controls were further conducted to PCR amplification targeting the V4 region with primers (515F, 5′-GTGCCAGCMGCCGCGGTAA-3′; 806R, 5′-GGACTACHVGGGTWTCTAAT-3′) to generate the sequencing library. Amplifications and purification of PCR products were performed as in the previous study [18], and final libraries were then sequenced on the NovaSeq 6000 platform (Illumina, Inc., CA, USA) with 250bp paired-end reads generated. Overall, the average clean reads of 16S rRNA gene sequencing of specimens in this study were 236,620 (> 90,000 total reads per sample), and the average number of bacterial reads available for downstream analysis was 58,318 (> 20,000 total reads per sample) (Additional file 2: Fig. S2A-B).

A set of 40 samples were subjected to metagenomic sequencing. The libraries were generated using NEB Next® Ultra™ DNA Library Prep Kit for Illumina® (New England Biolabs, MA, USA) with index codes added. The library quality was further assessed on the Qubit 4.0 Fluorometer (Life Technologies, NY, USA) and the Qsep400 High-Throughput Nucleic Acid Protein Analysis System (Houze Biological Technology Co, Hangzhou, China) system. Finally, high-quality libraries were acquired and sequenced on a single sequencing run on the NovaSeq 6000 platform (Illumina, Inc., CA, USA) with 150 bp paired-end reads (> 3 G raw data per sample) generated. The average clean reads of metagenomic sequencing of specimens in this study were 76,576,665 (> 30 million total reads per sample), and the average number of bacterial reads available for downstream analysis was 19,938,678 (> 6 million total reads per sample) (Additional file 2: Fig. S2C).

Data analysis

16S rRNA gene analysis

Raw “fastq” files were demultiplexed and processed using tools available in QIIME 2 (v2019.7) [32] with a Divisive Amplicon Denoising Algorithm (DADA) 2-based pipeline [33]. Afterward, a feature table containing amplicon sequence variants (ASVs) and associated abundances was generated based on forwarding reads. Representative sequence sets for each dada2 sequence variant were used for taxonomy classification with a Naive Bayes classifier. Sequences were then classified against SILVA (v132) [34]. Finally, alpha- and beta-diversity measurements were performed using an even sampling depth of 20,000 sequences per sample.

Negative and positive control assessment

We performed the PCR amplification of the V4 region within 16S rRNA gene for the sampling negative control (n = 1), DNA extraction controls (n = 5), and PCR amplification controls (n = 3). PCR target fragments were not detected in all negative controls and the concentration of purified products was well below the detection limit of the Qubit 4.0 fluorometer. We continued the process with library construction, but the library failed to pass the control samples and were therefore not sequenced. One set of customized mock community was used, which contained purified isolates of two species, a gram-negative bacillus Escherichia coli and a gram-positive coccus Enterococcus faecalis. Both strains were grown in culture until the OD value was equal to 1.0; 1 ml of each bacterial solution was taken and mixed well as a mock community. Among them, the proportion of E. coli in the mock community measured by qPCR was 55.84%, and the proportion of Enterococcus faecalis was 44.16%. In our 16S rRNA gene sequencing results, the reads were classified into genera Escherichia-Shigella accounting for 56.64%, Enterococcus accounting for 40.35%, and only 3.01% of unclassified or other genus sequences outside the mock community.

Metagenomics analysis

Sequencing raw data pre-processing was performed by KneadData (v0.7.4,, which integrates the tools FastQC (v0.11.9) [35], Trimmomatic (v0.39) [36], and Bowtie2 (v2.3.5.1) [37], to do quality check, quality filtering, and host sequences decontamination, respectively. Taxonomic classification of the resulting reads was performed using MetaPhlAn2 (v2.0) [38]. Metabolic pathway abundance was determined using the HUMAnN2 (v0.9.9) [39] pipeline, and 324 microbial pathways were mapped based on MetaCyc metabolic pathway database. Characterization of antibiotic resistance genes was determined by identifying the annotated markers of the Antibiotic Resistance Database (ARDB) [40] from the metagenomic database using the ShortBRED (v0.9.4) [41] pipeline. The Growth Rate InDex (GRiD) algorithm (v1.3) [42] was used to estimate the growth rates of target bacteria. The microbes in the human respiratory tract, such as opportunistic pathogenic including Streptococcus pneumoniae and Staphylococcus aureus, and potentially pathogenic bacteria including Escherichia coli, Fusobacterium nucleatum, Haemophilus sputorum, Mycoplasma pneumoniae, Porphyromonas gingivalis, Treponema denticola, and Treponema medium, which previously reported [25, 26, 43,44,45,46,47] were summarized as pathogenic species.

Microbial diversity

Shannon diversity index was used to characterize alpha-diversity using phyloseq (v1.40.0) [48] package in R (v4.0.5) [49]. Bray-Curtis distance dissimilarities were calculated, and principal coordinates analysis (PCoA) was generated to measure the beta-diversity among samples with phyloseq (v1.40.0).

Random Forest analysis

Random Forest regression was used to predict infections based on bacterial abundances using the randomForest package (v4.7-1.1) [50], and the receiver operating characteristic curves (ROC) for each categorical group were plotted. The area under the curve (AUC) was calculated using the pROC (v1.18.0) [51] package in R (v4.0.5). In the random forest model with 100 resamplings, taxa with significant Mean Decrease Accuracy (P < 0.05) were considered “important feature,” and the Top-10 taxa with higher value of Mean Decrease Accuracy for each categorical group were applied for statistical analysis.

Microbial interaction network

Spearman was used for network inference. The co-presence and mutually exclusive networks were analyzed using the Network analyzer tool of Cytoscape (v3.8.0) [52] and outputted for further analysis, together with node and edge tables for the whole association network.

Statistical analysis

The distance matrix between microbiota communities was tested by permutational multivariate analysis of variance (PERMANOVA) using vegan (v2.6.2) [53] package in R (v4.0.5). Linear discriminant analysis (LDA) effect size (LEfSe) was conducted on Galaxy (v1.0) [54] for each microbiota paired group to determine significantly different abundances of bacterial taxa, and the parameters were set with the default P value (P < 0.05) and an LDA score of 2.0. Principal component analysis (PCA) was used to reduce the data dimensionality with ggbiplot (v0.55) [55] package. Spearman correlation analysis was performed using dplyr (v1.1.2) [56] package for network inference among the abundance of bacterial antibiotic resistance genes (ARGs), microbial pathways, and microbes. Differences between categorical groups in the abundance of specific genera were determined by Kruskal-Wallis test as well as Wilcoxon rank-sum test with subsequent Bonferroni correction, with P < 0.05 considered significant, respectively, and the P-value generated by statistics between healthy and each infectious group was used for Bonferroni correction. Kruskal–Wallis test, Wilcoxon rank-sum test, and Bonferroni were performed with ggpubr (v0.6.0) [57] and R stats (v4.2.2) [49] package, and ggplot2 (v3.4.1) [58] package was used for data visualization.


Demographic characteristics of the participants with upper respiratory virus infection

From March to August 2019, a total of 735 individuals with influenza-like symptoms were included in this study. The multiplex RT-PCR was used for FluA, FluB, HRV, or RSV verification (Additional file 1: Table S1). As a result, we found that 255 cases were single positive for any of the four viruses, and 94 cases were identified as multiple infections. Also, we recruited 98 healthy participants as controls who were reviewed as negative for viral infection (Fig. 1A, B). The demographics and infection characteristics are summarized in Table 1.

Fig. 1
figure 1

Study flow and virus infection determination. A Diagram for the study design. B Venn diagram revealed clinical coinfections of different respiratory viruses

Viral infections are associated with disorders of specific characteristics of the upper respiratory microbiota

To assess the connection of viral infections with the URT microbiota, we performed PCoA and PERMANOVA analysis based on Bray-Curtis distances to reveal the compositional difference between these microbiomes. The result showed that all types of viruses, except RSV, exhibited significant changes in the composition of the oropharyngeal microbiota compared to healthy individuals (Fig. 2A). We further used LEfSe to identify the characteristics of the individual differences in microbiota compositions between the four single virus infections and healthy individuals (Additional file 1: Table S2). We found that Veillonella was significantly enriched in the participants with either FluA, FluB, HRV, or RSV infections compared to healthy individuals. In contrast, Granulicatella was less abundant in participants with any type of viral infection than in healthy controls (Fig. 2B-D). PCA analysis on LDA values of differential groups showed that the FluA infection exhibits unique features, such as the enrichment of Solobacterium and Ralstonia, compared to other types of viral infections. Notably, the oropharyngeal microbiota with FluA infection was distinct from that in RSV and HRV groups, which showed similar patterns, whereas the change in microbial compositional characteristics upon FluB infection alone fell between these two (Fig. 2E).

Fig. 2
figure 2

Viral infections correlate to specific upper respiratory tract microbiota disorders. A PCoA based on Bray-Curtis distances was analyzed, and PERMANOVA was used to evaluate the variation of microbial composition between the virus-positive and healthy groups. B LEfSe was used to calculate the difference between genera between virus-positive and healthy groups and take the union of the genera to show the enrichment features. C, D Differences in the relative abundance of Veillonella and Granulicatella between virus-positive and healthy groups (Kruskal-Wallis test as well as Wilcoxon rank-sum test with subsequent Bonferroni correction). E PCA was performed based on the LDA values of the different genera, and Biplot was used to evaluate the association between the infectious group (arrows) and genera (points)

These findings prompted us to ponder whether the dynamic characteristics of these microbiomes could indicate viral infection or even the type of infection. We then measured the effectiveness of the oropharyngeal microbiome in distinguishing participants with specific viral infections from healthy individuals using a random-forest classification model. We first found that a microbial classifier can effectively identify an infection state. In addition, the area under the receiver operating characteristic curve (AUC) of the FluB (AUC up to 92.64%) group was more extensive than that in the FluA (AUC = 85.51%), RSV (AUC = 85.08%), and HRV (AUC = 89.28%) groups (Fig. 3A, B), indicating that our model most readily distinguishes FluB infections among the four viruses included here. To further characterize the specific species of each microbial classifier, we made a heatmap to show ASVs with significant Mean Decrease Accuracy, which can distinguish infection from non-infection (healthy state) (Fig. 3C). The ASV from genera Granulicatella, which serve as a FluA-specific biomarker, was unable to distinguish patients with FluB, RSV, and HRV from participants in this study cohort (Fig. 3C and Additional file 1: Table S3). Our results reveal the value of oropharyngeal bacterial communities for studying a wide range of clinical viral infections; of particular note is that this microbial classifier can effectively distinguish influenza infection, which highlights the association between dysbiosis and specific viral infection.

Fig. 3
figure 3

Various viral infections are associated with specific signatures of microbiome disruption. A The random forest model was used to predict the participants infected with FluA (n = 22), FluB (n = 138), RSV (n = 34), HRV (n = 61), or any virus (n = 349) from healthy (n = 98), and the ROC was calculated. B All samples of the five types were selected respectively; the selected samples and healthy samples were used for random (n = 100) under-sampling to fit a random forest model, and the AUC was calculated to evaluate the prediction performance. C The Top10-taxa with a higher value of Mean Decrease Accuracy from each categorical group was used for analysis

Influenza virus infection is linked with increased diversity and functional changes in the upper respiratory microbiota

We analyzed the microbial characteristics of populations with gender, age, BMI, and viral infection types and found that only FluA-positive patients exhibited higher alpha diversity of URT microbiota (Fig. 4A). However, only FluB-positive patients showed a significant difference in beta diversity (Fig. 4B and Additional file 2: Fig. S3). These results above revealed a closer connection between influenza virus infections and URT microbiota, rather than HRV or RSV.

Fig. 4
figure 4

Influenza virus infection is associated with increased diversity and functional changes in the upper respiratory tract microbiota. A Differences in alpha-diversity between infected and healthy participants (Wilcoxon rank-sum test). B Differences in beta-diversity between infected and healthy participants (PERMANOVA in R with “vegan” package). C The Spearman interaction network of differential pathways, differential ARGs, ARGs count, and differential species among the indicated groups were constructed

We further investigated the functional and pathogenicity characteristics of the oropharyngeal microbiota associated with influenza virus infection based on the shotgun metagenomics sequencing data. We found that compared to healthy participants, FluA alone infection and FluA/FluB co-infection groups showed enrichment for multiple functional pathways, such as biotin biosynthesis II, L-arginine biosynthesis III (via N-acetyl-L-citrulline) (Additional file 2: Fig. S4A). We also found a significantly lower diversity of drug resistance genes in FluA-positive patients than in healthy participants; nevertheless, either FluA, FluB, or co-infection groups showed higher levels of the tetracycline resistance gene tetW (Additional file 2: Fig. S4B-C). Metagenomics sequencing results supported that FluA infection is associated with the significant enrichment of Veillonella and Atopobium, including species Veillonella unclassified, Veillonella parvula, and Atopobium rimae (Additional file 2: Fig. S4D). We further built an interaction network to demonstrate the association between those microbial species that differ in influenza-infected individuals compared to healthy individuals and ARG expression and functional pathways (Fig. 4C and Additional file 1: Table S4). We found, for example, the potential link of FluA-mediated V. parvula enrichment and its positive association with the functional pathways and tetracycline resistance genes, whereas the underlying mechanism needs further investigation.

Specific respiratory pathogens exhibit higher growth rates in influenza-infected individuals

Opportunistic pathogenic bacteria of the respiratory tract, such as S. pneumoniae, H. influenzae, and S. aureus, are associated with influenza-associated secondary infection [25, 26]. The potentially pathogenic bacteria of the human respiratory tract from previous reports were summarized [43,44,45,46,47], and nine species (in the blue circles) were identified in our oropharyngeal swabs based on metagenomics sequencing analysis (Additional file 2: Fig. S5A). However, the influence of influenza infection on the relative abundance of these pathogens was modest (Additional file 2: Fig. S5B).

To better understand the potential impact of viral infection on oropharyngeal flora, the Growth Rate InDex (GRiD) algorithm was used to estimate the growth rates of pathogens, including F. nucleatum, H. influenzae, H. sputorum, P. gingivalis, T. denticola and T. medium (Additional file 2: Fig. S6 and Additional file 1: Table S5). The results further showed that the growth rate of pathogenic bacteria was significantly lower in healthy controls compared to the commensal bacteria. However, this inhibition was not observed in the infected individuals, which indicates that the influenza virus infection may create an environment that promotes the growth of pathogens (Fig. 5A). Previous studies have suggested commensal bacteria as an alternative approach to the prevention of pathogen infection due to their acting on the host immune system to induce protective responses, or to inhibit the growth of respiratory pathogens by producing antimicrobial products or signals and competing for nutrients and adhesion sites [59, 60]. Our results confirmed that several commensal species such as Streptococcus infantis, Streptococcus mitis, and Corynebacterium durum were significantly negatively correlated with the growth rate of these listed pathogens, suggesting that these commensal microbes may play a crucial role in colonization resistance (Fig. 5B).

Fig. 5
figure 5

Specific respiratory pathogens exhibit higher growth rates in patients with influenza infection. A The growth rate values of all detected potential pathogens and commensals were used to analyze the difference between infected and healthy groups (Wilcoxon rank-sum test with subsequent Bonferroni correction). B The correlations of growth rates between pathogens and commensals were analyzed, and the heatmap showed objects with significant associations. C Spearman methods revealed the significant correlations between pathogens and bacterial pathways or antibiotics genes, and the networks were displayed with Cytoscape software

Subsequently, we constructed Spearman association networks among the identified six species of pathogens, functional pathways, and drug resistance genes (|r| > 0.5, P < 0.05). The results showed that the high growth rate of T. denticola was significantly correlated to formaldehyde assimilation III (dihydroxyacetone cycle) and other pathways, and the high growth rate of T. medium correlates to the enrichment of biotin biosynthesis II, 1,4-dihydroxy-2-naphthoate biosynthesis II, while H. sputorum is related to CMP-3-deoxy-D-manno-octulosonate biosynthesis I (Fig. 5C and Additional file 1: Table S6). In conclusion, FluA-positive participants exhibited significantly enhanced growth rates of specific respiratory pathogens, and colonization resistance, multiple bacterial pathways, and drug resistance may be involved in this process. These findings revealed a close association between FluA infection with microbial function and the underlying mechanisms of secondary bacterial infections warrant further investigation.


Morbidity and mortality in severe epidemics are commonly associated with secondary bacterial infections upon viral infections, which has attracted increasing attention to the relationships between dysbiosis and influenza virus, RSV, or other respiratory acute viruses [61,62,63,64]. Previous studies of respiratory viruses and human commensal bacteria had focused on a single identified viral infection, ignoring the complexity of an actual infection [65, 66], and it remained unclear whether this dysbiosis is virus-specific. Our study represented a large cohort study and revealed the virus-specific dysbiosis in URT with complex clinical viral infection types. Intriguingly, FluA, FluB, HRV, and RSV infections formed three distinct patterns of ecological dysbiosis, indicating that microbiota dysbiosis has specific virus-dependent features. We also observed significant changes in the composition of respiratory microbiota across four viral infection types. According to the characteristics of this virus-specific microbial disorder, we established the microbial classifiers upon different infection types using a machine learning model and revealed the unique co-occurrence ASVs from Atopobium and Prevotella 7 in FluA and ASVs from Bergeyella and Prevotella aurantiaca in FluB. These unique features were uncovered to help us for a better understanding of virus pathogenesis and precision therapy. In addition, influenza A was associated with changes in the diversity in the URT microbiome in our study. A previous study reported a significant increase in oropharyngeal microbiota diversity in patients with H7N9 influenza, in particular, H7N9 patients with secondary bacterial pneumonia [12]. This suggests an underlying mechanism for influenza-associated secondary bacterial infections.

A common dysbiosis feature was observed, in which, Veillonella was significantly enriched in all types of viral infection groups compared to healthy participants. Veillonella was previously reported as the causative agent of chronic anaerobic pneumonia [67]. Furthermore, Veillonella has been shown to be over-represented in the URT, bronchoalveolar lavage fluid, or lung tissues of COVID-19 patients [68,69,70], and the occurrence of Veillonella is closely associated with the enrichment of multiple pathogens such as S. aureus and Prevotella melaninogenica [71, 72]. These results suggest that Veillonella may serve as an indicator of dysbiosis following viral infections in the URT and that we should focus on its role and potential mechanisms in terms of its pro-inflammatory capacity and bacterial interactions.

Previous studies have reported significant differences in nasopharyngeal microbial profiles between HRV and RSV infections [73], and the oropharyngeal microbial differences among FluA, FluB, HRV, and RSV infections were observed in this study, which indicate universal changes in the microbiome of the URT during acute respiratory virus infections. Recent studies on the relationship between the upper respiratory microbiome and respiratory infections have mainly relied on 16S rRNA sequencing rather than metagenomic studies [62, 66, 74], which are limited by the lack of resolution on the strain or subspecies levels as functional genes (such as drug resistance genes) and pathways. In our study, we combined metagenomic sequencing with 16S rRNA analysis to improve the resolution of the analysis. Our metagenomic data further suggested that enriched V. parvula abundance in FluA and FluA/FluB co-infected patients is associated with higher biotin synthesis capacity. The biotin synthetic pathway is considered one of the unique biological pathways in bacteria and has been proposed as a promising antimicrobial target [75]. Our previous study showed that Veillonella is able to migrate from the oral cavity to the lungs and is associated with inflammation and impaired function in the LRT [76]. Together with this study, we can rationally infer that local airway microbiota dysbiosis may be impacted by the viral infection that initially occurs in the oral cavity and subsequently impacts distant microbiota through the oral-lung axis. Therefore, manipulation of the oropharyngeal microbiome may be a selective prevention strategy.

In addition to identifying the characteristics of respiratory dysbiosis associated with different viral infections, we further proposed that the increased growth rate of pathogenic bacteria is the crucial process of influenza viral infections, which may directly mediate secondary bacterial infections in the lower respiratory tract. Bioinformatic tools such as MetaPhlAn [38] and Kraken [77] are able to accurately quantify pathogenic bacteria abundance based on metagenomic sequencing data; however, the growth rates of species, particularly pathogens, are of greater concern, especially in the early stages of secondary infection, which may directly contribute to disease progression [42, 78]. More importantly, we can better understand microbial interactions by estimating the growth rate differences between commensals and pathogens. Our study aimed to reveal the growth patterns of commensal and pathogens upon different respiratory virus infections. We found that the growth rates of pathogenic bacteria were significantly increased only in patients infected with FluA or FluB. In addition, some commensals, such as S. mitis, may activate the colonization resistance to pathogenic bacteria during viral infection. Several previous studies have found that Streptococcus spp. abundance in URT was correlated with lower susceptibility to influenza A (H3N2) and influenza B infection [79]. Among healthy young adults vaccinated with live attenuated influenza vaccine, S. infantis was positively associated with influenza H1 immunoglobulin A (IgA) titers [80]. These results suggest that some members of Streptococcus may have probiotic potential against respiratory virus infections.


We found that respiratory viral infections correlate to viral type-specific disruption of the upper respiratory microbiota, and Veillonella was identified as a prominent biomarker for any type of the four respiratory virus infections. Intriguingly, influenza infections are uniquely associated with enhanced microbial diversity and growth rates of specific pathogens in URT. We proposed that pre-changes in the oropharyngeal microbiome are potentially indicative of respiratory viruses-induced secondary bacterial infections, which provides a better understanding of complications of upper respiratory tract infection. While further investigations are needed to translate these results into potential clinical and public health applications, our findings demonstrated a critical role of microbiome perturbations in bridging the link between upper respiratory viral infections and severe secondary infections.

Availability of data and materials

The datasets supporting the conclusions of this article are available in the PRJNA801796 repository, unique persistent identifier and hyperlink to datasets at The R scripts used to generate the results in this manuscript are available in the following Git repository:


  1. Li Y, Reeves RM, Wang X, et al. Global patterns in monthly activity of influenza virus, respiratory syncytial virus, parainfluenza virus, and metapneumovirus: a systematic analysis. Lancet Glob Health. 2019;7(8):e1031–45.

    Article  PubMed  Google Scholar 

  2. Qin T, Wang Y, Deng J, et al. Super Dominant Pathobiontic Bacteria in the Nasopharyngeal Microbiota Cause Secondary Bacterial Infection in COVID-19 Patients. Microbiol Spectr. 2022;10(3):e0195621.

    Article  PubMed  Google Scholar 

  3. Morris DE, Cleary DW, Clarke SC. Secondary bacterial infections associated with influenza pandemics. Front Microbiol. 2017;8:1041.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Liu Y, Ling L, Wong SH, et al. Outcomes of respiratory viral-bacterial coinfection in adult hospitalized patients. EClinicalMedicine. 2021;37:100955.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Almand EA, Moore MD, Jaykus L-A. Virus-bacteria interactions: an emerging topic in human infection. Viruses. 2017;9(3):58.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Finlay BB, McFadden G. Anti-immunology: evasion of the host immune system by bacterial and viral pathogens. Cell. 2006;124(4):767–82.

    Article  CAS  PubMed  Google Scholar 

  7. Consortium HMP. A framework for human microbiome research. Nature. 2012;486(7402):215–21.

    Article  Google Scholar 

  8. Wypych TP, Wickramasinghe LC, Marsland BJ. The influence of the microbiome on respiratory health. Nat Immunol. 2019;20(10):1279–90.

    Article  CAS  PubMed  Google Scholar 

  9. Hanada S, Pirzadeh M, Carver KY, et al. Respiratory viral infection-induced microbiome alterations and secondary bacterial pneumonia. Front Immunol. 2018;9:2640.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Man WH, de Steenhuijsen Piters WA, Bogaert D. The microbiota of the respiratory tract: gatekeeper to respiratory health. Nat Rev Microbio. 2017;15(5):259–70.

    Article  CAS  Google Scholar 

  11. Hammitt LL, Kazungu S, Welch S, et al. Added value of an oropharyngeal swab in detection of viruses in children hospitalized with lower respiratory tract infection. J Clin Microbiol. 2011;49(6):2318–20.

  12. Lu H-F, Li A, Zhang T, et al. Disordered oropharyngeal microbial communities in H7N9 patients with or without secondary bacterial lung infection. Emerging Infect Dis. 2017;6(1):1–11.

    CAS  Google Scholar 

  13. Kilbourne ED. Influenza pandemics of the 20th century. Emerging Infect Dis. 2006;12(1):9.

    Article  Google Scholar 

  14. Collaborators GI. Mortality, morbidity, and hospitalisations due to influenza lower respiratory tract infections, 2017: an analysis for the Global Burden of Disease Study 2017. Lancet Respir Med. 2019;7(1):69–89.

    Article  Google Scholar 

  15. Nair H, Nokes DJ, Gessner BD, et al. Global burden of acute lower respiratory infections due to respiratory syncytial virus in young children: a systematic review and meta-analysis. Lancet. 2010;375(9725):1545–55.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Klein EY, Monteforte B, Gupta A, et al. The frequency of influenza and bacterial coinfection: a systematic review and meta-analysis. Influenza Other Respir Viruses. 2016;10(5):394–403.

    Article  PubMed  PubMed Central  Google Scholar 

  17. MacIntyre CR, Chughtai AA, Barnes M, et al. The role of pneumonia and secondary bacterial infection in fatal and serious outcomes of pandemic influenza a (H1N1) pdm09. BMC Infect Dis. 2018;18:1–20.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Rouzé A, Martin-Loeches I, Povoa P, et al. Relationship between SARS-CoV-2 infection and the incidence of ventilator-associated lower respiratory tract infections: a European multicenter cohort study. Intensive Care Med. 2021;47(2):188–98.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Chertow DS, Memoli MJ. Bacterial coinfection in influenza: a grand rounds review. JAMA. 2013;309(3):275–82.

    Article  CAS  PubMed  Google Scholar 

  20. Lin HC, Liu YC, Hsing TY, et al. RSV pneumonia with or without bacterial coinfection among healthy children. J Formos Med Assoc. 2022;121(3):687–93.

    Article  CAS  PubMed  Google Scholar 

  21. Satav A, Crow R, Potdar V, et al. The burden of respiratory syncytial virus in children under 2 years of age in a rural community in Maharashtra. India Clin Infect Dis. 2021;73:S238–247.

    Article  PubMed  Google Scholar 

  22. Fan R, Fan C, Zhang J, et al. Respiratory syncytial virus subtype ON1/NA1/BA9 predominates in hospitalized children with lower respiratory tract infections. J Med Virol. 2017;89(2):213–21.

    Article  CAS  PubMed  Google Scholar 

  23. Blaschke AJ, Korgenski EK, Wilkes J, et al. Rhinovirus in febrile infants and risk of bacterial infection. Pediatrics. 2018;141(2).

  24. Yuan XH, Li YM, Shen YY, Yang J, Jin Y. Clinical and Th1/Th2 immune response features of hospitalized children with human rhinovirus infection. J Med Virol. 2020;92(1):26–33.

    Article  CAS  PubMed  Google Scholar 

  25. Farias JA, Fernández A, Monteverde E, et al. Critically ill infants and children with influenza A (H1N1) in pediatric intensive care units in Argentina. Intensive Care Med. 2010;36(6):1015–22.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Estenssoro E, Ríos FG, Apezteguía C, et al. Pandemic 2009 influenza A in Argentina: a study of 337 patients on mechanical ventilation. Am J Respir Crit Care Med. 2010;182(1):41–8.

    Article  PubMed  Google Scholar 

  27. Rattanaburi S, Sawaswong V, Chitcharoen S, et al. Bacterial microbiota in upper respiratory tract of COVID-19 and influenza patients. Exp Biol and Med. 2022;247(5):409–15.

    Article  CAS  Google Scholar 

  28. Yildiz S, Mazel-Sanchez B, Kandasamy M, et al. Influenza A virus infection impacts systemic microbiota dynamics and causes quantitative enteric dysbiosis. Microbiome. 2018;6(1):1–17.

    Article  CAS  Google Scholar 

  29. Ding T, Song T, Zhou B, et al. Microbial composition of the human nasopharynx varies according to influenza virus type and vaccination status. MBio. 2019;10(4):e01296–e1319.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Qin T, Geng T, Zhou H, et al. Super-dominant pathobiontic bacteria in the nasopharyngeal microbiota as causative agents of secondary bacterial infection in influenza patients. Emerg Microbes Infect. 2020;9(1):605–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Pretorius MA, Madhi SA, Cohen C, et al. Respiratory viral coinfections identified by a 10-plex real-time reverse-transcription polymerase chain reaction assay in patients hospitalized with severe acute respiratory illness–South Africa, 2009–2010. J Infect Dis. 2012;206(Suppl 1):S159–65.

    Article  CAS  PubMed  Google Scholar 

  32. Bolyen E, Rideout JR, Dillon MR, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Callahan BJ, McMurdie PJ, Rosen MJ, et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Quast C, Pruesse E, Yilmaz P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012;41(D1):D590–6.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Andrews S. FastQC a quality control tool for high throughput sequence dat. 2016.

    Google Scholar 

  36. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Truong DT, Franzosa EA, Tickle TL, et al. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat Methods. 2015;12(10):902–3.

    Article  CAS  PubMed  Google Scholar 

  39. Franzosa EA, McIver LJ, Rahnavard G, et al. Species-level functional profiling of metagenomes and metatranscriptomes. Nat Methods. 2018;15(11):962–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Liu B, Pop M. ARDB–Antibiotic Resistance Genes Database. Nucleic Acids Res. 2009;37:D443–7.

    Article  CAS  PubMed  Google Scholar 

  41. Kaminski J, Gibson MK, Franzosa EA, et al. High-Specificity Targeted Functional Profiling in Microbial Communities with ShortBRED. PLoS Comput Biol. 2015;11(12):e1004557.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Emiola A, Oh J. High throughput in situ metagenomic measurement of bacterial replication at ultra-low sequencing coverage. Nat Commun. 2018;9(1):4956.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Smith AM, McCullers JA. Secondary bacterial infections in influenza virus infection pathogenesis. Curr Top Microbiol Immunol. 2014;385:327–56.

    PubMed  PubMed Central  Google Scholar 

  44. Urban M, Pant R, Raghunath A, et al. The Pathogen-Host Interactions database (PHI-base): additions and future developments. Nucleic Acids Res. 2015;43:D645–55.

    Article  CAS  PubMed  Google Scholar 

  45. Li ZJ, Zhang HY, Ren LL, et al. Etiological and epidemiological features of acute respiratory infections in China. Nat Commun. 2021;12(1):5026.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Zeng H, Chan Y, Gao W, et al. Diversity of Treponema denticola and Other Oral Treponeme Lineages in Subjects with Periodontitis and Gingivitis. Microbiol Spectr. 2021;9(2):e0070121.

    Article  PubMed  Google Scholar 

  47. Nørskov-Lauritsen N. Classification, identification, and clinical significance of Haemophilus and Aggregatibacter species with host specificity for humans. Clin Microbiol Rev. 2014;27(2):214–40.

    Article  PubMed  PubMed Central  Google Scholar 

  48. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.

    Google Scholar 

  50. Liaw A, Wiener M. “Classification and Regression by randomForest. R News. 2002;2(3):18–22

    Google Scholar 

  51. Robin X, Turck N, Hainard A, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14(6):927–30.

    Article  Google Scholar 

  54. Segata N, Izard J, Waldron L, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Vincent QV. ggbiplot: a ggplot2 based biplot. R package version 0.55. 2011.

  56. Wickham H, François R, Henry L, et al. dplyr: A Grammar of Data Manipulation. R package version 1.1.2. 2023.

  57. Kassambara A. ggpubr: 'ggplot2' Based Publication Ready Plots. R package version 0.6.0. 2023.

  58. Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer-Verlag; 2016. ISBN 978-3-319-24277-4.

  59. Khan R, Petersen FC, Shekhar S. Commensal Bacteria: An Emerging Player in Defense Against Respiratory Pathogens. Front Immunol. 2019;10:1203–1203.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Keith JW, Pamer EG. Enlisting commensal microbes to resist antibiotic-resistant pathogens. J Exp Med. 2019;216(1):10–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Rosas-Salazar C, Tang ZZ, Shilts MH, et al. Upper respiratory tract bacterial-immune interactions during respiratory syncytial virus infection in infancy. J Allergy Clin Immunol. 2022;149(3):966–76.

    Article  CAS  PubMed  Google Scholar 

  62. Rosas-Salazar C, Kimura KS, Shilts MH, et al. SARS-CoV-2 infection and viral load are associated with the upper respiratory tract microbiome. J Allergy Clin Immunol. 2021;147(4):1226–1233.e2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Kaul D, Rathnasinghe R, Ferres M, et al. Microbiome disturbance and resilience dynamics of the upper respiratory tract during influenza A virus infection. Nat Commun. 2020;11(1):2537.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Brundage JF, Shanks GD. Deaths from bacterial pneumonia during 1918–19 influenza pandemic. Emerging Infect Dis. 2008;14(8):1193.

    Article  Google Scholar 

  65. Grier A, Gill AL, Kessler HA, et al. Temporal Dysbiosis of Infant Nasal Microbiota Relative to Respiratory Syncytial Virus Infection. J Infect Dis. 2021;223(9):1650–8.

    Article  CAS  PubMed  Google Scholar 

  66. Ramos-Sevillano E, Wade WG, Mann A, et al. The effect of influenza virus on the human oropharyngeal microbiome. Clin Infect Dis. 2019;68(12):1993–2002.

    Article  CAS  PubMed  Google Scholar 

  67. Shah A, Panjabi C, Nair V, et al. Veillonella as a cause of chronic anaerobic pneumonitis. Int J Infect Dis. 2008;12(6):e115–7.

    Article  PubMed  Google Scholar 

  68. Battaglini D, Robba C, Fedele A, et al. The role of dysbiosis in critically ill patients with covid-19 and acute respiratory distress syndrome. Front Med. 2021;8:671714.

    Article  Google Scholar 

  69. Ma S, Zhang F, Zhou F, et al. Metagenomic analysis reveals oropharyngeal microbiota alterations in patients with COVID-19. Signal Transduct Target Ther. 2021;6(1):1–11.

    Google Scholar 

  70. Dereschuk K, Apostol L, Ranjan I, et al. Identification of Lung and Blood Microbiota Implicated in COVID-19 Prognosis. Cells. 2021;10(6).

  71. de Steenhuijsen Piters WA, Huijskens EG, Wyllie AL, et al. Dysbiosis of upper respiratory tract microbiota in elderly pneumonia patients. ISME J. 2016;10(1):97–108.

    Article  PubMed  Google Scholar 

  72. Accorsi EK, Franzosa EA, Hsu T, et al. Determinants of Staphylococcus aureus carriage in the developing infant nasal microbiome. Genome Biol. 2020;21(1):301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Rosas-Salazar C, Shilts MH, Tovchigrechko A, et al. Differences in the Nasopharyngeal Microbiome During Acute Respiratory Tract Infection With Human Rhinovirus and Respiratory Syncytial Virus in Infancy. J Infect Dis. 2016;214(12):1924–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Hurst JH, McCumber AW, Aquino JN, et al. Age-Related Changes in the Nasopharyngeal Microbiome Are Associated With Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infection and Symptoms Among Children, Adolescents, and Young Adults. Clin Infect Dis. 2022;75(1):e928–37.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Sirithanakorn C, Cronan JE. Biotin, a universal and essential cofactor: synthesis, ligation and regulation. FEMS Microbiol Rev. 2021;45(4):fuab003.

  76. Zhang J, Wu Y, Liu J, et al. Differential Oral Microbial Input Determines Two Microbiota Pneumo-Types Associated with Health Status. Adv Sci. 2022;9(32):2203115.

  77. Wood DE, Salzberg SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014;15(3):R46.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Gao Y, Li H. Quantifying and comparing bacterial growth dynamics in multiple metagenomic samples. Nat Methods. 2018;15(12):1041–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Tsang TK, Lee KH, Foxman B, et al. Association Between the Respiratory Microbiome and Susceptibility to Influenza Virus Infection. Clin Infect Dis. 2020;71(5):1195–203.

    Article  CAS  PubMed  Google Scholar 

  80. Salk HM, Simon WL, Lambert ND, et al. Taxa of the Nasal Microbiome Are Associated with Influenza-Specific IgA Response to Live Attenuated Influenza Vaccine. PLoS ONE. 2016;11(9):e0162803.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We thank Zhipan Ke of the Center for Disease Control and Prevention of Nanhai District, Foshan, China, for the assistance with sample collection. We thank the volunteers who participated in this study for providing samples.


This project was supported by the National Natural Science Foundation of China (grant number 81971906 to T.D.) and the Guangdong Key Area Research and Development Program (grant number 2019B020228001 to T.D.).

Author information

Authors and Affiliations



H.L. and T.D. created the study concept and designed the experiments. H.L., X.W., and H.Z. performed the sample collection and metadata compilation. H.L., B.C. performed the experiments and analyzed the data. H.L., Y.C., J.Z., R.W. generated research materials. H.L. and T.D. generated figures and legends from analyzed data and interpreted the experiments. H.L. and T.D. wrote and revised the manuscript. T.D. conceived, designed, and supervised the overall study and procured funding.

Corresponding author

Correspondence to Tao Ding.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of Zhongshan School of Medicine, Sun Yat-sen University (ZSSOM-SYSU IRB protocol 2022–054) and was conducted under strict regulations regarding ethics in science and research. Written informed consent was obtained from all participants.

Consent for publication

All authors consent to the publication of this manuscript and declare no conflicts of interest.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

List of primers used for multiplex qPCR. Table S2. The differential genera between virus-positive and negative groups. Table S3. In the random forest model with 100 resamplings, the Mean Decrease Accuracy from group 1 and group 2 was used for the Wilcoxon test. Table S4. The Spearman interaction of differential pathways, differential ARGs, ARGs count, and differential species in the FluA group and healthy controls. Table S5. The growth rate of potential pathogens in situ calculated by GRID algorithm. Table S6. The Spearman interaction of all pathways, ARGs, and pathogens.

Additional file 2: Figure S1.

Multiple qPCR standard curves for FluA, FluB, RSV, and HR. Figure S2. Reads count and proportions of microbiota in sequencing profiles. (A) Microbial reads count in 16S rRNA gene sequencing data. (B) The proportions of microbiota in 16S rRNA gene sequencing data. (C) Microbial and human reads count in metagenomic sequencing data. Figure S3. The variation of microbial composition between virus-positive and negative groups. PCoA based on Bray–Curtis distances was analyzed, and PERMANOVA was used to evaluate the variation. Figure S4. The differential genera, pathway, and ARGs between health and FluA, FluB or FluAB group. (A) Compared to healthy controls, the functional pathways (MetaCyc database) with significant differences in virus infection were calculated. The number and color represent the fold change relative to the control group. Blank blocks indicate no significant difference. (B) The number of antibiotic resistance genes (ARGs) based on metagenomic data was used to evaluate the diversity of ARGs, and the variation between the infected and healthy groups was counted by Kruskal–Wallis test and Wilcoxon test with subsequent Bonferroni correction. (C) Significantly different resistance genes compared to healthy controls were shown. (D) LEfSe was used to calculate the different genera between healthand FluA, FluB, orFluAB participants, respectively. Figure S5. The relative abundance difference of potential pathogens between infected and healthy subjects. (A) Nine pathogenic species (in the blue circle) were identified in the oropharyngeal samples based on metagenomics sequencing analysis. (B) Kruskal-Wallis test and Wilcoxon test was used to calculate the relative abundance difference of nine pathogens between the infected and the healthy subjects. Figure S6. The growth rate of potential pathogens (in the blue circle) was detected based on metagenomics sequencing analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, H., Wu, X., Zeng, H. et al. Unique microbial landscape in the human oropharynx during different types of acute respiratory tract infections. Microbiome 11, 157 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: