Nasopharyngeal Microbiota Profiling of SARS-CoV-2 Infected Patients
Biological Procedures Online volume 22, Article number: 18 (2020)
We analyzed the bacterial communities of the nasopharynx in 40 SARS-CoV-2 infected and uninfected patients. All infected patients had a mild COVID-19 disease. We did not find statistically significant differences in either bacterial richness and diversity or composition. These findings suggest a nasopharyngeal microbiota at least early resilient to SARS-CoV-2 infection.
The human upper respiratory tract is the major portal of entry for infectious droplet or aerosol-transmitted microorganisms, including the 2019 emerged severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) . Consistent with an ongoing rise in the number of coronavirus disease 2019 (COVID-19) cases worldwide, the updated World Health Organization estimates on 29 April 2020 reported 3,023,788 confirmed cases, including 208,112 deaths, over time (https://covid19.who.int/). Because of apparently absent cross-protective immunity from related viral infections, SARS-CoV-2 transmissibility is high, hence facilitating widespread person-to-person transmission . Both asymptomatic and symptomatic patients may transmit the virus , consistent with similar viral loads detected in the nasal (sampled from the nasopharynx) and throat swabs from both patients .
As with influenza virus, respiratory samples provide the greatest yield of SARS-CoV-2 nucleic acid shedding [4, 5]. In one study analyzing nine patients with mild courses of COVID-19, swab samples taken during the first week of symptoms displayed high loads and successful isolation of the virus and, importantly, the presence of viral replicative RNA intermediates confirmed the active replication of SARS-CoV-2 in the throat . The virus uses the host receptor angiotensin-converting enzyme 2 (ACE2), predominantly expressed in the lung, to interact with its spike protein . This interaction would be responsible of the extension of tropism to multiple tissues , which include the throat tissue with anyhow low ACE2 expression , following a fusion activity gain of SARS-CoV-2 spike protein at the S1-S2 junction that is not present in SARS-CoV .
Unlike influenza virus [9,10,11,12], no studies have explored the relationship between SARS-CoV-2 infection and the bacterial community (also referred to as microbiota) within the nasopharynx. Therefore, we investigated whether the presence of the virus in the nasopharynx might reflect alterations of the resident microbiota by comparing the bacterial communities of SARS-CoV-2 infected and uninfected patients.
We collected nasopharynx samples from patients who were (n = 18) or were not (n = 22) diagnosed with COVID-19 based on the SARS-CoV-2 RNA detection using a real-time reverse transcriptase-polymerase chain reaction (RT-PCR) assay . This study received appropriate ethical review committee approval, with a waiver of informed consent.
Patients with samples collected based on the clinical COVID-19 suspicion (i.e., who presented symptoms of acute respiratory infection in the emergency department of the Fondazione Policlinico Universitario A. Gemelli (FPG) IRCCS Hospital of Rome, Rome, Italy)  from March 6 through March 9, 2020, were included. Nasopharyngeal swabs were tested for RT-PCR based COVID-19 diagnosis  using the Korean Ministry of Food and Drug Safety approved Allplex 2019-nCoV assay (Arrow Diagnostics S.r.l., Genova, Italy). The same samples were used to profile the nasopharyngeal microbiota in all 40 patients concomitantly.
Genomic DNA was extracted using the QIAGEN EZ1 Advanced XL system and the V5–V6 hypervariable region of the bacterial 16S rRNA gene was sequenced using an Illumina MiSeq V2 chemistry (2 × 250 bp) as described . Sequencing data were processed for alignment and quality filtering in QIIME2 v2019.1 [15, 16], and representative amplicon sequence variants were obtained by the DADA2 algorithm  available at https://doi.org/10.1038/s41587-019-0209-9. Samples were rarefied to the minimum number of sequence reads obtained totally (range, 107,933–165,427) to perform subsequent analyses. Taxonomic annotation was performed using both VSEARCH  and SILVA database v132 . Statistical tests, including Wilcoxon signed rank test and permutational analysis of variance (PERMANOVA) were performed in R software (R Foundation for Statistical Computing, version 3.6.0).
No patients developed pneumonia as documented in the absence of significant abnormalities (i.e., multifocal ground-glass opacities) on chest computed tomography (CT) images.
The microbiota of the nasopharynx was not different in patients positive for SARS-CoV-2 RNA compared to the microbiota of patients negative for SARS-CoV-2 RNA (Fig. 1). No significant differences between the patient groups in either bacterial richness and diversity (observed species, Shannon index, and inverse Simpson index were assessed; P > 0.05 by Wilcoxon signed rank test) or composition (no samples clustering within each group, as visualized by multidimensional scaling; P > 0.05 by PERMANOVA) was noticed.
We analyzed the relative abundance of bacterial taxa most represented in both patient groups and did not observe significant differences between the groups. Most sequences in all samples (98% in both SARS-CoV-2 infected and uninfected patients) belonged to five phyla, namely Firmicutes (42 and 51%, respectively), Bacteroidetes (25 and 20%, respectively), Proteobacteria (18 and 16%, respectively), Actinobacteria (8 and 6%, respectively), and Fusobacteria (5 and 5%, respectively). These findings were consistent with the high similarity shared by the bacterial communities in both SARS-CoV-2 infected and uninfected patients.
In this study, SARS-CoV-2 infection was not associated with a different profile of the nasopharynx bacterial community analyzed at COVID-19 diagnosis, which is when the patients had a mild disease. This suggests that no microbiota compositional alterations occurred in response to SARS-CoV-2 or that SARS-CoV-2 was unable to induce these alterations in our patients. Both suggestions lend support to a nasopharyngeal microbiota at least early resilient to SARS-CoV-2 infection.
As nasopharyngeal swabs can be a reasonable proxy for lung samples , resilience of the nasopharynx bacterial community would also imply robustness of lung bacterial community in restricting SARS-CoV-2 growth or attachment. This would mirror a lack of effect on the oropharyngeal microbiota in volunteers challenged with H3N2 influenza virus, many of which developed very mild disease . Thus, we did not find bacterial taxa for which relative abundance was significantly higher in SARS-CoV-2 infected patients. Conversely, taxa such as Dolosigranulum, Moraxella, Staphylococcus, and Streptococcus (three of which also detected in our study) were found to be enriched in the nasopharynx microbiota of H3N2 influenza virus infected patients compared to healthy control individuals .
In one study that analyzed elderly pneumonia patients in comparison with healthy elderly , oropharyngeal microbiota profiles differed significantly between the groups, with three of them being significantly associated with pneumonia. Therefore, it is surprising that no published studies until now reported association between upper respiratory tract microbiota and COVID-19, especially in severe disease cases. Here, we studied COVID-19 patients at a time when (mild) symptoms were typical of upper respiratory tract infection and were not affecting the lungs. Consistently, no patients had multifocal ground-glass opacities on chest CT , as mentioned above.
The limitations of this study are that SARS-CoV-2 infections were defined by RT-PCR positivity, with no culture results, likely leading to include infections with low infectiousness and the number of analyzed samples was small. Although RT-PCR performed on a nasopharyngeal swab remains the reference diagnostic test , patients with negative RT-PCR results could be highly probable COVID-19 cases based on positive chest CT findings . We did not use chest CT to screen for SARS-CoV-2 infection  because of awareness that CT findings can be absent, especially in patients with early and/or mild disease . Therefore, we cannot exclude the possibility that false negative RT-PCR results have biased the grouping of patients in our study. Further investigation of SARS-CoV-2 infected patients with detailed viral shedding data and sequential respiratory samples is warranted.
Availability of Data and Materials
All data analyzed in this study are available from the corresponding author on request.
Zhu N, Zhang D, Wang W, et al. A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med. 2020;382:727–33.
Chan JF, Yuan S, Kok KH, et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet. 2020;395:514–23.
Li R, Pei S, Chen B, et al. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). Science. 2020;368:489–93.
Zou L, Ruan F, Huang M, et al. SARS-CoV-2 viral load in upper respiratory specimens of infected patients. N Engl J Med. 2020;382:1177–9.
Wölfel R, Corman VM, Guggemos W, et al. Virological assessment of hospitalized patients with COVID-2019. Nature. 2020;581:465–9.
Hoffmann M, Kleine-Weber H, Schroeder S, et al. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 2020; 181:271–280.e8.
Sungnak W, Huang N, Bécavin C, et al. SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes. Nat Med. 2020;26:681–7.
Letko M, Marzi A, Munster V. Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses. Nat Microbiol. 2020;5:562–9.
Greninger AL, Chen EC, Sittler T, et al. A metagenomic analysis of pandemic influenza a (2009 H1N1) infection in patients from North America. PLoS One. 2010;5:e13381.
Hanada S, Pirzadeh M, Carver KY, Deng JC. Respiratory viral infection-induced microbiome alterations and secondary bacterial pneumonia. Front Immunol. 2018;9:2640.
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:e01296–19.
Ramos-Sevillano E, Wade WG, Mann A, et al. The effect of influenza virus on the human oropharyngeal microbiome. Clin Infect Dis. 2019;68:1993–2002.
WHO. Laboratory testing for coronavirus disease (COVID-19) in suspected human cases. Interim guidance. Geneva, World Health Organization, 2020. Available at: https://apps.who.int/iris/bitstream/handle/10665/331501/WHO-COVID-19-laboratory-2020.5-eng.pdf?sequence=1&isAllowed=y. Accessed 26 April 2020.
WHO. Global surveillance for COVID-19 caused by human infection with COVID-19 virus. Interim guidance. Geneva, World Health Organization, 2020. Available at: https://www.who.int/docs/default-source/coronaviruse/global-surveillance-for-covid-v-19-final200321-rev.pdf. Accessed 26 April 2020.
Mohsen A, Park J, Chen YA, Kawashima H, Mizuguchi K. Impact of quality trimming on the efficiency of reads joining and diversity analysis of Illumina paired-end reads in the context of QIIME1 and QIIME2 microbiome analysis frameworks. BMC Bioinformatics. 2019;20:581.
Hall M, Beiko RG. 16S rRNA gene analysis with QIIME2. Methods Mol Biol. 1849;2018:113–29.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584.
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. 2013;41:D590–6.
Charlson ES, Bittinger K, Haas AR, et al. Topographical continuity of bacterial populations in the healthy human respiratory tract. Am J Respir Crit Care Med. 2011;184:957–63.
El Homsi M, Chung M, Bernheim A, et al. Review of chest CT manifestations of COVID-19 infection. Eur J Radiol Open. 2020;7:100239.
Ai T, Yang Z, Hou H, et al. Correlation of chest CT and RT-PCR testing in coronavirus disease 2019 (COVID-19) in China: a report of 1014 cases. Radiology. 2020;200642.
Long C, Xu H, Shen Q, et al. Diagnosis of the coronavirus disease (COVID-19): rRT-PCR or CT? Eur J Radiol. 2020;126:108961.
We thank Franziska Lohmeyer, PhD (Scientific Direction, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy) for English revision of the manuscript, and clinical staff members in Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy for collection of samples. None of these individuals received compensation for their contributions.
This work was funded by donations from Reale Group and Fondazione Valentino Garavani & Giancarlo Giammetti to support the COVID-19 Research in our Institution. The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Ethics Approval and Consent to Participate
This study was approved by the Ethics Committee of the FPG (approval number 17057/20), with a waiver of informed consent.
Consent for Publication
The authors declared that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
De Maio, F., Posteraro, B., Ponziani, F.R. et al. Nasopharyngeal Microbiota Profiling of SARS-CoV-2 Infected Patients. Biol Proced Online 22, 18 (2020). https://doi.org/10.1186/s12575-020-00131-7