Unique microbial landscape in the human oropharynx during different types of acute respiratory tract infections
Microbiome volume 11, Article number: 157 (2023)
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.
Respiratory viruses, such as influenza (Flu) and human respiratory syncytial virus (RSV), are major threats to global public health . 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 . Microbes form a complex system (i.e., the microbiome) closely involved in human health and various conditions . Viral infections disrupt the microbial-host balance maintained by this complex system . 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 . 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 . 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% , 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 , 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).
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 . 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 , 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).
16S rRNA gene analysis
Raw “fastq” files were demultiplexed and processed using tools available in QIIME 2 (v2019.7)  with a Divisive Amplicon Denoising Algorithm (DADA) 2-based pipeline . 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) . 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.
Sequencing raw data pre-processing was performed by KneadData (v0.7.4, http://huttenhower.sph.harvard.edu/kneaddata), which integrates the tools FastQC (v0.11.9) , Trimmomatic (v0.39) , and Bowtie2 (v188.8.131.52) , to do quality check, quality filtering, and host sequences decontamination, respectively. Taxonomic classification of the resulting reads was performed using MetaPhlAn2 (v2.0) . Metabolic pathway abundance was determined using the HUMAnN2 (v0.9.9)  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)  from the metagenomic database using the ShortBRED (v0.9.4)  pipeline. The Growth Rate InDex (GRiD) algorithm (v1.3)  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.
Shannon diversity index was used to characterize alpha-diversity using phyloseq (v1.40.0)  package in R (v4.0.5) . 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) , 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)  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)  and outputted for further analysis, together with node and edge tables for the whole association network.
The distance matrix between microbiota communities was tested by permutational multivariate analysis of variance (PERMANOVA) using vegan (v2.6.2)  package in R (v4.0.5). Linear discriminant analysis (LDA) effect size (LEfSe) was conducted on Galaxy (v1.0)  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)  package. Spearman correlation analysis was performed using dplyr (v1.1.2)  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)  and R stats (v4.2.2)  package, and ggplot2 (v3.4.1)  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.
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).
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.
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.
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).
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 . 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 . 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 , 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 . 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 . 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  and Kraken  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 . Among healthy young adults vaccinated with live attenuated influenza vaccine, S. infantis was positively associated with influenza H1 immunoglobulin A (IgA) titers . 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 http://www.ncbi.nlm.nih.gov/sra. The R scripts used to generate the results in this manuscript are available in the following Git repository: https://github.com/TaoDing/URTI_microbiota.
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.
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.
Morris DE, Cleary DW, Clarke SC. Secondary bacterial infections associated with influenza pandemics. Front Microbiol. 2017;8:1041.
Liu Y, Ling L, Wong SH, et al. Outcomes of respiratory viral-bacterial coinfection in adult hospitalized patients. EClinicalMedicine. 2021;37:100955.
Almand EA, Moore MD, Jaykus L-A. Virus-bacteria interactions: an emerging topic in human infection. Viruses. 2017;9(3):58.
Finlay BB, McFadden G. Anti-immunology: evasion of the host immune system by bacterial and viral pathogens. Cell. 2006;124(4):767–82.
Consortium HMP. A framework for human microbiome research. Nature. 2012;486(7402):215–21.
Wypych TP, Wickramasinghe LC, Marsland BJ. The influence of the microbiome on respiratory health. Nat Immunol. 2019;20(10):1279–90.
Hanada S, Pirzadeh M, Carver KY, et al. Respiratory viral infection-induced microbiome alterations and secondary bacterial pneumonia. Front Immunol. 2018;9:2640.
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.
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.
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.
Kilbourne ED. Influenza pandemics of the 20th century. Emerging Infect Dis. 2006;12(1):9.
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.
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.
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.
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.
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.
Chertow DS, Memoli MJ. Bacterial coinfection in influenza: a grand rounds review. JAMA. 2013;309(3):275–82.
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.
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.
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.
Blaschke AJ, Korgenski EK, Wilkes J, et al. Rhinovirus in febrile infants and risk of bacterial infection. Pediatrics. 2018;141(2).
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.
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.
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.
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.
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.
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.
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.
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.
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.
Callahan BJ, McMurdie PJ, Rosen MJ, et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.
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.
Andrews S. FastQC a quality control tool for high throughput sequence dat. 2016. http://www.bioinformatics.babraham.ac.uk/projects/fastqc.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923.
Truong DT, Franzosa EA, Tickle TL, et al. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat Methods. 2015;12(10):902–3.
Franzosa EA, McIver LJ, Rahnavard G, et al. Species-level functional profiling of metagenomes and metatranscriptomes. Nat Methods. 2018;15(11):962–8.
Liu B, Pop M. ARDB–Antibiotic Resistance Genes Database. Nucleic Acids Res. 2009;37:D443–7.
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.
Emiola A, Oh J. High throughput in situ metagenomic measurement of bacterial replication at ultra-low sequencing coverage. Nat Commun. 2018;9(1):4956.
Smith AM, McCullers JA. Secondary bacterial infections in influenza virus infection pathogenesis. Curr Top Microbiol Immunol. 2014;385:327–56.
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.
Li ZJ, Zhang HY, Ren LL, et al. Etiological and epidemiological features of acute respiratory infections in China. Nat Commun. 2021;12(1):5026.
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.
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.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217. https://doi.org/10.1371/journal.pone.0061217.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018. https://www.R-project.org/.
Liaw A, Wiener M. “Classification and Regression by randomForest. R News. 2002;2(3):18–22 https://CRAN.R-project.org/doc/Rnews/.
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. https://doi.org/10.1186/1471-2105-12-77.
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.
Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14(6):927–30.
Segata N, Izard J, Waldron L, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
Vincent QV. ggbiplot: a ggplot2 based biplot. R package version 0.55. 2011. https://github.com/vqv/ggbiplot/tree/experimental.
Wickham H, François R, Henry L, et al. dplyr: A Grammar of Data Manipulation. R package version 1.1.2. 2023. https://CRAN.R-project.org/package=dplyr.
Kassambara A. ggpubr: 'ggplot2' Based Publication Ready Plots. R package version 0.6.0. 2023. https://rpkgs.datanovia.com/ggpubr/.
Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer-Verlag; 2016. ISBN 978-3-319-24277-4. https://ggplot2.tidyverse.org.
Khan R, Petersen FC, Shekhar S. Commensal Bacteria: An Emerging Player in Defense Against Respiratory Pathogens. Front Immunol. 2019;10:1203–1203.
Keith JW, Pamer EG. Enlisting commensal microbes to resist antibiotic-resistant pathogens. J Exp Med. 2019;216(1):10–9.
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.
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.
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.
Brundage JF, Shanks GD. Deaths from bacterial pneumonia during 1918–19 influenza pandemic. Emerging Infect Dis. 2008;14(8):1193.
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.
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.
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.
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.
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.
Dereschuk K, Apostol L, Ranjan I, et al. Identification of Lung and Blood Microbiota Implicated in COVID-19 Prognosis. Cells. 2021;10(6).
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.
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.
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. https://doi.org/10.1093/infdis/jiw456.
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. https://doi.org/10.1093/cid/ciac184.
Sirithanakorn C, Cronan JE. Biotin, a universal and essential cofactor: synthesis, ligation and regulation. FEMS Microbiol Rev. 2021;45(4):fuab003.
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.
Wood DE, Salzberg SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014;15(3):R46.
Gao Y, Li H. Quantifying and comparing bacterial growth dynamics in multiple metagenomic samples. Nat Methods. 2018;15(12):1041–4.
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.
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.
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.).
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.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
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.
About this article
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). https://doi.org/10.1186/s40168-023-01597-9